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HIGH THROUGHPUT DIRECTED EVOLUTION BY RATIONAL 
. MUTAGENESIS 

RELATED APPLICATIONS 

Benefit of priority is claimed to U.S. application Serial No. 
10/022,249, filed December 17, 2001 , to Manuel Vega and Lila Drittanti, 
entitled "HIGH THROUGHPUT DIRECTED EVOLUTION BY RATIONAL 

5 MUTAGENESIS." Benefit of priority is also claimed to U.S. provisional 
patent application serial No. 60/315,382, filed August 27, 2001, to 
Manuel Vega, Lila Drittanti and Marjorie Flaux, entitled "HIGH THROUGH- 
PUT DIRECTED EVOLUTION BY RATIONAL MUTAGENESIS." Where 
permitted, the subject matter of each of these applications is incorporated 
10 in its entirety by reference thereto. 
FIELD OF INVENTION 

Processes and systems for the high throughput directed evolution 
of peptides and proteins, particularly those that act in complex biological 
settings, are provided. The proteins and peptides include, but are not 
15 limited to, intracellular proteins, messenger/signaling/hormone proteins 
and viral proteins. 
BACKGROUND 

Directed evolution refers to biotechnological processes for 
optimizing the activity of proteins by means of random changes 
20 introduced into selected respective genes. Directed evolution involves the 
creation of a library of mutated genes, and then selection of the mutants 
that encode proteins having desired properties. The process can be an 
iterative one in which gene products that have improvement in a desired 
property are subjected to further cycles of mutation and screening. 
25 Directed evolution provides a way to adapt natural proteins to work in 
new chemical or biological environments, and/or to elicit new functions. 
The .potential plasticity of proteins is such that chances exist that for 
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every new challenge, such as a new environment and desired new or 
altered activity, it should be possible, given a sufficient pool of modified 
proteins {or encoding nucleic acids), that an appropriately 'evolved' 
protein could be found that would have a desired activity. The problem is 
5 in generating and then identifying the appropriate sequence. 

There have been practical approaches to this problem (see, e.g., 
U.S. Patent Nos. 6,096,548; 6,117,679; 6,165,793; 6,180,406; 
6,132,970; 6,171,820; 6,238,884; 6,174,673; 6,057,103; 6,001,574; 
5,763,239; 5,837,500; 5,571,698; 6,156,509; 5,723,323; 5,862,514; 

10 5,871,974; 5,779,434 and others). Typically theses approaches are of 
two types. One is a purely "rational" approach that is based on the 
assumption that the optimized proteins can be rationally designed. This, 
however, requires sufficient information regarding the laws that govern 
protein folding, molecular interactions, intra-molecular forces and other 

15 dynamics of protein activity. This rational approach is extremely 

dependent on a number of variables and parameters that are not known. 
Consequently, although useful in some specific cases and applications, 
the rational approach intended to 'predict' protein structure remains 
limited in applicability. 

20 In contrast to the rational approach, random approaches have also 

been employed. One random approach requires synthesis of all possible 
protein sequences or a statistically sufficient large number of proteins and 
then screening them to identify proteins having the desired activity or 
property. Since the resources to synthesize all possible theoretical 

25 sequences of a single protein is not possible, this approach is 

impracticable. Other random approaches are based on gene shuffling 
methods, which are PCR-based methods that generate random 
rearrangements between two or more sequence-related genes to 
randomly generate variants of the gene. 
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The development and scope of directed evolution, thus, has been 
limited, and its potential remains to be exploited. In order to exploit the 
potential of directed evolution, alternative approaches for generating and 
identifying evolved proteins are needed. It is an object herein to provide 
5 methods and products to exploit the potential of directed evolution. 
SUMMARY 

Provided herein are methods for performing directed evolution for 
the optimization of proteins that function in complex biological settings. 
Methods of high throughput directed evolution of proteins are provided. 
10 In practicing the methods, each molecule is individually designed, 

produced, processed, screened and tested in a high throughput format. 
Neither random or combinatorial methods nor mixtures of molecules are 
used. 

The methods provided herein include the steps of identifying a 
15 protein target of interest; obtaining nucleic acids that encode the target, 
which may be from any source, such as a natural library, a collection 
generated by known gene shuffling techniques and related methods, and, 
then creating variants of the proteins using methods for rational 
mutagenesis provided herein. Whatever method is used to select or 
20 generate the nucleic acids encoding the protein targets, each molecule is 
processed and screened separately in a high throughput format. 

The nucleic acids encoding each variant are individually screened. 
They can be screened in any suitable assay, including cell-based assays 
and biochemical assays. For cell based assays, each nucleic acid 
25 molecule is introduced into an expression vector for expression in a 

bacterial cell or into a vector for expression in a eukaryotic host cell. In all 
instances, the nucleic acids of interest are introduced into host cells in an 
expression vector, such as by transfection for bacterial hosts and 
transduction with viral vectors into eukaryotic hosts with viral vectors. 
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Each variant is introduced into a host and the resulting cells are 
maintained separately, such as in an addressable array of wells in a 
microtiter plate or other substrate with discrete locations for performing 
reactions or retaining molecules of interest. Typical formats are 96 loci, 
5 and multiples thereof (384, 1536, 3072, . . . 96 x n, where n is 1 to any 
number desired, such as 10, 20, 30, 50 . . .100), although any 
convenient number of loci may be employed. 

Since the process is conducted in a high throughput format, for 
many embodiments, it is often important to assess the relative numbers 

10 of transformed, transduced or transfected cells. Hence the relative (or 
actual) titer of the vector, such as the recombinant viral vector, must be 
known to permit analysis of results. For high throughput formats, it is 
important to assess the relative or actual concentration of the viral vector 
(or plasmid) so that results can be compared among ail cells and variants. 

15 Methods for titering (determining the concentration) of the nucleic acid 
encoding the variant and/or the recombinant virus are also provided. 

The processes require accurate titering of the viruses in a collection 
or among collections (libraries) so that the activities of the screened 
mutant proteins can be compared. Provided are general methods for the 

20 quantitative assessment of the parameters of activity corresponding to 

the individual variants in the library, based upon intracellular serial dilution 
generated by precise titering with the gene transfer viral vectors. Any 
method that permits accurate titering may be used, including that 
described in International PCT application No. PCT/FR01/01 366, based on 

25 French application n° 0005852, filed 9 May 2000, and published as 

International PCT application No. WO 01/186291 . A method of titering, 
designated Tagged Replication and Expression Enhancement Technology 
(TREE") is provided herein. 
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Each of the different cells is separately screened by a suitable 
assay, and the results analyzed. Methods for assessing the interactions 
in biological systems, such as a Hill-based analysis (see, published 
International PCT application No. WO 01/44809 based on PCT n° 
5 PCT/FR00/O35O3, Dec, 2000, and the description herein), or a second 
order polynomial or other algorithm that describes the interaction between 
cells and biological agents to select variants that have a desired property 
are employed in the processes herein. 

A semi-rational method for evolution of proteins that is particularly 

10 designed for use in the methods herein or in any method that uses 

"evolved" proteins is also provided. The method, which is based on an 
amino-acid scanning protocol, is for rationally designing the variants for 
use in the directed evolution and selection method, and can employ 
iterative processing of the steps of the high throughput methods provided 

1 5 herein. In this method, once the target protein or domain is identified, 
nucleic acid molecules encoding variants are prepared. Each variant 
encoded by the nucleic acid molecules has a single amino acid replaced 
with another selected amino acid, such as alanine (Ala), glycine (Qly), 
serine (Ser) or any other suitable amino acid, typically one selected to 

20 have a neutral effect on secondary and tertiary structure. The resulting 
series of variants are separately screened in the high throughput format 
provided herein, and those that have a change in the target activity are 
selected and the modified amino acids are designated "hits." Nucleic acid 
molecules encoding proteins in which each hit position is replaced by the 

25 eighteen remaining amino acids then are synthesized and the resulting 
collection of molecules is screened, such by introduction into host cells, 
and the proteins that result in an improvement of a targeted activity, are 
identified. Such proteins are designated "leads." Leads may be further 
modified by producing proteins that have combinations of the mutations 
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identified in the leads. This method, which does not require any 
knowledge of the structure of a target protein, permits precise control of 
locations where changes are introduced and also the amount of change 
that is introduced. . .... 

5 The high throughput directed evolution processes provided herein 

include the use of virus libraries containing mutant versions of a gene; 
viral libraries of such mutant genes are also provided. 

Reporter cells are infected with the titered viruses that encode the 
mutant genes. The mutant genes are expressed and read-out data from 

10 either biochemical or cell-based assays, while isolating each mutant/virus 
physically from the others (i.e. one-by-one analysis), are collected and 
analyzed. Serial dilution assays (i.e. a series of dilutions for each 
individual mutant/virus in the library) are used and the biochemical/cell- 
based assays are performed on each single dilution for each individual 

15 mutant/virus. Analysis of the serial dilution readout-data can be 
performed using any method of analysis that permits one-by-one 
comparisons. Hill-based analysis (see, published International PCT 
application No. WO 01/44809 based on PCT n° PCT/FROO/03503, Dec, 
2000, and the description herein) are employed for analysis of the data. 

20 Protein/protein domain variants identified using the methods are 

also provided. Also provided are nucleic acid molecules and proteins and 
polypeptides produced by the methods and viruses and cells that contain 
the nucleic acid molecules and proteins. 

In an exemplary embodiment of methods provided herein, the 

25 process of rational directed evolution provided herein is applied to the 
AAV rep gene. The resulting recombinant rep protein variants and rAAV 
are also provided. Among the rep proteins are those that result in 
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increased rAAV production in rAAV that encode such mutants, thereby, 
among a variety of advantages, offering a solution to the need in the gene 
therapy industry to increase the production therapeutic vectors without 
up-scaling manufacturing. 
5 Thus, for exemplification, some methods provided herein have been 

used to identify amino acid "hit" positions in adeno-associated virus 
(AAV) rep proteins that are relevant for AAV or rAAV production. Those 
amino acid positions are such that a change in the amino acid leads to a 
change in protein activity either to lower activity or to higher activity 
10 compared to native-sequence Rep proteins. The hit positions were then 
used to generate further mutants designated 'Meads." Provided herein are 
the resulting mutant rep proteins that result in either higher or lower 
levels of AAV or rAAV virus compared to the wild-type (native) Rep 
protein (s). 

15 In addition to enhancing AAV production, among the rep mutants 

are those that inhibit papillomavirus <PV) and PV-associated diseases, 
including certain cancers and human immunodeficiency virus (HIV) and 
HIV-associated diseases. 

Systems and computer controlled systems for performing the high 

20 throughput processes are also provided. 
DESCRIPTION OF THE FIGURES 

FIGURES 1 A - 1E summarize various exemplary embodiments of 
the high throughput processes provided herein. FIGURE 1 A depicts an 
embodiment of the process in which an amino acid scan is employed to 

25 generate a library of mutants, which are then introduced into viral 

vectors, such as an adeno-associated viral vector (AAV), a herpes virus, 
such as herpes simplex virus (HSV) and other herpes virus vectors, a 
vaccinia virus vector, retroviral vectors, such as MuMLV, MoMLV, feline 
leukemia virus, and HIV and other lentiviruses, adenovirus vectors and 
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other suitable viral vector, each member of the library is individually 
tested and phenotypically characterized to identify HITS. FIGURE 1 B 
summarizes round 2 in which LEADS are developed by mutagenesis at 
and/or surrounding the positions identified as HITS; FIGURE 1C 
5 summarizes the optional next round in which recombination among 
LEADS is performed to further optimize the LEADS; FIGURE 1D depicts 
the process in mammalian cells; and FIGURE 1E depicts the process in 
bacterial cells. 

FIGURE 2A depicts an exemplary titering process (in this instance 
10 the TREE" for titering AAV) in a 96 well format; FIGURE 2B shows the 
results and analysis of a titering process performed using the TREE" 
procedure; and FIGURE 2C shows an exemplary calibration curve for the 
calculation of the titer using the TREE" method. 

FIGURE 3A and 3B depict "HITS" and "LEADS" respectively for . 
15 identification of AAV rep mutants "evolved" for increased activity. 

FIGURE 4 shows the genetic map of AAV, including the location of 
promoters, and transcripts; amino acid 1 of the Rep 78 gene is at 
nucleotide 321 in the AAV-2 genome. 

FIGURES 5A and 5B show the alignment of amino acid sequences 
20 of Rep78 among AAV-1; AAV-6; AAV-3; AAV-3B; AAV-4; AAV-2; AAV- 
5 sequences, respectively; the hit positions with 100 percent homology 
among the serotypes are bolded italics, where the position is different 
(compared to AAV-2, no. 6 in the Figure) in a particular serotype, it is in 
bold; a sequence indicating relative conservation of sequences among 
25 the serotypes is labeled "C". 
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Legend: 

1 is AAV-1; 2 is AAV-6, 3 is AAV-3, 4 is AAV-3B, 
5 is AAV-4, 6 is AAV-2, and 7 is AAV-5; 
V* where the amino acid is present > 20%; 
5 ":" where the amino acid is present ^ 40%; 

" + M where the amino acid is present ^ 60%; 
"*" where the amino acid is present ^ 80%; and 
where the amino acid is the same amongst all 
serotypes depicted it is represented by its single letter 
-10 code. 
DETAILED DESCRIPTION 
A. Definitions 

Unless defined otherwise, all technical and scientific terms used 
herein have the same meaning as is commonly understood by one of skill 
15 in the art to which this invention belongs. All patents, patent 

applications, published applications and publications, Genbank sequences, 
websites and other published materials referred to throughout the entire 
disclosure herein are, unless noted otherwise, incorporated by reference 
in their entirety. In the event that there are a plurality of definitions for 
20 terms herein, those in this section prevail. 

As used herein, directed evolution refers to methods that adapt 
natural proteins or protein domains to work in new chemical or biological 
environments and/or to elicit new functions. It is more a more broad- 
based technology than DNA shuffling. 
25 As used herein, high-throughput screening (HTS) refers to 

processes that test a large number of samples, such as samples of test 
proteins or cells containing nucleic acids encoding the proteins of interest 
to identify structures of interest or to identify test compounds that 
interact with the variant proteins or cells containing them. HTS 
30 operations are amenable to automation and are typically computerized to 
handle sample preparation, assay procedures and the subsequent 
processing of large volumes of data. 
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As used herein, DNA shuffling is a PCR-based technology that 
produces random rearrangements between two or more sequence-related 
genes to generate related, although different, variants of given gene. 

As used herein, "hits" are mutant proteins that have an alteration in 
5 any attribute, chemical, physical or biological property in which such 
alteration is sought. In the methods herein, hits are generally generated 
by systematically replacing each amino acid in a protein or a domain 
thereof with a selected amino acid, typically Alanine, Glycine, Serine or 
any amino acid, as long as each residue is replaced with the same 

10 residue. Hits may be generated by other methods known to those of skill 
in the art and tested by the high throughput methods herein. For 
purposes herein a Hit typically has activity with respect to the function of 
interest that differs by at least 10%, 20%, 30% or more from the wild 
type or native protein. The desired alteration, which is generally a 

15 reduction in activity, will depend upon the function or property of interest. 
As used herein, "leads" are "hits" whose activity has been 
optimized for the particular attribute, chemical, physical or biological 
property. In the methods herein, leads are generally produced by 
systematically replacing the hit loci with all remaining 1 8 amino acids, and 

20 identifying those among the resulting proteins that have a desired activity. 
The leads may be further optimized by replacement of a plurality of "hit" 
residues. Leads may be generated by other methods known to those of 
skill in the art and tested by the high throughput methods herein. For 
purposes herein a lead typically has activity with respect to the function 

25 of interest that differs from the native activity, by a desired amount and is 
at by at least 10%, 20%, 30% or more from the wild type or native 
protein. Generally a Lead will have an activity that is 2 to 10 or more 
times the native protein for the activity of interest. As with hits, the 
change in the activity is dependent upon the activity that is "evolved." 
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The desired alteration will depend upon the function or property of 
interest. 

As used herein, MOI is multiplicity of infection. 
As used herein, ip, with reference to a virus or recombinant vector, 
5 refers to a titer of infectious particles. 

As used herein, pp refers to the total number of vector (or virus) 
physical particles 

As used herein, biological and pharmacological activity includes any 
activity of a biological pharmaceutical agent and includes, but is not 

10 limited to, biological efficiency, transduction efficiency, gene/transgene 
expression, differential gene expression and induction activity, titer, 
progeny productivity, toxicity, cytotoxicity, immunogenicity, cell 
proliferation and/or differentiation activity, anti-viral activity, 
morphogenetic activity, teratogenetic activity, pathogenetic activity, 

15 therapeutic activity, tumor suppressor activity, ontogenetic activity, 
oncogenetic activity, enzymatic activity, pharmacological activity, 
cell/tissue tropism and delivery. 

As used herein, "output signal" refers to parameters that can be 
followed over time and, if desired, quantified. For example, when a virus 

20 infects a cell, the infected cell undergoes a number of changes. Any 
such change that can be monitored and used to assess infection, is an 
"output signal," and the cell is referred to as a "reporter cell." Output 
signals include, but are not limited to, enzyme activity, fluorescence, 
luminescence, amount of product produced and other such signals. 

25 Output signals include expression of a viral gene or viral gene product, 
including heterologous genes (transgenes) inserted into the virus. Such 
expression is a function of time ("t") after infection, which in turn is 
related to the amount of virus used to infect the cell, and, hence, the 
concentration of virus ("s") in the infecting composition. For higher 
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concentrations the output signal is higher. For any particular 
concentration, the output signal increases as a function of time until a 
plateau is reached. Output signals may also measure the interaction 
between cells, expressing heterologous genes, and biological agents 
5 As used herein, adeno-associated virus (AAV) is a defective and 

non-pathogenic parvovirus that requires co-infection with either 
adenovirus or herpes virus for its growth and multiplication, able of 
providing helper functions. A variety of serotypes are known, and 
contemplated herein. Such serotypes include, but are not limited to: 

10 AAV-1 (Genbank accession no. NC002077; accession no. VR-645); AAV- 
2 (Genbank accession no. NC001401; accession no. VR-680); AAV-3 
(Genbank accession no. NC001729; accession no. VR-681); AAV-3b 
(Genbank accession no. NC001863); AAV-4 (Genbank accession no. 
NC001829; ATCC accession no. VR-646 ); AAV-6 (Genbank accession 

15 no.NC001862); and avian associated adeno-virus (ATCC accession no. 
VR-1449). The preparation and use of AAVs as vectors for gene 
expression in vitro and for in vivo use for gene therapy are well known 
(see, e.g., U.S. Patent Nos. 4,797,368, 5,139,941, 5,798,390 and 
6,127,175; Tessier era/. (2001) J. Virol. 75:375-383; Salvetti et ai. 

20 (1 998) Hum Gene Ther 20:695-706; Chadeuf et a/. (2000) J Gene Med 
2:260-268). 

As used herein, the activity of a Rep protein or of a capsid protein 
refers to any biological activity that can be assessed. In particular, 
herein, the activity assessed for the rep proteins is the amount {i.e., titer) 
25 of AAV produced by a cell. 

As used herein, the Hill equation is a mathematical model that 
relates the concentration of a drug (i.e., test compound or substance) to 
the response being measured 

30 
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£naxIQT 

y = [D] n + [D 50 ] n 

where y is the variable being 
5 measured, such as a response, 

signal, y max is the maximal response achievable, [D] is the molar 
concentration of a drug, [D 50 ] is the concentration that produces a 50% 
maximal response to the drug, n is the slope parameter, which is 1 if the 
drug binds to a single site and with no cooperativity between or among 
10 sites. A Hill plot is log 10 of the ratio of ligand-occupied receptor to free 
receptor vs. log [D] (M). The slope is n, where a slope of greater than 1 
indicates cooperativity among binding sites, and a slope of less than 1 
can indicate heterogeneity of binding. This general equation has been 
employed for assessing interactions in complex biological systems (see, 
15 published International PCT application No. WO 01/44809 based on PCT 
n° PCT/FR00/03503, see, also, EXAMPLES). 

As used herein, in the Hill-based analysis (published International 
PCT application No. WO 01/44809 based on PCT n° PCT/FR00/O3503), 
the parameters, n,K,T,e,nfi, are as follows: 
20 /7 potency of the biological agent acting on the assay (cell- 

based) system; 

k constant of resistance of the assay system to elicit a 
response to a biological agent; 

€ is global efficiency of the process or reaction triggered by the 
25 biological agent on the assay system; 

t is the apparent titer of the biological agent; 
0 is the absolute titer of the biological agent; and 
n is the heterogeneity of the biological process or reaction. 
In particular, as used herein, the parameters n (potency) or k 
30 (constant of resistance) are used to respectively assess the potency of a 
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test agent to produce a response in an assay system and the resistance 
of the assay system to respond to the agent. 

As used herein, e(efficiency), is the slope at the inflection point of 
the Hill curve (or, in general, of any other sigmoidal or linear 
5 approximation), to assess the efficiency of the global reaction (the 
biological agent and the assay system taken together) to elicit the 
biological or pharmacological response. 

As used herein, r (apparent titer) is used to measure the limiting 
dilution or the apparent titer of the biological agent. 
10 As used herein, 6 (absolute titer), is used to measure the absolute 

limiting dilution or titer of the biological agent. 

As used herein, n (heterogeneity) measures the existence of 
discontinuous phases along the global reaction, which is reflected by an 
abrupt change in the value of the Hill coefficient or in the constant of 
15 resistance. 

As used herein, a library of mutants refers to a collection of 
plasmids or other vehicles that carry (encode) the gene variants, such that 
individual plasmids or other vehicles carry individual gene variants. 
When a library of proteins is contemplated, it will be so-stated. 
20 As used herein, a "reporter cell" is the cell that "reports", i.e., 

undergoes the change, in response to introduction of the nucleic acid 
infection and, therefore, it is named here a reporter cell. 

As used herein, "reporter" or "reporter moiety" refers to any moiety 
that allows for the detection of a molecule of interest, such as a protein 
25 expressed by a cell. Typical reporter moieties include, for example, 
fluorescent proteins, such as red, blue and green fluorescent proteins. 
For expression in cells, nucleic acid encoding the reporter moiety can be 
expressed as a fusion protein with a protein of interest or under the 
control of a promoter of interest. 
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As used herein, a titering virus increases or decreases the output 
signal from a reporter virus, which is a virus that can be detected, such 
as by a detectable label or signal. 

As used herein, phenotype refers to the physical or other 
5 manifestation of a genotype (a sequence of a gene). In the methods 

herein, phenotypes that result from alteration of a genotype are assessed. 

As used herein, activity refers to the function or property to be 
evolved. An active site refers to a site(s) responsible or that participates 
in conferring the activity or function. The activity or active site evolved 
10 (the function or property and the site conferring or participating in 

conferring the activity) may have nothing to do with natural activities of 
a protein. For example, it could be an 'active site' for conferring 
immunogenicity (immunogenic sites or epitopes) on a protein. 

As used herein, the amino acids, which occur in the various amino 
15 acid sequences appearing herein, are identified according to their known, 
three-letter or one-letter abbreviations (see, Table 1). The nucleotides, 
which occur in the various nucleic acid fragments, are designated with 
the standard single-letter designations used routinely in the art. 

As used herein, amino acid residue refers to an amino acid formed 
20 upon chemical digestion (hydrolysis) of a polypeptide at its peptide 

linkages. The amino acid residues described herein are presumed to be in 
the "L" isomeric form. Residues in the "D" isomeric form, which are so- 
designated, can be substituted for any L-amino acid residue, as long as 
the desired functional property is retained by the polypeptide. NH 2 refers 
25 to the free amino group present at the amino terminus of a polypeptide. 
COOH refers to the free carboxy group present at the carboxyl terminus 
of a polypeptide. In keeping with standard polypeptide nomenclature 
described in J. Biol. Chem., 243:3552-59 (1969) and adopted at 37 
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C.F.R. § § 1.821 - 1.822, abbreviations for amino acid residues are 
shown in the following Table: 

Table 1 

Table of Correspondence 


SYMBOL 


1 -Letter 

3-Letter 

AMINO ACID 

Y 

Tyr 

tyrosine 

G 

Gly 

glycine 

F 

Phe 

phenylalanine 

M 

Met 

» methionine 

A 

Ala 

alanine 

S 

Ser 

serine 

I 

He 

isoleucine 

L 

Leu 

leucine 

T 

Thr 

threonine 

V 

Val 

valine | 

P 

Pro 

proline 

K 

Lys 

lysine 

H 

His 

histidine 

Q 

Gin 

glutamine 

E 

Glu 

glutamic acid 

Z 

Glx 

Glu and/or Gin 

w 

Trp 

tryptophan 

R 

Arg 

arginine 

D 

Asp 

aspartic acid 

N 

Asn 

asparagine 

B 

Asx 

Asn and/or Asp 

C 

Cys 

cysteine 
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SYMBOL 


X 

Xaa 

Unknown or other 


It should be noted that all amino acid residue sequences 
represented herein by formulae have a left to right orientation in the 
5 conventional direction of amino-terminus to carboxyl-terminus. In 

addition, the phrase "amino acid residue" is broadly defined to include the 
amino acids listed in the Table of Correspondence and modified and 
unusual amino acids, such as those referred to in 37 C.F.R. § § 1.821- 
1.822, and incorporated herein by reference. Furthermore, it should be 

10 noted that a dash at the beginning or end of an amino acid residue 

sequence indicates a peptide bond to a further sequence of one or more 
amino acid residues or to an amino-terminal group such as NH 2 or to a 
carboxyl-terminal group such as COOH. 

In a peptide or protein, suitable conservative substitutions of amino 

15 acids are known to those of skill in this art and may be made generally 
without altering the biological activity of the resulting molecule. Those of 
skill in this art recognize that, in general, single amino acid substitutions 
in non-essential regions of a polypeptide do not substantially alter 
biological activity (see, e.g., Watson et al. Molecular Biology of the Gene, 

20 4th Edition, 1987, The Benjamin/Cummings Pub. co., p. 224). 

Such substitutions are preferably made in accordance with those 
set forth in TABLE 2 as follows: 


WO 03/023032 


PCTYIB02/03921 


-18- 


TABLE 2 


Original residue Conservative substitution 

Ala (A) Gly; Ser 

Arg (R) Lys 

5 Asn (N) Gin*' His 

Cys (C) Ser 

Gin (Q) Asn 

Glu (E) Asp 

Gly (G) Ala; Pro 

10 His (H) Asn; Gin 

He (!) Leu; Val 

Leu (L) He; Val 

Lys (K) Arg; Gin; Glu 

Met (M) Leu; Tyr; He 

1 5 Phe (F) Met; Leu; Tyr 

Pro (P) Ala; Gly 

Ser (S) Thr 

Thr (T) Ser 

Trp (W) Tyr 

20 Tyr (Y) Trp; Phe 

Val (V) He; Leu 
Other substitutions are also permissible and may be determined 

empirically or in accord with known conservative substitutions. 

As used herein, nucleic acids include DNA, RNA and analogs 

25 thereof, including protein nucleic acids (PNA) and mixture thereof. 

Nucleic acids can be single or double stranded. When referring to probes 
or primers, optionally labeled, with a detectable label, such as a 
fluorescent or radiolabel, single-stranded molecules are contemplated. 
Such molecules are typically of a length such that they are statistically 

30 unique of low copy number (typically less than 5, preferably less than 3) 
for probing or priming a library. Generally a probe or primer contains at 
least 14, 1 6 or 30 contiguous of sequence complementary to or identical 
to a gene of interest. Probes and primers can be 1 0, 1 4, 1 6, 20, 30, 50, 
100 or more nucleic acid bases long. 

35 As used herein, homologous means about greater than 25% nucleic 

acid sequence identity, preferably 25% 40%, 60%, 80%, 90% or 95%. 
The intended percentage will be specified. The terms "homology" and 
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identity" are often used interchangeably. In general, sequences are 
aligned so that the highest order match is obtained (see, e.g.: 
Computational Molecular Biology, Lesk, A.M., ed., Oxford University 
Press, New York, 1 988; Biocomputing: Informatics and Genome Projects, 
5 Smith, D.W., ed.. Academic Press, New York, 1993; Computer Analysis 
of Sequence Data, Part l t Griffin, A.M., and Griffin, H.G., eds., Humana 
Press, New Jersey, 1994; Sequence Analysis in Molecular Biology, von 
Heinje, G., Academic Press, 1987; and Sequence Analysis Primer, 
Gribskov, M. and Devereux, J., eds., M Stockton Press, New York, 1991; 

10 Carilio et aL (1988) S/AM J Applied Math 45:1073). By sequence 

identity, the number of conserved amino acids are determined by standard 
alignment algorithms programs, and are used with default gap penalties 
established by each supplier. Substantially homologous nucleic acid 
molecules would hybridize typically at moderate stringency or at high 

15 stringency all along the length of the nucleic acid of interest. Also 

contemplated are nucleic acid molecules that contain degenerate codons 
in place of codons in the hybridizing nucleic acid molecule. 

As used herein, a nucleic acid homolog refers to a nucleic acid that 
includes a preselected conserved nucleotide sequence, such as a 

20 sequence encoding a therapeutic polypeptide. By the term "substantially 
homologous", it is meant having at least 80%, preferably at least 90%, 
most preferably at least 95% homology therewith or a less percentage of 
homology or identity and conserved biological activity or function. 
The terms "homology" and "identity" are often used 

25 interchangeably. In this regard, percent homology or identity may be 

determined, for example, by comparing sequence information using a GAP 
computer program. The GAP program uses the alignment method of 
Needleman and Wunsch (J. MoL Biol. 48:443 (1970), as revised by Smith 
and Waterman (Adv. Appl. Math. 2:482 (1981). Briefly, the GAP program 
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defines similarity as the number of aligned symbols (i.e., nucleotides or 
amino acids) which are similar, divided by the total number of symbols in 
the shorter of the two sequences. The preferred default parameters for 
the GAP program may include: (1) a unitary comparison matrix 
5 (containing a value of 1 for identities and 0 for non-identities) and the 
weighted comparison matrix of Gribskov and Burgess, Nucl. Acids Res. 
14:6745 (1986), as described by Schwartz and Dayhoff, eds., ATLAS OF 
PROTEIN SEQUENCE AND STRUCTURE, National Biomedical Research 
Foundation, pp. 353-358 (1979); (2) a penalty of 3.0 for each gap and an 
10 additional 0.10 penalty for each symbol in each gap; and (3) no penalty 
for end gaps. 

Whether any two nucleic acid molecules have nucleotide sequences 
that are, for example, at least 80%, 85%, 90%, 95%, 96%, 97%, 98% 
or 99% /'identical" can be determined using known computer algorithms 

15 such as the "FAST A" program, using for example, the default parameters 
as in Pearson and Lipman, Proc. NatL Acad. ScL USA 55:2444 (1988). 
Alternatively the BLAST function of the National Center for Biotechnology 
Information database may be used to determine identity 

In general, sequences are aligned so that the highest order match 

20 is obtained. "Identity" perse has an art-recognized meaning and can be 
calculated using published techniques. (See, e.g.: Computational 
Molecular Biology, Lesk, A.M., ed., Oxford University Press, New York, 
1988; Biocomputing: Informatics and Genome Projects, Smith, D.W., ed., 
Academic Press, New York, 1993; Computer Analysis of Sequence Data, 

25 Part I, Griffin, A.M., and Griffin, H.G., eds., Humana Press, New Jersey, 
1 994; Sequence Analysis in Molecular Biology, von Heinje, G., Academic 
Press, 1987; and Sequence Analysis Primer, Gribskov, M. and Devereux, 
J., eds., M Stockton Press, New York, 1991). While there exist a number 
of methods to measure identity between two polynucleotide or 
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polypeptide sequences, the term "identity" is well known to skilled 
artisans (Carillo, H. & Lipton, D., SIAM J Applied Math 4S:1073 (1988)). 
Methods commonly employed to determine identity or similarity between 
two sequences include, but are not limited to, those disclosed in Guide to 
5 Huge Computers, Martin J. Bishop, ed., Academic Press, San Diego, 
1994, and Carillo, H. & Lipton, D., SIAM J Applied Math 45:1073 
(1988). Methods to determine identity and similarity are codified in 
computer programs. Preferred computer program methods to determine 
identity and similarity between two sequences include, but are not limited 

10 to, GCG program package (Devereux, J., et aL, Nucleic Acids Research 
72(l):387 (1984)), BLASTP, BLASTIM, FASTA (Atschul, S.F., et aL, J 
Molec Biol 2/5:403 (1990)), and CLUSTALW. For sequences displaying 
a relatively high degree of homology, alignment can be effected manually 
by simply lining up the sequences and manually or visually matching the 

15 conserved portions. 

Therefore, as used herein, the term "identity" represents a 
comparison between a test and a reference polypeptide or polynucleotide. 
For example, a test polypeptide may. be defined as any polypeptide that 
is 90% or more identical to a reference polypeptide. 

20 For the alignments presented herein (see, Fig. 5) for the AAV 

serotype, the CLUSTALW program was employed with parameters set as 
follows: scoring matrix BLOSUM, gap open 10, gap extend 0.1, gap 
distance 40% and transitions/transversions 0.5; specific residue penalties 
for hydrophobic amino acids (DEGKNPQRS), distance between gaps for 

25 which the penalties are augmented was 8, and gaps of extremities 
penalized less than internal gaps. 

As used herein, a "corresponding" position on a protein, such as 
the AAV rep protein, refers to an amino acid position based upon 
alignment to maximize sequence identity. For AAV Rep proteins an 
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alignment of the Rep 78 protein from AAV-2 and the corresponding 
protein from other AAV serotypes (AAV-1, AAV-6, AAV-3, AAV-3B, 
AAV-4, AAV-2 and AAV-5) is shown in Figure 5. The "hit" positions are 
shown in italics. 

5 As used herein, the term at least "90% identical to" refers to 

percent identities from 90 to 100% relative to the reference polypeptides. 
Identity at a level of 90% or more is indicative of the fact that,* assuming 
for exemplification purposes a test and reference polynucleotide length of 
100 amino acids are compared. No more than 10% (i.e., 10 out of 100) 

10 amino acids in the test polypeptide differs from that of the reference 
polypeptides. Similar comparisons may be made between a test and 
reference polynucleotides. Such differences may be represented as point 
mutations randomly distributed over the entire length of an amino acid 
sequence or they may be clustered in one or more locations of varying 

15 length up to the maximum allowable, e.g. 10/100 amino acid difference 
(approximately 90% identity). Differences are defined as nucleic acid or 
amino acid substitutions, or deletions. 

As used herein, it is also understood that the terms substantially 
identical or similar varies with the context as understood by those skilled 

20 in the relevant art. 

As used herein, genetic therapy involves the transfer of 
heterologous nucleic acids to the certain cells, target cells, of a mammal, 
particularly a human, with a disorder or conditions for which such therapy 
is sought. The nucleic acid, such as DNA, is introduced into the selected 

25 target cells in a manner such that the heterologous nucleic acid, such as 
DNA, is expressed and a therapeutic product encoded thereby is 
produced. Alternatively, the heterologous nucleic acid, such as DNA, 
may in some manner mediate expression of DNA that encodes the 
therapeutic product, or it may encode a product, such as a peptide or 


WO 03/023032 


PCT7IB02/03921 


-23- 

RNA that in some manner mediates, directly or indirectly, expression of a 
therapeutic product. Genetic therapy may also be used to deliver nucleic 
acid encoding a gene product that replaces a defective gene or 
supplements a gene product produced by the mammal or the cell in which 
5 it is introduced. The introduced nucleic acid may encode a therapeutic 
compound, such as a growth factor or inhibitor thereof, or a tumor 
necrosis factor or inhibitor thereof, such as a receptor therefor, that is not 
normally produced in the mammalian host or that is not produced in 
therapeutically effective amounts or at a therapeutically useful time. The 

10 heterologous nucleic acid, such as DNA, encoding the therapeutic product 
may be modified prior to introduction into the cells of the afflicted host in 
order to enhance or otherwise alter the product or expression thereof. 
Genetic therapy may also involve delivery of an inhibitor or repressor or 
other modulator of gene expression. 

15 As used herein, heterologous or foreign nucleic acid, such as DNA 

and RNA, are used interchangeably and refer to DNA or RNA that does 
not occur naturally as part of the genome in which it is present or which 
is found in a location or locations in the genome that differ from that in 
which it occurs in nature. Heterologous nucleic acid is generally not 

20 endogenous to the cell into which it is introduced, but has been obtained 
from another cell or prepared synthetically. Generally, although not 
necessarily, such nucleic acid encodes RNA and proteins that are not 
normally produced by the cell in which it is expressed. Any DNA or RNA 
that one of skill in the art would recognize or consider as heterologous or 

25 foreign to the cell in which it is expressed is herein encompassed by 

heterologous DNA. Heterologous DNA and RNA may also encode RNA or 
proteins that mediate or alter expression of endogenous DNA by affecting 
transcription, translation, or other regulatable biochemical processes. 
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Examples of heterologous nucleic acid include, but are not limited to, 
nucleic acid that encodes traceable marker proteins, such as a protein 
that confers drug resistance, nucleic acid that encodes therapeutically 
effective substances, such as anti-cancer agents, enzymes and hormones, 
5 and DNA that encodes other types of proteins, such as antibodies. 

Hence, herein heterologous DNA or foreign DNA, includes a DNA 
molecule not present in the exact orientation and position as the 
counterpart DNA molecule found in the genome. It may also refer to a 
DNA molecule from another organism or species {/.e., exogenous). 

10 As used herein, a therapeutically effective product introduced by 

genetic therapy is a product that is encoded by heterologous nucleic acid, 
typically DNA, that, upon introduction of the nucleic acid into a host, a 
product is expressed that ameliorates or eliminates the symptoms, 
manifestations of an inherited or acquired disease or that cures the 

15 disease. 

As used herein, isolated with reference to a nucleic acid molecule 
or polypeptide or other biomolecule means that the nucleic acid or 
polypeptide has separated from the genetic environment from which the 
polypeptide or nucleic acid were obtained. It may also mean altered from 

20 the natural state. For example, a polynucleotide or a polypeptide naturally 
present in a living animal is not "isolated," but the same polynucleotide or 
polypeptide separated from the coexisting materials of its natural state is 
"isolated", as the term is employed herein. Thus, a polypeptide or 
polynucleotide produced and/or contained within a recombinant host cell 

25 is considered isolated. Also intended as an "isolated polypeptide" or an 
"isolated polynucleotide" are polypeptides or polynucleotides that have 
been purified, partially or substantially, from a recombinant host cell or 
from a native source. For example, a recombinantly produced version of 
a compounds can be substantially purified by the one-step method 
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described in Smith and Johnson, Gene 57:31-40 (1988). The terms 
isolated and purified are sometimes used interchangeably. 

Thus, by "isolated" it is meant that the nucleic is free of the coding 
sequences of those genes that, in the naturally-occurring genome of the 
5 organism (if any) immediately flank the gene encoding the nucleic acid of 
interest. Isolated DNA may be single-stranded or double-stranded, and 
may be genomic DNA, cDNA, recombinant hybrid DNA, or synthetic 
DNA. It may be identical to a native DNA sequence, or may differ from 
such sequence by the deletion, addition, or substitution of one or more 

10 nucleotides. 

Isolated or purified as it refers to preparations made from biological 
cells or hosts means any cell extract containing the indicated DNA or 
protein including a crude extract of the DNA or protein of interest. For 
example, in the case of a protein, a purified preparation can be obtained 

15 following an individual technique or a series of preparative or biochemical 
techniques and the DNA or protein of interest can be present at various 
d egrees Q f purity in these preparations. The procedures may include for 
example, but are not limited to, ammonium sulfate fractionation, gel 
filtration, ion exchange change chromatography, affinity chromatography, 

20 density gradient centrifugation and electrophoresis. 

A preparation of DNA or protein that is "substantially pure" or 
"isolated" should be understood to mean a preparation free from naturally 
occurring materials with which such DNA or protein is normally 
associated in nature. "Essentially pure" should be understood to mean a 

25 "highly" purified preparation that contains at least 95% of the DNA or 
protein of interest. 

A cell extract that contains the DNA or protein of interest should be 
understood to mean a homogenate preparation or cell-free preparation 
obtained from cells that express the protein or contain the DNA of 


WO 03/023032 


PCT/IB02/03921 


-26- 

interest. The term "cell extract" is intended to include culture media, 
especially spent culture media from which the cells have been removed. 

As used herein, receptor refers to a biologically active molecule 
that specifically binds to (or with) other molecules. The term "receptor 
5 protein" may be used to more specifically indicate the proteinaceous 
nature of a specific receptor. 

As used herein, recombinant refers to any progeny formed as the 
result of genetic engineering. 

As used herein, a promoter region refers to the portion of DNA of a 
10 gene that controls transcription of the DNA to which it is operatively 

linked. The promoter region includes specific sequences of DNA that are 
sufficient for RNA polymerase recognition, binding and transcription 
initiation. This portion of the promoter region is referred to as the 
promoter. In addition, the promoter region includes sequences that 
15 modulate this recognition, binding and transcription initiation activity of 
the RNA polymerase. These sequences may be cis acting or may be 
responsive to trans acting factors. Promoters, depending upon the nature 
of the regulation, may be constitutive or regulated. 

As used herein, the phrase "operatively linked" generally means the 
20 sequences or segments have been covalently joined into one piece of 
DNA, whether in single or double stranded form, whereby control or 
regulatory sequences on one segment control or permit expression or 
replication or other such control of other segments. The two segments 
are not necessarily contiguous. For gene expression a DNA sequence and 
25 a regulatory sequence(s) are connected in such a way to control or permit 
gene expression when the appropriate molecules, e.g., transcriptional 
activator proteins, are bound to the regulatory sequence(s). 

As used herein, production by recombinant means by using 
recombinant DNA methods means the use of the well known methods of 


WO 03/023032 


PCTAB02/03921 


-27- 

molecular biology for expressing proteins encoded by cloned DNA, 
including cloning expression of genes and methods, such as gene 
shuffling and phage display with screening for desired specificities. 

As used herein, a splice variant refers to a variant produced by 
5 differential processing of a primary transcript of genomic DNA that results 
in more than one type of mRNA. 

As used herein, a composition refers to any mixture of two or more 
products or compounds. It may be a solution, a suspension, liquid, 
powder, a paste, aqueous, non-aqueous or any combination thereof 
10 As used herein, a combination refers to any association between 

two or more items. 

As used herein, substantially identical to a product means 
sufficiently similar so that the property of interest is sufficiently 
unchanged so that the substantially identical product can be used in place 

15 of the product. 

As used herein, the term "vector" refers to a nucleic acid molecule 
capable of transporting another nucleic acid to which it has been linked. 
One type of preferred vector is an episome, i.e., a nucleic acid capable of 
extra-chromosomal replication. Preferred vectors are those capable of 

20 autonomous replication and/or expression of nucleic acids to which they 
are linked. Vectors capable of directing the expression of genes to which 
they are operatively linked are referred to herein as "expression vectors". 
In general, expression vectors of utility in recombinant DMA techniques 
are often in the form of "plasmids" which refer generally to circular 

25 double stranded DNA loops which, in their vector form, are not bound to 
the chromosome. "Plasmid" and "vector" are used interchangeably as the 
plasmid is the most commonly used form of vector. Other such other 
forms of expression vectors that serve equivalent functions and that 
become known in the art can be used subsequently hereto. 
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As used herein, vector is also used interchangeable with "virus 
vector" or "viral vector". In this case, which will be clear from the 
context, the "vector" is not self-replicating. Viral vectors are engineered 
viruses that are operatively linked to exogenous genes to transfer (as 
5 vehicles or shuttles) the exogenous genes into cells. 

As used herein, transduction refers to the process of gene transfer 
and expression into mammalian and other cells mediated by viruses. 
Transfection refers to the process when mediated by plasmids. 

* As used herein, "polymorphism" refers to the coexistence of more 
10 than one form of a gene or portion thereof. A portion of a gene of which 
there are at least two different forms, i.e., two different nucleotide 
sequences, is referred to as a "polymorphic region of a gene". A 
polymorphic region can be a single nucleotide, referred to as a single 
nucleotide polymorphism (SNP), the identity of which differs in different 
15 alleles. A polymorphic region can also be several nucleotides in length. 

As used herein, "polymorphic gene" refers to a gene having at least 
one polymorphic region. 

As used herein, "allele", which is used interchangeably herein with 
"allelic variant" refers to alternative forms of a gene or portions thereof. 
20 Alleles occupy the same locus or position on homologous chromosomes. 
When a subject has two identical alleles of a gene, the subject is said to 
be homozygous for the gene or allele. When a subject has two different 
alleles of a gene, the subject is said to be heterozygous for the gene. 
Alleles of a specific gene can differ from each other in a single nucleotide, 
25 or several nucleotides, and can include substitutions, deletions, and 

insertions of nucleotides. An allele of a gene can also be a form of a gene 
containing a mutation. 

As used herein, the term "gene" or "recombinant gene" refers to a 
nucleic acid molecule comprising an open reading frame and including at 
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least one exon and (optionally) an intron sequence. A gene can be either 
RNA or DNA. Genes may include regions preceding and following the 
coding region {leader and trailer). 

As used herein, "intron" refers to a DNA sequence present in a 
5 given gene which is spliced out during mRNA maturation. 

As used herein, "nucleotide sequence complementary to the 
nucleotide sequence set forth in SEQ ID NO: x" refers to the nucleotide 
sequence of the complementary strand of a nucleic acid strand having 
SEQ ID NO: x. The term "complementary strand" is used herein 

10 interchangeably with the term "complement". The complement of a 
nucleic acid strand can be the complement of a coding strand or the 
complement of a non-coding strand. When referring to double stranded 
nucleic acids, the complement of a nucleic acid having SEQ ID NO: x 
refers to the complementary strand of the strand having SEQ ID NO: x or 

15 to any nucleic acid having the nucleotide sequence of the complementary 
strand of SEQ ID NO: x. When referring to a single stranded nucleic acid 
having the nucleotide sequence SEQ ID NO: x, the complement of this 
nucleic acid is a nucleic acid having a nucleotide sequence which is 
complementary to that of SEQ ID NO: x. 

20 As used herein, the term "coding sequence" refers to that portion 

of a gene that encodes an amino acid sequence of a protein. 

As used herein, the term "sense strand" refers to that strand of a 
double-stranded nucleic acid molecule that has the sequence of the 
mRNA that encodes the amino acid sequence encoded by the double- 

25 stranded nucleic acid molecule. 

As used herein, the term "antisense strand" refers to that strand of 
a double-stranded nucleic acid molecule that is the complement of the 
sequence of the mRNA that encodes the amino acid sequence encoded 
by the double-stranded nucleic acid molecule. 
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As used herein, an array refers to a collection of elements, such as 
nucleic acid molecules, containing three or more members. An 
addressable array is one in which the members of the array are 
identifiable, typically by position on a solid phase support or by virtue of 
5 an identifiable or detectable label, such as by color, fluorescence, 

electronic signal [i.e. radiofrequency (RF), microwave or other frequency 
that does not substantially alter the interaction of the molecules of 
interest), bar code or other symbology, chemical or other such label. 
Hence, in general the members of the array are immobilized to discrete 

10 identifiable loci on the surface of a solid phase or directly or indirectly 
linked to or otherwise associated with the identifiable label, such as 
affixed to a microsphere or other particulate support (herein referred to as 
beads) and suspended in solution or spread out on a surface. 

As used herein, a support (also referred to as a matrix support, a 

15 matrix, an insoluble support or solid support) refers to any solid or 

semisolid or insoluble support to which a molecule of interest, typically a 
biological molecule, organic molecule or biospecific ligand is linked or 
contacted. Such materials include any materials that are used as affinity 
matrices or supports for chemical and biological molecule syntheses and 

20 analyses, such as, but are not limited to: polystyrene, polycarbonate, 
polypropylene, nylon, glass, dextran, chitin, sand, pumice, agarose, 
polysaccharides, dendrimers, buckyballs, polyacrylamide, silicon, rubber, 
and other materials used as supports for solid phase syntheses, affinity 
separations and purifications, hybridization reactions, immunoassays and 

25 other such applications. The matrix herein can be particulate or can be 
in the form of a continuous surface, such as a microtiter dish or well, a 
glass slide, a silicon chip, a nitrocellulose sheet, nylon mesh, or other 
such materials. When particulate, typically the particles have at least one 
dimension in the 5-10 mm range or smaller. Such particles, referred 
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colledtively herein as "beads", are often, but not necessarily, spherical. 
Such reference, however, does not constrain the geometry of the matrix, 
which may be any shape, including random shapes, needles, fibers, and 
elongated. Roughly spherical "beads", particularly microspheres that can 
5 be used in the liquid phase, are also contemplated. The "beads" may 
include additional components, such as magnetic or paramagnetic 
particles (see, e.g.,, Dyna beads (Dynal, Oslo, Norway)) for separation 
using magnets, as long as the additional components do not interfere with 
the methods and analyses herein. 

10 As used herein, matrix or support particles refers to matrix 

materials that are in the form of discrete particles. The particles have any 
shape and dimensions, but typically have at least one dimension that is 
100 mm or less, 50 mm or less, 10 mm or less, 1 mm or less, 100 //m or 
less, 50 //m or less and typically have a size that is 1 00 mm 3 or less, 50 

15 mm 3 or less, 10 mm 3 or less, and 1 mm 3 or less, 100//m 3 or less and may 
be order of cubic microns. Such particles are collectively called "beads." 

As used herein, the abbreviations for any protective groups, amino 
acids and other compounds, are, unless indicated otherwise, in accord 
with their common usage, recognized abbreviations, or the IUPAC-IUB 

20 Commission on Biochemical Nomenclature (see, (1972) Biochem. 
/ 7:942-944). 

B. High Throughput Process 

Provided herein are high throughput processes for the generation of 
and identification of proteins that exhibit desired phenotypes. The 
25 processes include methods that are particularly adapted for high 

throughput protocols, which require accurate methods for identifying 
modified proteins. 

A general directed evolution process includes the following steps: 


1 
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1 . Generation of diversity at the nucleic acid level, on the gene 
to be 'evolved' 

2. phenotypic characterization of the gene variants generated; 

and 

5 3. identification of optimized gene variants. 

The processes provided herein effect these steps that can be 
performed in a high throughput format (see, FIGURES 1) that is optionally 
automated. A distinguishing characteristic of the processes provided 
herein, is that each candidate nucleic acid molecule is separately 

10 generated and screened. In an automated process at least some of the 
steps are performed without human intervention and are generally 
controlled by software. Most, if not all steps, are performed in 
addressable formats, such as at discrete locations in or on solid supports, 
such as microtiter plates or in other addressable formats, such as linked 

15 to coded supports. The supports can be electronically, physically, 

chemically or otherwise identifiable, such as by an identifiable symbology, 
including a bar code, or can be color coded. 

1 . Generation of Diversity using a semi-rational approach 

A semi-rational approach to creating diversity or evolving genes is 

20 provided herein. The goal is to create diversity but to decrease the 
number of molecules to be screened. By reducing the numbers, the 
molecules can be screened in high throughput format molecule-by- 
molecule (or groups thereof). 

Generation of diversity at the nucleic acid level, in principle, can be 

25 accomplished by a number of diverse technologies like mutagenesis 
(either site-directed or random), recombination, shuffling and de-novo 
synthesis. These different technologies differ in the degree of diversity 
they generate as well as in the minimal length of the unitary change they 
can introduce (from single base to large domains). The outcome of step 1 
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is a collection of diverse, although highly related, molecules that 

constitutes what is called a 'library'. 

This step is crucial, since it provides the initial conditions for the 

entire process and is determinative of the outcome. The chances of 

5 finding an optimized gene version in a library is a function of the total 

diversity present in the library. In addition, the type of diversity 

introduced (such as, but not limited to, single point mutations, multiple 

point mutations, scarce small rearrangements, recombination of large 

domains, multiple shuffling) condition the outcome, particularly with 

10 respect to the generation of new variants compared to the original gene, 

and the probability that the new variants, not only exhibit the "evolved" 

function or property, but also work in their natural biological networks 

where they are expected to act by interacting, recognizing, and/or being 

recognized, by a large panoply of other proteins and other molecules. 

1 5 Rapid discovery of protein variants at the amino acid level by 

rational mutagenesis (aa-scan) 

A method, referred to herein as an amino-acid scan method for 
directed evolution, is provided herein for generating protein variants. This 
method can be performed on an entire protein or selected domains 

20 thereof, or can be used to identify benchmark sequences, such as 

functional domains, and, for example, recombine them as exchangeable 
units or restrict the diversity to limited or specific regions of the protein. 
Not only can this method be used with the processes provided herein, but 
it also has applications for any methods that use such variants or require 

25 generation of such variants, such as, but not limited to, searches for 
consensus sequences and homology regions that are used in functional 
genomics, functional proteomics; comparative modeling in protein 
crystallography and protein modeling; searches for natural diversity, 
(e.g., directed evolution methods in 6,171,820, 6,238,884, 6,174,673, 
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6,057,103, 6,001,574, 5,763,239, ); exon- or family-shuffling-based 
diversity {e.g., directed evolution using gene shuffling (see, e.g., U.S. 
Patent Nos. 6,096,548, 6,117,679, 6,165,793, 6,180,406, 6,132,970); 
the optimization of only the CDRs regions {e.g., directed evolution of 
5 antibodies see., e.g., U.S. Patent Nos. 5,723,323, 6,258,530, 5770,434, 
5,862,514) and other methods (see, e.g., U.S. Patent Nos. 5,837,500, 
5,571,698, 6,156,509). 

The amino-acid scanning-based method provided herein has 
advantages that prior methods do not have. For example, prior methods 

10 are based upon the underlying assumption that there are parts of the 
molecule (gene or protein) that are sufficiently adapted to perform their 
respective function, and further changes are not advantageous. Such 
methods do not look at total potential plasticity of a given molecule, but 
at the plasticity still permitted while keeping some basic functions in 

15 place. By choosing this route, however, additional potential variation is 
missed. The potential in the intrinsic plasticity of those regions that are 
presumed 'preserved' is lost. For instance, methods {e.g., those in U.S. 
Patent Nos. 6,171,820, 6,238,884, 6,174,673, 6,057,103, 6,001,574, 
5,763,239) that use natural diversity can miss the potential plasticity 

20 within those regions that are naturally 'conserved', i.e. there where there 
is no natural diversity. Methods that rely on exon- or family-shuffling- 
based diversity can miss the potential plasticity within regions contained 
in the shuffled fragments, i.e. within the fragments exchanged as a block. 
The method provided herein in contrast is sufficiently flexible to 

25 create mutants at a variety of levels, including at the single amino acid 
level; i.e. the method can generate mutants that differ from each other at 
a single amino and not at a larger block level. The challenge solved by 
the method herein is to generate diversity at the single amino acid level, 
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without moving too close to a pure 'random' approach, which results in 
an intractable number of mutants. 

The method provided herein is based on the premise that there are 
single amino acids or small blocks of sequence of amino acids that are 
5 either (1) directly involved in the activities that the methods 'evolve' 
(these amino acids would be at or close to the 'active sites' of the 
protein); or (2) directly involved in maintaining within the protein the intra- 
molecular environment that allows the active site{s) to stay active. 

Potential plasticity at the amino acid level can be exploited if amino 

10 acids or blocks of amino acids directly involved in the active sites for the 
activity to be evolved are known. Often they are not known. The 
problem that is solved herein, however, is how to exploit the potential 
plasticity at the amino acid level when nothing is known about the 
structure of the protein in question or about the position of its single or 

15 numerous active sites. 

The technology referred to herein as amino acid-scanning has been 
used to precisely identify those amino acids directly involved in the active 
sites of some enzymes and receptors (see, e.g., Becl-Sickinger et al. 
(1994) Eur. J. Biochem. 223:947-958; Gibbs et al. (1991U Biol. Chem. 

20 255:8923-8931; Matsushita et al. (2000) J. Biol. Chem. 275:11044- 
1 1049) but has not been employed for directed evolution or for the 
generation of diversity. The amino acid scan as practiced in the prior art 
is used to produce a set of protein mutants, often within the region 
suspected to contain the active site(s), such that in each individual 

25 mutant a selected residue, such as Ala, replaces a different amino acid. 
Ala or other neutral amino acids generally, although not necessarily, is 
selected as a replacement amino acid since, except in instances in which 
the replaced amino acid is directly involved in an active site, it should 
have a neutral effect on the protein activity and not disturb the native 
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secondary structure of the protein. In instances in which the replaced 
amino acid is directly involved in an active site the activity of that site will 
be lost or altered. Amino-acid scanning, such as Ala-scanning, has been 
successfully applied for the identification of active sites in a number of 
5 proteins, and has been performed in computer-based rational drug design 
methods. Other amino acids, particularly amino acids that have a neutral 
effect, such as Glycine, can also be used. 

The amino acid scanning method is employed herein for the 
generation of the mutant proteins for screening for identification of sites 

10 or loci in a target protein or regions in a protein that alter a selected 

activity. In performing this method, the amino acids are each replaced, 
one-by-one along the full-length of the protein, or one-by-one. in pre- 
selected domains, such as domains that possess a desired activity or 
exhibit a particular function. Once sites of interest are identified other 

15 methods for generating diversity from the resulting molecules can be 
employed or the further steps of the method provided herein can be 
performed. 

The method includes the following steps: x 
(1) Identification of the active site(s) on the full 

20 length protein sequence. In one embodiment a full-length amino acid- 
scan, typically, although not necessarily, an Ala-scan, or the identification 
and positioning of the active site(s) on proteins of either known or 
unknown function. For purposes herein, an active site is not necessarily 
the natural active site involved in the natural activity of a target protein, 

25 but those amino acids involved in the activities of the proteins under 

'directed evolution' with the purpose of either gain, improvement or loss 
of function. 

The whole process of the 'identification of the active site(s) on the 
full length protein sequence requires the following sub-steps: 
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a. Generation of a mutant library (on the 
gene to be evolved) in which each individual mutant contains a single 
mutation located at a different amino acid position and that includes a 
systematic replacement of the native amino acid by Ala or any other 

5 amino acid (always the same throughout the entire protein sequence); 

b. phenotypic characterization of the 
individual mutants, one-by-one and assessment of mutant protein activity; 

c. identification of those mutants that 

10 display an alteration, typically a decrease, in the selected protein activity, 
thus, indicating that amino acids directly involved in the active site(s) 
have been hit. The aa positions whose aa-scan mutations display an 
alteration, typically a loss or decrease, in activity are named HITS. 
The identification of the active site(s) (HITS) is thus, by this 

15 method, made in a completely unbiased manner. There are no 

assumptions about the specific structure of the protein in question nor 
any knowledge or assumptions about the active site(s). The results of the 
amino acid scan identify such sites. 

Once the active site(s) (the HITS) has(ve) been identified, those 

20 amino acids either at or surrounding the active sites, such as within 1 , 2, 
3, ... 10, 20 or any selected regions, as the unitary elements of 
exchange and generate diversity either at or around one of those sites or 
as a combined diversity at several sites at a time can be assessed. 
This process includes the following steps: 

25 a. Generation of a new mutant library (on the gene to be 

evolved) in which each individual mutant contains either single or multiple 
mutations located at (or surrounding) a specific active site (a HIT) position 
detected by the precedent aa-scan process. In the example these 
mutations include replacement, in each individual mutant, of the native 
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amino acid located either at (or surrounding) the HIT position by one of all 
other possible amino acids, such that, in the library, and at (or 
surrounding) each HIT position the native amino acid has been replaced 
by all possible amino acids. 
5 b. Identification of those mutants that display an 

increase in protein activity, thus indicating that a new sequence at or 
surrounding an active site has been identified with higher activity 
compared to the native sequence. These optimized sequences are named 
LEADS. 

10 The process can be repeated as many times as desired, in search 

for new combinations of optimal amino acids at (or surrounding) the 
different HIT sites. Each time, the process includes the steps of 
generating of a new mutant library (of the gene to be evolved) in which 
each individual mutant contains either single or multiple mutations located 

15 at (or surrounding) a specific active site (a HIT) position; phenotypic 

characterization of the individual mutants, one-by-one and assessment of 
mutant protein activity; and identification of those mutants that display an 
increase in protein activity, thus indicating that a new sequence at or 
surrounding an active site has been identified with higher activity 

20 compared to the native sequence. These optimized sequences are again 
named LEADS (second generation LEADS). 

2. Phenotypic characterization of the gene variants 
This step requires the expression of the gene variants in order to 
allow them to manifest their respective phenotypes. Gene expression can 

25 be accomplished by different means: in vitro, in reconstituted systems or 
in vivo in cellular systems, including bacterial and eukaryotic cells. For all 
exemplification purposes, reference is made to in vivo systems. Those 
of skill in the art can readily adapt these methods for in vitro systems, 
including those using biochemical assays. 


WO 03/023032 


PCT/IB02/03921 


-39- 

This step is a crucial step for several reasons: 

(a) Expression system and protein processing. 
Depending on the system used (either bacteria or eukaryotic cells), 

as well as on the specific gene to be 'evolved', the variant proteinsmay 
5 or may not be appropriately processed, especially when post-translational 
modifications are involved, and therefore be able or not to elicit their 
potential activity. Consequently, the expression system (bacteria vs. 
eukaryotic cells) has to be carefully chosen. 

(b) Standardization of the expression system. 

10 The technologies available for gene transfer and expression into 

either bacteria or eukaryotic (let's consider mammalian) cells widely vary 
in their intrinsic efficiencies. While it is very easy to efficiently transfer 
and express genes in bacteria by chemical/physical methods 
(transformation), that is not the case for mammalian cells, where the 

15 transformation (here called transfection) process is inefficient and 

unreliable, especially when reproducibility and robustness are necessary in 
miniature, large number and small scale high throughput settings like 
those necessary to analyze gene variant libraries. Therefore, when 
transfection is used on mammalian cells, the specific activity measured 

20 for the individual variants in the library most probably does not accurately 
reflect the real specific activity of the molecules involved. As provided 
herein, transduction, the process of gene transfer and expression into 
mammalian and other cells mediated by viruses, overcomes the limitations 
of transfection. 

25 (3) Characterization. 

A distinction must be made between the 'expression' of the gene 
variants and their 'phenotypic characterization'. The expression system 
(either bacteria or mammalian cells) is only the vehicle to convert the 
gene variants into protein variants. The phenotypic characterization is 
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performed on the protein variants, and may have nothing to do, 
depending on the specific system under study, with the cellular system 
used to express the variants. The phenotypic characterization requires the 
use of specific assays (either biochemical (cell-free) or cell-based assays) 
5 in which the activity of the different cell mutants is challenged and 
assessed. In addition to the implications discussed below, these assays 
must be designed in such a way that they reflect the final environment in 
which the 'evolved' protein is expected to act. As an example, when 
optimizing an enzyme to be used in an artificial industrial setting, the 

10 assay should reproduce those conditions {temperature, pH, media 

composition...) of the real-life industrial reaction mixture, which may be 
relatively easy to do. When the final destination of the 'evolved' protein 
is a complex biological setting, such as the intracellular environment, the 
extra-cellular milieu (example: circulating proteins) or the structure of a 

15 virus, the necessary assay (s) may be quite difficult to setup. With a few 
exceptions, most of the work done so far on directed evolution has been 
made on simple enzymes for which all the necessary settings are 
relatively easy to implement. 

Methods for accurately titering viruses 

20 Much progress in gene therapy, genomics, biotechnology and, in 

g enera | # biomedical sciences, depends on the ability to generate and 
analyze large numbers and small amounts of specific viruses. High 
throughput technologies are employed in disciplines such as functional 
genomics and gene therapy in which the use of viruses plays a key role 

25 for the efficient transfer and analysis of large collections of genes or 
libraries. Also, virus samples and biomedical samples containing viruses 
are routinely analyzed in thousands of hospitals, health centers, 
academic labs and biotech settings. 
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Furthermore in processes herein, accurate titration can be important 
in at least two steps in the process. After preparation of the viruses with 
the mutated variant, and prior to screening, it is necessary to know the 
concentration of titer of the viruses in the sample, so that results among 
5 the samples can be compared. The methods in this section designated 
Real Time Virus Titering (RTVT") and (TREE) are advantageously used. 

The methods in this section are also used in data analysis when 
measuring the output signal. As described below, output signal can be 
assessed by a Hill analysis or a second order polynomial or other 

10 algorithm that describes the interaction of biological molecules in complex 
system. In addition, where the output signal is actually the number of 
viral particles or ip produced, the methods in this section RTVT and TREE 
are advantageously used. 

Prior art virus titration methods (RCA, dRA...),for determining the 

15 amount of virus present in a biological sample, are based on the 

assessment of some kind of output signal, such as cytopathic effect, lysis 
or plaques and cell fusion focuses, induced in a reporter cell following a 
fixed time after infection with varying concentrations of the sample 
containing the virus. The lowest concentration of the sample at which no 

20 signal can be measured is taken as the titer of the virus in that sample. 
These approaches are known as "serial dilution" or "limiting dilution" 
methods. In limiting dilution methods, one single virus concentration, 
measured at a given time end-point gives rise to a single measurement of 
the output signal. These methods tend to be destructive in that 

25 assessment destroys the reaction so that only a single measurement can 
be taken on a sample.. 
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Real Time Virus Titering (RTVT") 
When a virus infects a cell, the infected cell undergoes a number 
of changes that can be followed over time and quantified. Such changes 
are designated herein as the "output signal". The cell reports an output 
5 signal in response to the infection and, therefore, it is named here a 
reporter cell. One such output signal, is, for example, the expression of 
the genes carried by the virus (whether they are viral genes or 
exogenous genes (transgenes)). The output signal (for instance the level 
of expression of those genes) develops over time and depends, mainly, 

10 on two factors: i) the time point ("t") at which its level is measured after 
infection and ii) the amount of virus infecting the cell; i.e. the 
concentration of the virus preparation used to infect the cell ("s"). 
The output signal, at a given time point after infection, will be higher for 
higher concentrations of the virus infecting the reporter cells; and for any 

15 given concentration of virus, the output signal increases with time after 
infection until it generally reaches a plateau level. 

Real Time Virus Titering (RTVT™) published as International PCT 
application No. WO 01/186291, which is based on PCT/FR01/01 366 and 
EXAMPLES below) uses non-destructive methods for the assessment of 

20 output signal. Real Time Virus Titering is a viral titration method based on 
the kinetic analysis of the development of the output signal in 
virus-infected cells,, tested at a single concentration of virus or biological 
sample, instead of fixing the time point after infection and varying the 
concentration of the sample as is done in limiting dilution methods, in the 

25 Real Time Virus Titering RTVT" method, a fixed concentration of virus is 
used and the generation of a signal over time is assessed. Hence the 
signal is measured as a function of time at a single virus concentration. In 
this situation, a single virus sample (concentration), whose output signal 
is measured at a number of time points, can give rise to as many 
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measurements of the output signal as needed, and, eventually to a 
continuous, over time, assessment of the signal in real time. 

Real Time Virus Titering RTVT" can be advantageously used in high 
throughput methods in which large numbers of biological samples are 
5 analyzed at the same time. This is the case, for instance, when titering 
viruses in a virus library. Limiting dilution methods rely on the output 
signal generated by a number of dilutions of each individual sample. If, for 
example, 1 0 dilutions (or experimental points) of each virus are used for a 
titration using a limiting dilution method, the analysis of a library 

10 containing 10,000 viruses require analysis of 10 5 (i.e., 10 x 10,000) 

experimental points. The Real Time Virus Titering RTVT - method requires 
only one dilution per sample, thereby requiring 10-fold fewer experimental 
points than a limiting dilution method. For a Real Time Virus Titering 
RTVT" titering system, the time (t£) necessary for the output signal to 

1 5 reach a reference value (fi) is a direct function of the concentration of 
virus. Thus, Xfi can be used to directly determine the concentration of the 
virus. 

A limitation of the Real Time Virus Titering RTVT" limiting dilution 

titering method, however, is that not ail the viruses (nor the genes carried 

20 by the viruses) generate a readily measured output signal that can be 

followed over time using non-destructive methods. 

Tagged Replication and expression enhancement 
(TREE) 

A method for titering designated Tagged Replication and Expression 
25 Enhancement Technology (TREE") is provided herein. This system 
includes: i) a cell, ii) a reporter virus carrying a reporter gene, whose 
activity can be followed over time by a non- destructive method (i.e., 
fluorescence), iii) the virus (or virus library to be titered), herein referred 
to as the "titering virus". The elements are employed such that the virus 
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to be titered interferes with any output signal generated by the reporter 
virus, leading to either decrease or increase in the amount of that signal. 
The higher the amount of virus to be titered, the higher is the interference 
with the reporter virus and output signal. In the absence of virus to be 
5 titered, the kinetics of the output signal generated by the reporter virus 
are followed using the Real Time Virus Titering RTVT" titering method. In 
the presence of increasing amounts of the virus to be titered the output 
signal takes longer (or shorter) to develop as a function of the amount of 
virus to be titered. 

10 Using the TREE" titering method, tj3, the time necessary for the 

output signal to reach a reference value (0) is a direct function of the 
concentration of virus and, therefore, tJ3 can be used to determine the 
concentration of the virus to be titered. It is demonstrated herein (see the 
EXAMPLES) that when using the TREE system for titering, once an 

15 appropriate reference value (/?) is determined for the output signal 
generated from the reporter virus, the time tJ3 is a function of the 
concentration of the virus being titered (see Example). Therefore, the 
concentration (titer) of the virus to be titered, can be assessed by 
assessing the change induced in tfi by an aliquot of the virus to be titered. 

20 In a calibrated TREE titering assay, only one aliquot virus to be titered is 
needed to determine its titer, which is determined by measuring the shift 
in the tfi of the system. The only condition is that the virus to be titered 
must "interfere" [i.e., increase or decrease) the output signal of the 
reporter virus. 

25 A calibration curve representing t/ff vs. the amount of virus to be 

titered is obtained using aliquots of a reference batch of virus of known 
titer (previously determined using any titering procedure). The calibration 
curve can then be used to determine the amount of virus in a sample of 
unknown titer, based on the change caused by an aliquot of the sample 
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on the t£ of the system and the corresponding titer read from the 
calibration curve. 

3. Identification of gene variants 
There are at least two considerations in this step: 
5 (a) Selection vs. screening. 

Depending on the specific protein involved, and under certain and 
very specific assay conditions, those variants that have been 'evolved' 
may elicit a selective advantage over the native version. This situation 
represents the most simple case: the cells (bacteria or mammalian) 

10 expressing the library of protein variants, as a pool or mixture, can be 

simply exposed to the selective conditions which by themselves will allow 
to put in evidence the best optimized variants. This situation is however 
very rare and difficult to achieve. It can be hardly believed that for any 
protein that one may want to optimize, a suitable 'selective' assay could 

15 be set up. For the vast majority of the cases, selection will not be 
possible. Therefore pools of molecules cannot be used, because the 
specific readouts of the assays could not be attributed to individual 
variants. When the simplistic selection approach is not possible, then two 
things are absolutely needed: (a) a 'one-by-one' approach, i.e. each 

20 individual variant must be physically separated from the others and its 
activity tested independently; (b) an accurate and quantitative analysis 
that can distinguish slight differences in activity among the different 
variants along a wide range of performance values. 

(b) Accurate quantitative analytics 

25 When selection is not possible, the optimized variants must be 

distinguished from the native variant otherwise. The different degrees of 
optimization among the different variants in the library should, in addition, 
be distinguished if those variants showing the highest optimization level 
are to be identified. A powerful quantitative analytical protocol is then 
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mandatory. These analytics should be able to attribute quantitative 
features {on the activity tested in the specific assay) to each of the 
variants tested and to rank them according to their individual 
performance. This requires, in addition, that each variant in the library is 
5 assayed individually; the use of pools or mixtures of molecules would 
hamper the ability to identify the right variants. 

For such analysis, the output signal can be assessed by a Hill 
analysis (see Examples and (published International PCT application No. 
WO 01/44809 based on PCT n° PCT/FROO/03503),.) or a second order 

10 polynomial (see, Examples and (Drittanti eta/. (2000) Gene Then 7: 924- 
929)) or other algorithm that describes the interaction of biological 
molecules in complex system, such as the interaction between cells and 
biological agents. In addition, where the output signal is actually the 
number of viral particles or ip produced, the methods in this designated 

15 Real Time Virus Titering (RTVT") and Tagged Replication and expression 
enhancement (TREE™) are advantageously used (for a discussion of 
RTVr, see, International PCT application No. PCT/FR0 1/01 366 published 
as International PCT application No. WO 01/186291 and the EXAMPLES 
below) or a refinement of that method provide herein and designated 

20 Tagged Replication and expression enhancement (TREE™) described 
above and in the examples. 
C. Practice of the process 

In one embodiment, the process provided herein includes the 
following steps. 

25 1 . Generation of diversity or source of existing diversity 

Generation of a plasmid library containing the genetic variants. The 
genetic variants are physically separated from each other. Any model 
such as, but not limited to, amino acid scanning, mutagenesis, or 
recombination may be used to generate the plasmid library. 
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2. Expression of the genetic variants 
Any method for expression of variants is contemplated. In 
particular the following alternatives are particularly suitable for high 
throughput performance. 
5 a. Expression in bacterial hosts 

The mutated forms of the nucleic acid are prepared or introduced 
into plasmids for expression in bacterial cells. The genetic variants are 
expressed from suitable bacterial cells, which are prepared by 
transformation of aliquots of the cells with each member of the plasmid 
10 library (each genetic variant continues to be physically separated from 
each other). 

b. Expression in eukaryotic host cells 
A virus library is generated from the plasmid library. The virus 
library, in which each different member is separately maintained, is 
15 prepared by: 

(1) Transaction of the plasmid library into 
appropriate virus-producer cells (viruses produced, each one carrying a 
different genetic variant present in the original plasmid library, are 
physically separated from each other); 

20 (2) Titration of the virus library (of each individual 

virus present in the library, separately). Titration is effected by any 
method, but generally by either a method designated Real Time Virus 
Titering (RTVT™) (see, International PCT application No. 
PCT/FR01/01366 published as International PCT application No. 

25 WO 01/186291 and the EXAMPLES below) or a refinement of that 
method provided herein and designated Tagged Replication and 
expression enhancement (TREE™) described above and in the examples; 
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{3) Standardization of the virus library to equal 
concentrations of all the individual viruses in the library (individual viruses 
continue to be physically separated from each other); 

(4) Expression of the genetic variants from 
5 appropriate mammalian cells by transduction with the virus library (each 
genetic variant continues to be physically separated from each other and 
each individual virus is handled separately from the others). 

3. Phenotypic characterization of the variant proteins. 
The variant proteins are expressed (from either plasmids in bacterial 
10 cells (step 2) or viruses in mammalian cells (step 4)) and their activity is 
assessed in one or more appropriate specific assays. The assays can be 
both types: biochemical (cell-free) assays and/or cell-based assays. The 
variant proteins in the library are physically separated from each other and 
their activity is individually assessed on a one-by-one basis. 
15 The assays can be performed in one of a variety of ways, 

including, but are not limited to: 

a. Using a single-point dilution for each individual variant 
protein, followed by a kinetic analysis (multiple time points) of the read- 
out by technologies like Tagged Replication and expression enhancement 

20 (TREE™), or any other appropriate technology 

b. Using serial dilutions of each individual variant protein, 
followed by, for example, the Hill-based analysis of the read-out by 
technologies or any other appropriate technology. Hill based analyses 
assess the interaction between cells and biological agents (see, published 

25 International PCT application No. WO 01/44809 based on PCT n° 
PCT/FR0O/O3503). 

The goal of these methods is to identify proteins having an evolved 
function or property. 
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Lead identification 
Based on the results obtained from the assays described above, 
each individual protein variant is individually tested for the parameters 
that assess the activity, property, function or structure of interest. 
5 Variants are ranked out according to their activity features. Those variant 
proteins best suited for the specificities of each individual project and 
system under study are then selected. The selected leads can be used for 
the desired purpose or further evolved or mutated to achieve desired 
activities. 

10 Typically, as for most directed evolution methods, the process is an 

iterative one, in which mutated variants are produced, screened, the best 
identified and then selected. The selected variants are then subjected to 
further evolution and the screening process repeated. This is repeated 
until the desired goal is achieved. 

15 This further evolution may employ the methods herein or any 

directed evolution method or combinations thereof. The methods for 
variant production will include the amino-acid scan method herein, which 
provides a rational approach to variant generation. Other rounds can 
include combinations of any other method for directed evolution known 

20 and/or combinations thereof. 

D. Directed evolution of a viral gene 

Recombinant viruses have been developed for use as gene therapy 
vectors. Gene therapy applications are hampered by the need for 
development of vectors with traits optimized for this application. The 

25 high throughput methods provided herein are ideally suited for 

development of such vectors. In addition to use for development of 
recombinant viral vectors for gene therapy, these methods can also be 
used to study and modify the viral vector backbone architecture, trans- 
complementing helper functions, where appropriate, regulatable and 
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tissue specific promoters and transgene and genomic sequence analyses. 
Recombinant AAV (rAAV) is a gene therapy vector that can serve as a 
model for application of the methods herein for these and other purposes. 
Adeno-associated virus (AAV) is a defective and non-pathogenic 
5 parvovirus that requires co-infection with either adenovirus or herpes 
virus, which provide helper functions, for its growth and multiplication. 
There is an extensive body of knowledge regarding AAV biology and 
genetics (see, e.g., Weitzman eta/. (1996) J. Virol. 70: 2240-2248 
(1996); Walker eta/. (1997) J. Virol. 77:2722-2730; Urabe et al. (1999) 

10 J. Virol. 23:2682-2693; Davis etaf. (2000) J. Virol. 23:74:2936-2942; 
Yoon etal. (2001) J. Virol. 75:3230-3239; Deng et al. (1992) Anal 
Biochem 200:81-85; Drittanti etal. (2000) Gene Therapy 7:924-929; 
Srivastava etal. (1983) J. Virol. 45:555-564; Hermonat etal. (1984) J. 
Virol. 57:329-339; Chejanovsky et al. (1989) Virology 773:120-128; 

1 5 Chejanovsky et al. ( 1 990) J. Virol. 64: 1 764-1 770; Owens et al. (1991) 
Virology 754:14-22; Owens era/. (1992) J. Virol. 55:1236-1240; 
Qicheng Yang etal. (1992) J. Virol. 55:6058-6069; Qicheng Yang etal. 
(1993) J. Virol. 57:4442-4447; Owens etal. (1993) J. Virol. 52:997- 
1005; Sirkka etal. (1994) J. Virol. 55:2947-2957; Ramesh etal. (1995) 

20 Biochem. Biophy. Res. Com. Vol 210 (3), 717-725; Sirkka (1995) J. 
Virol. 55:6787-6796; Sirkka etal. (1996) Biochem. Biophy. Res. Com. 
220:294-299; Ryan etal. (1996) J. Virol. 70:1542-1553; Weitzman etal. 
(1996) J. Virol. 70:2440-2448; Walker etal. (1997) J. Virol. 77:2722- 
2730; Walker etal. (1997) J. Virol. 77:6996-7004; Davis etal. (1999) J. 

25 Virol. 73:2084-2093; Urabe etal. (1999) J. Virol. 73:2682-2693; Gavin 
etal. (1999) J. Virol. 73:9433-9445; Davis etal. (2000) J. Virol. 
74:2936-2942; Pei Wu etal. (2000) J. Virol. 74:8635-8647; Alessandro 
Marcello et al. (2000) J. Virol. 74:9090-9098). AAV are members of the 
family Parvoviridae and are assigned to the genus Dependovirus . 
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Members of this genus are small, non-enveloped, icosahedral with linear 
and single-stranded DNA genomes, and have been isolated from many 
species ranging from insects to humans. 

AAV can either remain latent after integration into host chromatin 
5 or replicate following infection. Without co-infection, AAV can enter host 
cells and preferentially integrate at a specific site on the q arm of 
chromosome 1 9 in the human genome. 

The AAV genome contains 4975 nucleotides and the coding 
sequence is flanked by two inverted terminal repeats (ITRs) on either side 

10 that are the only sequences in cis required for viral assembly and 

replication. The ITRs contain palindromic sequences, which form a hairpin 
secondary structure, containing the viral origins of replication. The ITRs 
are organized in three segments: the Rep binding site (RBS), the terminal 
resolution site (TRS), and a spacer region separating the RBS from the 

15 TRS. 

Regulation of AAV genes is complex and involves positive and 
negative regulation of viral transcription. For example, the regulatory 
. proteins Rep 78 and Rep 68 interact with viral promoters to establish a 
feedback loop (Beaton et aL (1989) J. Virol 53:4450-4454; Hermonat 

20 (1994) Cancer Lett 81: 129-1 36). Expression from the p5 and p19 

promoters is negatively regulated in trans by these proteins. Rep 78 and 
68, which are required for this regulation, have bind to inverted terminal 
repeats (ITRs; Ashktorab eta/. (1989) J. Virol. 63:3034-3039) in a site- 
and strand-specific manner, in vivo and in vitro. This binding to ITRs 

25 induces a cleavage at the TRS and permits the replication of the hairpin 
structure, thus, illustrating the Rep helicase and endonuclease activities 
(Im et aL (1990) Cell 61:447-457; and Walker et aL (1997) J. Virol. 
77:6996-7004), and the role of these non-structural proteins in the initial 
steps of DNA replication (Hermonat et aL (1984) J. Virol. 52:329-339). 
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Rep 52 and 40, the two minor forms of the Rep proteins, do not bind to 
ITRs and are dispensable for viral DNA replication and site-specific 
integration (Im et aL (1 992) J. Virol. 66 A 1 1 9-1 1 2834; Ni et at. (1 994) J. 
Virol. 53:1128-1138. 
5 The genome {see, FIG. 4) is organized into two open reading 

frames (ORFs, designated left and right) that encode structural capsid 
proteins (Cap) and non-structural proteins (Rep). There are three 
promoters: p5 (from nucleotides 255 to 261: TATTTAA), p19 (from 
nucleotide 843 to 849: TATTTAA) and p40 (from nucleotides 1822 to 

10 1827: ATATAA). The right-side ORF (see FIG. 4) encodes three capsid 
structural proteins (Vp 1-3). These three proteins, which are encoded by 
overlapping DNA, result from differential splicing and the use of an 
unusual initiator codon (Cassinoti eta/. (1988) Virology 767:176-184). 
Expression of the capsid genes is regulated by the p40 promoter. Capsid 

15 proteins VP1, VP2 and VP3 initiate from the p40 promoter. VP1 uses an 
alternate splice acceptor at nucleotide 2201; whereas VP2 and VP3 are 
derived from the same transcription unit, but VP2 use an ACG triplet as 
an initiation codon upstream from the start of VP3. On the left side of 
the genome, two promoters p5 and p19 direct expression of four 

20 regulatory proteins. The left flanking sequence also uses a differential 
splicing mechanism (Mendelson et aL (1986) J. Virol 50:823-832) to 
encode the Rep proteins, designated Rep 78, 68, 52 and 40 on the basis 
molecular weight. Rep 78 and 68 are translated from a transcript 
produced from the p5 promoter and are produced from the unspliced and 

25 spliced form, respectively, of the transcript. Rep 52 and 40 are the 
translation products of unspliced and spliced transcripts from the p19 
promoter. 

The rep protein is a adeno-associated virus protein involved in a 
number of biological processes necessary to AAV replication. The 
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production of the rRep proteins enables viral DNA to replicate, 
encapsulate and integrate (McCarty et al. (1992) J. Virol 55:4050-4057; 
Horer et al. (1995) J. Virol 53:5485-5496, Berns et al. (1996) Biology of 
Adeno-associated virus, in Adeno-associated virus (AAV) Vectors in Gene 
5 Therapy, KJ. Berns and C. Giraud, Springer (1996); and Chiorini et al. 
(1996) The Roles of AAV Rep Proteins in gene Expression and Targeted 
Integration, from Adeno-associated virus (AAV) Vectors in Gene Therapy, 
K.I. Berns and C. Giraud, Springer (1996)). A rep protein with improved 
activity could lead to increased amounts of virus progeny thus allowing 

10 higher productivity of rAAV vectors. 

AAV and rAAV have many applications, including use as a gene 
transfer vector, for introducing heterologous nucleic acid into cells and for 
genetic therapy. Advances in the production of high-titer rAAV stocks to 
the transition to human clinical trials have been made, but improvement of 

15 rAAV production will be complemented with special attention to clinical 
applications of rAAV vectors as a successful gene therapy approach. 
Productivity of rAAV (i.e. the amount of vector particles that can be 
obtained per unitary manufacturing operation) is one of the rate limiting 
steps in the further development of rAAV as gene therapy vector. 

20 Methods for high throughput production and screening of rAAV have 
been developed (see, e.g., Drittanti et al. (2000) Gene Therapy 7:924- 
929). Briefly, as with the other steps in methods provided herein, the 
plasmid preparation, transfection, virus productivity and titer and 
biological activity assessment are intended to be performed in an 

25 automatable high throughput format, such as in a 96 well (or other 
number or multiples thereof, such as 384, 1536 . . . 9600, 9984 . . •) 
formats. 

Since the Rep protein is involved in replication it can serve as a 
target for increasing viral production. Since it has a variety of functions 
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and its role in replication is complex, it has heretofore been difficult to 
identify mutations that result in increase viral production. The methods 
herein, which rely on in vivo screening methods, permit optimization of its 
activities as assessed by increases in viral production. Provided herein 
5 are R e p proteins and viruses and viral vectors containing the mutated Rep 
proteins that provide such increase. The amino acid positions on the rep 
proteins that are relevant for rep proteins activities in terms of AAV or 
rAAV virus production are provided. Those amino acid positions are such 
that a change in the amino acid leads to a change in protein activity either 

•10 to lower activity or increase activity. As shown herein, the alanine or 
amino acid scan revealed the amino acid positions important for such 
activity (i.e. hits). Subsequent mutations produced by systematically 
replacing the amino acids at the hit positions with the remaining 18 amino 
acids produced so-called "leads" that have amino acid changes and result 

15 in higher virus production. In this particular example, the method used 
included the following specific steps. 
Amino acid scan 

In order to first identify those amino acid (aa) positions on the rep 
protein that are involved in rep protein activity, an Ala-scan was 

20 performed on the rep sequence. For this, each aa in the rep protein 

sequence was individually changed to Alanine. Each resulting mutant rep 
protein was then expressed and the amount of virus it could produced 
measured as indicated below. The relative activity of each individual 
mutant compared to the native protein is indicated in FIG 3A. HITS are 

25 those mutants that produce a decrease in the activity of the protein (in 
the example: all the mutants with activities below about 20 % of the 
native activity). 

In a second experimental round, which included a new set of 
mutations and phenotypic analysis, each amino acid position hit by the 
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Ala-scan step, was mutated by amino acid replacement of the native 
amino acid by the remaining 18 amino acids, using site-directed 
mutagenesis. 

In both rounds, each mutant was individually designed, generated 
5 and processed separately, and optionally in parallel with the other 
mutants. Neither combinatorial generation of mutants nor mixtures 
thereof were used in any step of the method. 

A plasmid library was thus generated in which each plasmid contained 
a different mutant bearing a different amino acid at a different hit 
10 position. Again, each resulting mutant rep protein was then expressed 
and the amount of virus it could produce measured as indicated below. 
The relative activity of each individual mutant compared to the native 
protein is indicated in FIGURE 3B. LEADS are those mutants that lead to 
an increase in the activity of the protein (in the example: the ten mutants 
15 with activities higher, typically between 6 to 10 times, than the native 
activity) . 

Expression of the genetic variants and phenotypic characterization. 
The rep protein acts as an intracellular protein through complex 
interaction with a molecular network composed by cellular proteins, DNA, 

20 AAV proteins and adenoviral proteins (note: some adenovirus proteins 
have to be present for the rep protein to work). The final outcome of the 
rep protein activity is the virus offspring composed by infectious rAAV 
particles. It can be expected that the activity of rep mutants would affect 
the titer of the rAAV virus coming out of the cells. 

25 As the phenotypic characterization of the rep variants can only be 

accomplished by assaying its activity from inside mammalian cells, a 
mammalian cell-based expression system as well as a mammalian cell- 
based assay was used. The individual rep protein variants were expressed 
in human 293 HEK cells, by transfection of the individual plasmids 
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constituting the diverse plasmid library. All necessary functions were 
provided as follows: 

(a) the cellular proteins present in the permissive specific 293 HEK 

cells; 

5 (b) the AAV necessary proteins and DNA were provided by co- 

transfection of the AAV cap gene as well as a rAAV plasmid vector 
providing the necessary signaling and substrate ITRs sequences; 

(c) the adenovirus (AV) proteins were provided by co-transfection 
with a plasmid expressing all the AV helper functions. 

10 A library of recombinant viruses with mutant rep encoding genes 

was generated. Each recombinant, upon introduction into a mammalian 
cell and expression resulted in production of rAAV infectious particles. 
The number of infectious particles produced by each recombinant was 
determined in order to assess the activity of the rep variant that had 

15 generated that amount of infectious particles. 

The number of infectious particles produced was determined in a 
cell-based assay in which the activity of a reporter gene, in the 
exemplified embodiment, the bacterial lacZ gene, or virus replication (Real 
time PCR) was performed to quantitatively assess the number of viruses. 

20 The limiting dilution (titer) for each virus preparation (each coming from a 
different rep variant) was determined by serial dilution of the viruses 
produced, followed by infection of appropriate cells (293 HEK or HeLa 
rep/cap 32 cells) with each dilution for each virus and then by 
measurement of the activity of the reporter gene for each dilution of each 

25 virus. Hill plots (NAUTSCAN") as described herein (published as 
International PCT application No. WO 01/44809 based on PCT n° . 
PCT/FR00/03503, Dec, 2000; see EXAMPLES) or a second order 
polynomial function (Drittanti et al. (2000) Gene Ther. 7: 924-929) was 
used to analyze the readout data and to calculate the virus titers. Briefly, 
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the titer was calculated from the second order polynomial function by 
non-linear regression fitting of the experimental data. The point where 
the polynomial curve reaches its minimum is considered to be the titer of 
the rAAV preparation. A computer program for calculation of titers has 
5 been developed {see Drittanti et al. (2000) Gene Then 7: 924-929) to 
assess the minima. 

The TREE method described herein can be used to analyze the 
readout data and to calculate the virus titers. The results are shown in the 
EXAMPLE below. 

10 Comparison between results of full-length Hit position analysis 

reporter here and the literature 

The experiments identified a number of heretofore unknown 
mutation loci, which include the hits at positions: 4, 20, 22, 28, 32, 38, 
39, 54, 59, 124, 125, 127, 132, 140, 161, 163, 193, 196, 197, 221, 

15 228, 231, 234, 258, 260, 263, 264, 334, 335, 341, 342, 347, 350, 
354, 363, 364, 367, 370, 376, 381, 389, 407, 41 1, 414, 420, 421, 
422, 428, 429, 438, 440, 451, 460, 462, 484, 488, 495, 497, 498, 
499, 503, 511, 512, 516, 517 and 518 with reference to the amino 
acids in Rep78 and Rep 68. Rep 78 is encoded by nucleotides 321- 

20 2,186; Rep 68 is encoded by nucleotides 321-1 906 and 2228-2252; Rep 
52 is encoded by nucleotides 993-2186, and Rep 40 is encoded by amino 
acids 993-1906 and 2228-2252 of wildtype AAV. 

Also among these are mutations that may have multiple effects. 
Since the Rep coding region is quite complex, some of the mutations may 

25 have several effects. Amino acids 542, 598, 600 and 601, which are in 
the Rep 68 and 40 intron region, are also in the coding region of Rep 78 
and 52. Codon 630 is in the coding region of Rep 68 and 40 and non 
coding region of Rep 78 and 52. 
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Mutations at 10, 86, 101, 334 and 519 have been previously 
identified, and mutations, at loci 64, 74, 88, 175, 237, 250 and 429, but 
with different amino acid substitutions, have been previously reported. In 
all instances, however, the known mutations reportedly decrease the 
5 activity of Rep proteins. Among mutations described herein, are 
mutations that result in increases in the activity the Rep function as 
assessed by detecting increased AAV production. 
Lead identification. 

Based on the results obtained from the assays described herein (i.e. 
10 titer of virus produced by each rep variant), each individual rep variant 
was assigned a specific activity. Those variant proteins displaying the 
highest titers were selected as leads and are used to produce rAAV. 

In further steps, rAAV and Rep proteins that contain a plurality of 
mutations based on the hits (see Table in the EXAMPLES, listing the hits 
. 15 and lead sites), are produced to produce rAAV and Rep proteins that 
have activity that is further optimized. Examples of such proteins and 
AAV containing such proteins are described in the EXAMPLES. 

The rAAV rep mutants are used as expression vectors, which, for 
example, can be used transiently for the production of recombinant AAV 
20 stocks. Alternatively, the recombinant plasmids may be used to generate 
stable packaging cell lines. To create a stable producer cell line, the 
recombinant vectors expressing the AAV with mutant rep genes, for 
example, are cotransfected into host cells with a plasmid expressing the 
neomycin phosphotransferase gene (neor) by transfection methods well 
25 known to those skilled in the art, followed by selection for G418 
resistance. 

Also among the uses of rAAV, particularly the high titer stocks 
produced herein, is gene therapy for the purpose of transferring genetic 
information into appropriate host cells for the management and correction 
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of human diseases including inherited and acquired disorders such as 
cancer and AIDS. The rAAV can be administered to a patient at 
therapeutically effective doses. A therapeutically effective dose refers to 
that amount of the compound sufficient to result in amelioration of 
5 symptoms of disease. 

Gene therapy 

Toxicity and therapeutic efficacy of the rAAV can be determined by 
standard pharmaceutical procedures in cell cultures or experimental 
animals, e.g., for determining the LDS 50 (the dose lethal to 50% of the 

10 population) and the ED50 (the dose therapeutically effective in 50% of 
the population). The dose ratio between toxic and therapeutic effects is 
the therapeutic index and it can be expressed as the ratio LD 50 /ED 5O . 
Doses that exhibit large therapeutic indices are preferred. Doses that 
exhibit toxic side effects may be used; care should be taken to design a 

1 5 delivery system that targets rAAV to the site of treatment in order to 
minimize damage to untreated ^ells and reduce side effects. 

The data obtained from cell culture assays and animal studies can 
be used in formulating a range of dosage for use in humans. The dosage 
of such rAAV lies preferably within a range of circulating concentrations 

20 that include the ED 50 with little or no toxicity. The dosage may vary 
within this range depending upon the dosage form employed and the 
route of administration utilized. A therapeutically effective dose can be 
estimated initially from cell culture assays. A dose may be formulated in 
animal models to achieve a circulating plasma concentration range that 

25 includes the IC 50 (ie., the concentration of the test compound which 
achieves a half-maximal infection or a half- maximal inhibition) as 
determined in cell culture. Such information can be used to more 
accurately determine useful doses in humans. Levels in plasma may be 
measured, for example, by high performance liquid chromatography. 
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The following examples are included for illustrative purposes only 

and are not intended to limit the scope of the invention. The specific 

methods exemplified can be practiced with other species. The examples 

are intended to exemplify generic processes. 

5 EXAMPLE 1 

Titering or assessment of concentration by a method designated Real 
Time Vector Titering (RTVT™) 

This Example is based on the method described in International PCT 
application No. PCT/FR01/01 366, based on French application 

10 n° 0005852, filed 9 May 2000, and published as International PCT 
application No. WO 01/186291. This method assesses the 
titer or concentration of a biological agent (virus, gene transfer vector) in 
a sample, by measurement of the kinetics of change of a reporter 
parameter following the exposure of cells to the biological agent. 

15 As noted above, reporter parameters may include, but are not 

limited to,: gene / transgene expression related to the gene/transgene 
products, such as enzymatic activity, fluorescence, luminescence, antigen 
activity, binding to receptors or antibodies, and regulation of gene 
expression), differential gene expression, viral/vector progeny 

20 productivity, toxicity, cytotoxicity, cell proliferation and/or differentiation 
activity, anti-viral activity, morphogenetic activity, pathogenetic activity, 
therapeutic activity, tumor suppressor activity, oncogenetic activity, 
pharmacological activity. 

Serial dilution methods 

25 The assessment of the concentration or titer of biological agents 

using current approaches needs serial dilutions of the agent. 
Serial dilutions of the agent are applied to a cell-based reporter system, 
that elicits an output signal in response to the exposure to the agent. The 
intensity of the signal is a function of the concentration of the agent. The 
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titer or concentration of the agent is determined as the highest dilution 
that still elicits a measurable response in the output. The higher the 
number of dilutions tested, the higher the accuracy of the value obtained 
for the titer. 

5 This approach requires a set of serial dilutions for every biological 

agent whose titer needs to be determined. Thus, the application of this 

approach to the simultaneous titration of a large number of different 

biological agents is limited by the number of experimental points needed 

(example: for 30 biological samples: 20 serial dilutions x 30 biological 

10 agents: 600 experimental points). 

The approach in International PCT application No. 
PCT/FR0 1/0 1366 published as International PCT application 
No. WO 01/186291 

The intensity of the output signal (after exposure of reporter cells 
15 to the biological agent) is not only dependent on the concentration of the 
agent but also on the time after exposure. As time increases, the intensity 
of the signal increases. The kinetics of change of the intensity over time 
depends upon the concentration of the agent. Thus, lower concentrations 
of the agent will require longer times for the intensity to reach a given 
20 value that would be reached in shorter times after exposure to higher 
concentrations of the same agent. 

This approach (designated Real Time Virus Titering (RTVT"*) uses 
the following: a reference plot representing the relationship between the 
concentration of the agent and the time necessary for the intensity to 
25 reach a given threshold value is obtained using a reference preparation of 
biological agent, whose concentration or titer is known. This plot is then 
used to obtain the concentration of the biological agent under study by 
entering the time that a dilution of that agent needed for the intensity to 
reach the threshold value. 
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Using this approach, there is no need for serial dilutions of the 
biological agent(s) under study. Once the reference plot (tfi vs c) is 
obtained, it can be used for the determination of the concentration or titer 
of any number of biological agents. Only one dilution of the biological 
5 agent under study is necessary to obtain the corresponding value of tfi 
that is then used to obtain the concentration or titer using the reference 
plot. 

Thus, the application of this approach to the simultaneous titration 
of large numbers of different biological agents is facilitated by the fact 
10 that only one dilution of each sample is needed (example: for 30 

biological samples: 7 dilution x 30 biological agents: 30 experimental 
points (compared to 600 needed with the current approach). 
This approach is specially suited for the high throughput assessment of 
concentration or titer of large numbers of biological agents. 
1 5 The system 

The system includes the following elements: 
a preparation of the biological agent (virus, gene transfer vector, 
protein,...) whose concentration or titer is unknown and has to be 
determined. 

20 a reporter system including culture of a cell line (or a mixture 

of cell lines) that reacts to the exposure to the biological 
agent by displaying a specific output signal, 
a master preparation of a reference biological agent, of 
known concentration or titer, that is able to generate the 

25 output signal when the reporter cells are exposed to it. 

Practice of the method 

When the reporter cells are exposed to either the biological agent 
under study or the reference biological agent, an output signal is 
generated, that can be measured. 
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The intensity of the output signal is called /; the concentration of the 
biological agent used is called c, the time of exposure of the cells to the 
biological agent is called t. The intensity of the output signal (/) is a 
function of c and t : 
5 / increases as the concentration (c ) of the biological agent applied 

to the cells increases; 

/ increases as the time of exposure of the cells to the biological 
agent (f) increases. 

If the time t after exposure of the cells to the biological agent is 
10 kept constant, then, / will change as a direct function of c . 

If the concentration c of the biological agent is kept constant, 
then, / will be a direct function of f. 

P is defined here as a threshold value of the intensity of the output 
signal, arbitrarily defined for every system under study. 
15 XJ3 is defined here as the time necessary to reach the threshold 0. 

Use of fi and t0 to determine the concentration or titer 
of a biological agent. 

The reporter cells are exposed to serial dilutions of a reference 
biological agent, whose concentration (or titer) is previously known. The 
20 intensity of the output signal (/) is measured at several time points (t) for 
every concentration (serial dilutions) c of the reference biological agent. 

/ is plotted vs t , and that, for every concentration c used of the 
reference biological agent. 

Using the plots obtained above, and for every concentration c of 
25 the reference biological agent, the time {t/3) necessary for the intensity of 
the output signal to reach a threshold value JS is obtained. 

With the data obtained above, xp is plotted vs c. 

This plot represents the time necessary for the intensity of the 
output signal to. reach the threshold value jJasa function of the 
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concentration of the biological agent used. This is a standard plot and 
will be used to determine the unknown concentration of the biological 
agent under study by measuring the time that a given dilution of it needs 
to give an output signal whose intensity equals the threshold /?. 
5 The reporter cells are exposed to a dilution of the biological agent 

under study (whose concentration or titer is to be determined). The 
intensity of the output signal (/) is measured over time until it reaches the 
threshold value p. The time necessary for / to reach the value fi is 
recorded as ty?. 

10 The tfi value recorded above is entered into the standard plot 

obtained above and the corresponding c value is obtained. 

This c value represents the concentration or titer of the biological 

agent under study. 

Example of the Real Time Virus Titering RTVT" 
15 titering method 

Rat-2 cells were infected with serial dilutions of a reference 

preparation of a retroviral vector carrying the green fluorescent protein 

(GFP) gene (vector pSI-EGFPI see, Ropp eta/. (1995) Cytometry 27:309- 

317). At increasing times after infection, the level of expression of the 

20 transgene was determined (as the level of fluorescence due to the GFP 

gene) as the output signal. 

Table 3 represents the values obtained: 


Concentration 

Time after infection 

Output signal 

(1 = 10 6 particles/ml) 

(hrsl 

fluorescence 

0.1 

16 

20.4 


24 

30.1 


40 

95.1 


48 

138.7 


64 

157.3 

0.25 

16 

26.8 


24 

48.5 


40 

173.3 
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48 

228.2 


64 

191.7 

0.5 

16 

38.1 


24 

72 


40 

198.7 


48 

296.2 


64 

203.7 


The threshold value of £ = 100 was arbitrarily selected for this example. 
10 The time {tfi) necessary for the output signal to reach the threshold P , for 
every concentration is shown in table 4. 
Table 4 


15 


Concentration 
(1 - 10 6 particles/ml) 

\ (hrs) 

0.1 

42 

0.25 

31 


A plot of tfi versus concentration for the reference virus shows that 
the concentration and Xfi exhibit a clearly defined relationship, that allows 
20 for the calculation of the concentration (c) of a sample, if the 
corresponding tfi of that sample is known. 

EXAMPLE 2 

Tagged Replication and Expression Enhancement (TREE) for titering 
As discussed above, TREE is a method for titering and 

25 standardization of preparations of viruses, vectors, antibodies, libraries, 
proteins, genes and any other moiety that is detectable based upon a 
output signal, such as fluorescence. The TREE method is an improvement 
of the Real Time Virus Titering (RTVT) method (see, International PCT 
application No. PCT/FR01/01 366 published as International PCT 

30 application No. WO 01/186291). It is performed with a reporter moiety, 
such as a reporter virus (with a known titer) and the test sample (with 
unknown titer). The reporter, such as a reporter virus, has a readily 
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detectable output signal that can be measured as a function of time. The 
effect of the moiety, such as a virus, of unknown titer is assessed. The 
moiety whose titer is assessed either increases or decreases the output 
signal as a function of time. This change in signal is used to assess the 
5 amount or concentration of the moiety of unknown concentration, and 
hence its titer. 

The method is exemplified herein using an AAV system for the 
determination of the titer of an AAV vector and an AAV-reporter vector 
as a competitor and wild type Adenovirus as helper virus. One of skill in 

10 the art readily can adapt the method to other systems, including other 

viruses, and other moieties for which a reporter system can be developed. 
Other such moieties include, but are not limited to, viral vectors, 
plasmids, libraries, proteins, antibodies, vaccines, genes, and nucleic acid 
molecules. 

1 5 Materials and Methods 

1 . Cells and Viruses 
HeLa rep-cap32 cells, a HeLa derived cell line (kindly provided by P. 
Moullier, Laboratoire de Therapie Genique, CHU, Nantes; see, Salvetti et 
a/. (1998) Hum Gene Ther 20:695-706; Chadeuf eta/. (2000) J Gene 

20 Med 2:260-268) was grown in DMEM with 10% fetal calf serum. These 
cells were plated 24 h before infection at a density of 1 x 10 4 cells in 
single well of 96-well plates. rAAV-LacZ (10 10 ip/ml), rAAV-GFP (10 9 
ip/ml) vectors and Human Adenovirus type 5 (Ad5) (10 11 pfu/ml) were 
from CHU, Nantes. 

25 Hela rep-cap32 cells had been produced by cotransfecting plasmid 

pspRC, which harbors the AAV rep-cap genome with the ITRs deleted (bp 
190 to 4484 of wild-type AAV), with plasmid PGK-Neo, conferring 
resistance to G41 8 on Hela cells (see, Chadeuf et aL (2000) J. Gene 
Med. 2:260-268 and Salvetti eta/. (1998( Hum Gene Then 5:695-706). 
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Hela rep-cap 32 cells are a packaging line that harbor one copy of the 
genome with the ITRs deleted (see, also Tessier eta/. (2001) J. Virol. 
75/375-383). 

Plasmid pspRC contains the AAV genome (positions 190-4,484 bp) 
5 with the ITRs deleted and was obtained by excising the rep-cap fragment 
(Xba\ fragment) from the well-known vector psub201 (Samulski et aL 
(1987) J Virol 57:3096-31 01; also called pSSV9) by Xbal digestion and 
inserting it into the Xbal site of plasmid pSP72 (Promega). Plasmid 
psub201 (see, e.g., U.S. Patent No. 5,753,500) is a modified full-length 

10 AAV type 2 genomic clone contains all of the AAV type 2 wild-type 
coding regions and cis acting terminal repeats. 

2. Infection and measurement 
Four serial dilutions of a rAAV-LacZ (0.01, 0.0075, 0.005 and 
0.0025 //I, see Table 2 below, designated samples 1-4, respectively) were 

15 made and used for co-infection of HeLa rep-cap32 cells together with 8 
different Ad5 multiplicity of infection (MOI; from 0.1 to 100/cell) and with 
10" 3 ml (10 8 infectious particles (ip)) or 10* 4 ml (10 5 ip) rAAV-QFP viral 
vector. All the samples were done in triplicate. After infection, the plates 
were read at different times, from 34.5 h to 80 h (every 30 minutes). 

20 rAAV-GFP is an SSV9-derived vector; SSV9 is a clone containing 

the entire adeno-associated virus (AAV) genome inserted into the Pvull 
site of plasmid pEMBL (see, Du et al. (1996) Gene Ther 3:254-261). The 
rAAV-GFP and rAAV-LacZ plasmids are SSV9 with a GFP or LacZ gene 
under control of the cytomegalovirus (CMV) immediate-early promoter. 

25 All the samples were done in triplicate. 
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3. Process 

Figure 2 shows the overall procedure in 96 well format. Cells were 
plated 24 h before infection. Co-infection of rAAV-GFP with serial 
dilutions of rAAV-LacZ together with Ad5 (different MOD, were done. 
5 Then the plates were read at different times using the Analyst AD&HT 
micro plate-reader (LJL BioSystems). 

4. Analysis 

For this kinetic technique, Fluorescence Intensity (Fl) of the infected 
cells is measured as a function of the time. Serial dilutions of the AAV- 

10 competitor vector AAV-lacZ vector, which decreases the fluorescence 
signal, are performed. For this example, fluorescence was measured for 
AAV-GFP with 10 6 ip and 10 5 ip and then 10 s ip of the AAV-GFP reporter 
virus in the serial dilutions of the competitor virus, AAV-lacZ vector in a 
96-well format (samples 1-4, see Table below). Measurements were 

15 taken of each well and curves of Fl (of the GFP) versus time (hrs) were 
obtained (see Figure 2B). 

An arbitrary one value for Fl (see Fig. 2B, 6 x 10 6 Fl units), 
typically, though not necessarily, near the greatest separation among the 
curves so that the numbers are readily discernable, was selected. The 

20 point at which each of the curves intersect this value is beta time (tfi) for 
that combination of amounts of reporter plus dilution of the virus of 
unknown titer, tff, taken from the Fl vs. Time (hrs) curves, for each 
sample containing a dilution of the unknown plus 1 0 6 ip of the reporter 
virus is set forth in column 3 of Table 2 below. 

25 To determine the titer of the test virus, the Xp for the AAV GFP 

(reporter virus) is plotted versus quantity of ip (i.e a straight line between 
the X0 for the 10 6 ip and the 10 s ip) (Fig. 2C). For any tfi of the unknown 
virus, the quantity of ip can be determined from this curve. The beta 
time (tP) of each sample (in this case for the different dilutions of rAAV 
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10 


15 


20 


25 


LacZ mixed with 10 6 infectious particles of rAAV-GFP) is determined, and 
then the residual number of infectious particles of rAAV-GFP for each 
sample. The difference between 10 6 ip of rAAV-GFP put in each sample 
and the number of ip detected by fluorescence in the same well is the 
actual quantity of rAAV-GFP competed (consumed) by the unknown 
rAAV (in this case rAAV-LacZ). This number is determined for each 
dilution. The quantity rAAV-GFP consumed is the same quantity of 
unknown rAAV in the sample. This quantity is present in one volume of 
unknown rAAV, which in this example is 1 ml. Based upon this, the 
infectious titer of the unknown rAAV is determined. The results are 
shown in Table 2. 

TABLE 2 

AAV LacZ titration by TREE titration 


Sample 

volume 

Xfi (hrs) 

Residual 
AAV-GFP 

Consumed 
AAV-GFP 

AAV LacZ 
Concentration 
{i.p./10 2 //I) 

1 

10.0 x 10* 3 

66.5 

5.56 x 10 s 

4.44 x 10 s 

4.44 x 10* 

2 

7.5 x 10" 3 

66.5 

. 5.56 x 10 s 

4.44 x 10 5 

5.93 x 10 5 

3 

5.0 x 10' 3 

65 

7.06 x 10 s 

2.94 x 10 s 

5.88 x 10 s 

4 

2.5 x 10' 3 

64 

.8.06 x 10 5 

1.94 x 10 s 

7.76 x 10 s 


30 


ip//y| = 6.01 x 10 5 ip/0.01 /yl). The standard deviation was 1.37 x 10 
ip/ml with an error of ±23%. 

EXAMPLE 3 

Hill Analysis of the screening assay output 

It is important to have reliable methods for screening and/or 
evaluating the performance of a set of biological agents, such as a library 
of viral or non-viral recombinant vectors, vaccines, recombinant proteins 
and antibodies, in a complex biological system, such as living target cells 
When developing such agents, for example gene therapy vectors and 
other agents for therapeutic use, it is necessary to be able to evaluate 
and compare performance among candidates. 


WO 03/023032 


PCT/IB02/03921 


-70- 

The progress of gene transfer into gene therapy depends upon the 
capacity to develop gene transfer vectors into therapeutic drugs. Clinically 
relevant vectors need to be efficient and safe, in reaching and infecting 
target cells and in ensuring a persistent level of expression of the 
5 therapeutic gene with a minimum of adverse effects. The availability of 
standardized quantitative methods, suitable for an accurate and objective 
assessment of titer, performance and safety, is necessary for the 
pharmaceutical development of gene vectors as drugs. 

Any method for assessment is contemplated herein as long as it is 
10 adapted for use in a high throughput format. Of particular interest is the 
Hill equation based method of international PCT application No. 
WO 01/44809 (International PCT application No. PCT/FR00/03503, based 
on French application FR 9915884). 

Two widely used parameters that provide quantitative information 
15 about the potential performance of a gene transfer vector preparation are 
the titer of physical particles and the titer of infectious particles. Vector 
preparations with high titer of infectious particles and low physical 
particles/infectious particles ratio are considered to be of higher quality. 

The titer in physical particles(pp) (see, e.g., Mittereder et at. (1996) 
20 J. Virol. 70:7498-7509; Atkinson eta/. (1998) Nucl. Acids Res. 26:2821- 
2823; Kechli etaL (1998) Hum. Gene Ther. 3:587-590; and Nelson et al. 
(1998) Hum. Gene Ther. 3:2401-2405), which represents the total 
number of vector particles, is usually evaluated from the vector content 
by detecting the nucleic acid contents (nucleic acids hybridization and 
25 OD 260 respectively for AAV and AdV), detecting viral protein content (for 
example, reverse transcriptase (RT) activity and p24 content for MLV and 
HIV, respectively). 

Among the physical particles (pp), there are particles potentially 
active in performing transduction (ip, infectious particles), as well as 
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particles that are inactive (nip, non-infectious particles) (Ruffing et al. 
(1994) J.Gen. Virol. 75:3385-3392; Kechli et al. (1998) Hum. Gene Ther. 
S:587-590.). The pp and the ip/nip ratio, are features of the packaging 
system, the manufacturing process and the vector itself. 
5 The infectious particles (ip) (infectious units, transducing units, 

etc.) are evidenced by the changes observed in the infected cells (vector 
DNA replication, provirus integration, cell lysis, transgene expression and 
other observable parameters). Infectious particles (ip) measures the 
number of particles effective in performing a process whose output is 

10 being measured; not all particles participate or are capable of participating 
in all processes. 

The precise assessment of ip is not straightforward. Existing 
methods are mainly based on serial dilution experiments followed by 
either linear extrapolation or asymptote approximates. The titer of 

15 infectious particles Up: infectious unity, transduction unity) (see, e.g., 
Mittereder et a/. (1 996) J. Virol. 70:7498-7509; Weitzman et al. (1996) 
J. Virol. 70:1 845-1 854; Salvetti et al. (1998) Hum. Gene Ther. 5:695- 
706) is evaluated by the studying observed changes in infected cells, 
such as viral replication, provirus integration, cellular lysis and transgene 

20 expression, using methods based on serial dilutions, followed either by a 
linear extrapolation or an asymptotic approximation. Thus, ip measures 
the number of active particles in the measured process; it includes 
physical particle (pp) and inactive particles {nip or noninfectious 
particles) . 

25 In order to resolve the problem of the titer determination and the 

comparison of different recombinant viruses used in gene therapy, the 
variation of the particles/infectious power ratio has been used (see, e.g., 
Atkinson et al. (1998) Nucl. Acids Res. 26:2821-2823; and International 
PCT application No. WO 99/1 1764, which describe a method that uses 
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step of amplification viral genetic material in a host cellular line, 
preparation of vectors of unknown titer obtained by serial dilution and an 
internal check of known titer). In particular, the method uses cells 
infected with a viral preparation in the different wells of a microtiter plate, 
5 viral genome replication in the host cells, nucleic acid hybridization and 
determination of the relative amount of replicated viral nucleic acid in 
each well. 

All of these methods measure the physical particles (pp) titer and/or 
measure infectious particles Up) titer in order to evaluate a gene transfer 

10 vector. A high quality vector preparation is one with an high titer of 
infectious particles and a low pplip ratio. These parameters provide 
quantitative information on the performance of a gene transfer vector. 
Because of the inaccuracy of the procedures used for assessing pp and 
especially ip, these parameters are not informative enough to precisely 

15 define the features of a gene therapy vector nor those of a particular 

preparation thereof. The actual procedures used for pp and ip evaluation 
change with the vector type, are not very reproducible nor exact, so 
these parameters do not contain enough information to allow a very fine 
definition of vector characteristics and performances. 

20 Hill equation-based analyses 

In this method complex biological processes, including those 
involving the response of cells (in vitro and in vivo) to biological agents, 
such as, for example, cells, viruses, vaccines, viral and non-viral gene 
vectors, antibodies, antigens, proteins in general and plasmids, are 

25 characterized using the formal analysis first introduced by Hill (see, Hill 
(1910) J. Physiol. 40:4P; Hill (1913) /. Biochem.J. 7:471-480). 
International PCT application No. WO 01/44809 (based on- 
PCT/FROO/03503, priority claimed to French application FR 9915884) 
describes the use of the Hill equation (see, Hill (1910) J. Physiol. 40:4?) 
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Hill (1913) /. Biochem.J. 7:471-480; see, International PCT application 
No. PCT/FROO/03503) for analysis and characterization of the biological 
and/or pharmacological activity of biological agents (viruses, vectors or 
cells) on biological assay systems in vitro (cell-based) or in vivo. 
5 A number of useful parameters, derived from the Hill equation, are 

scored and used to quantify relevant features of the biological agent, of 
the cells, as well as of the biological process or reaction involved. 

In particular, methods for calculation and analysis of the 
parameters of biological and pharmacological activity of native, 

10 attenuated, modified or recombinant viruses, vaccines, recombinant viral 
and non-viral gene transfer vectors, cells, antibodies and protein factors 
in in vitro (cell-based) or in vivo assays are described. This method is 
adapted for high throughput processes and is sufficiently accurate to 
allow a very fine definition of vector characteristics and performances. 

15 International PCT application No. WO 01/44809 provides, a 

standard process for evaluating the interaction between any biological 
agent, such as a gene therapy vector, with a complex biological system 
(living target cells). It provides a screening process for a pool of complex 
biological agents, in order to select test agents that have a desired 

20 property, activity, structure or whatever is being sought. 

Different biological agents and assay systems (cells) are compared 
and ranked out on the basis of their performance, assessed through the 
Hill parameters. Thus, the accurate analysis and comparison of the 
biological response of complex assay systems (in vitro and in vivo) to 

25 complex biological agents is achieved experimentally. The Hill-based 

analysis (rt,K,T,€ r n,0) is used for a variety of purposes, including, but not 
limited to: 
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i) validation and optimization of the manufacturing processes used 
to obtain the biological agents; 

ii) development and optimization of the components of the 
biological agents (proteins, genomes, genetic units); 

5 iii) development and optimization of assays and analytical tests for 

the characterization of the biological agents. 
The method includes the steps of: 

(a) preparation, for each biological agent, of a sample scale, 
obtained by a serial dilution of the biological agent at a R1 concentration, 
10 (b) incubation of each sample of the dilution scale obtained in (a), 

with the target cells at a constant concentration R2, 

(c) determination of the P product from the reaction R1 + R2, at a 
t moment, in each the sample; and 

(d) realization of a theoretical curve H from the experimental points 
15 R1 and P, for each biological agent by iterative approximation of 

parameters of the reaction R1 + R2 - P, at the t moment, in accordance 

with this equation: 

P = P ma x <T R1 > r + ("R1> r > r=1 n {2) 

in which: 

20 R1 represents the biological agent concentration in a sample 

from the scale; 

R2 is concentration of target cells [in vitro or in vivo) 

P (output) represents the product from the reaction R1 + 

R2 at a t moment; 
25 P max represents the reaction maximal capacity; 

k represents, at a constant R2 concentration, the resistance 

of the biological system for responding to the biological agent (resistance 

constant R2); 
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r represents a coefficient that depends on R1 and 
corresponds to the Hill coefficient; and 

77 represents the intrinsic power of the R1 biological agent to 
induce a response in the biological system (P production at the t 
5 moment), and 

(e) sorting the k and rr values obtained in (d) for each 
biological agent and the biological agent, and then ranking according to 
the values thereof. 

Using the parameters (n t K,T t e,€,n) the activity of a biological agent 
10 on a complex biological system, as well as its intrinsic features can be 
fully characterized and compared. In addition, different biological 
systems either in vitro cell-based) or in vivo can be compared. 
Hill Equation 
The Hill equation: 
15 r = n 

P = I P max . R1 r / (K + (R1) r ) R2 constant (1) 

r =1 

where R1, P, P max and K represent, respectively, the concentration 
of the reagent R1, the concentration of the product, the maximal capacity 

20 of the reaction and the 'affinity' constant between R1 and R2. The Hill's 
coefficient (r) is a function of R1 . The coefficient r is equal to 1 when 
independent non-interactive binding sites are involved between R1 and 
R2, such as in reactions that follow kinetics described by Michaelis- 
Menten; and r varies from 1 to n for systems where the sites involved in 

25 the interaction between the R1 and R2 are not independent from each 

other, and the affinity for R1 at any R2 binding site varies as a function of 
either i) the degree of occupancy of other R2 sites; ii) the concentration 
of R1 itself or iii) the concentration of other (positive or negative) 
regulators. 
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The Hill equation, thus, is a general formalization that describes the 
interaction reaction between molecules. It expresses the amount of 
product formed as a function of the concentration of the reagents and of 
the affinity constant of the system. Originally developed for the study of 
5 the dissociation between haemoglobin and oxygen, the Hill equation 
covers the formal Michaelis-Menten analysis of enzyme kinetics, the 
analysis of ligand-receptor binding and of the allosteric protein systems. 

According to Hill, for a simple reaction like 

R1 + R2-P 
10 K 

where the affinity between R1 and R2 changes with concentration 

of each, the Hill equation describes the accumulation of the product P as 

a function of the concentration of one of the reagents (R1) and of the 

intrinsic properties (K) of the system. 

15 This equation can be applied to complex biological systems. For 

example, the response of the cells to infection (P), can be analyzed by 
applying an Hill-type equation. The amount of cells growing in vitro (R2) 
are infected with increasing concentrations of recombinant viruses (R1), 
and (P) is monitored. A Hill equation is iteratively fitted to the 

20 experimental data. 

For analyses of viral output as exemplified herein, 

*- K r* 

virus + cell -* transduced cell -* output (viral genome replication), 
Equation (1) is specifically reformulated as: 

25 P = P max (/7R1) r / (k + (/7R1) r ) r = 1 n (2) 

where P, P max , R1, tt, r and k , as described above, represent, 
respectively, the output signal (P) (the level of viral gene expression, or 
the level of virus replication), the maximal output signal (Pmax), the initial 
concentration (R1) of infectious viral particles (those susceptible to trigger 

30 the process leading to P), the potency of the vector {rr; a factor that 


WO 03/023032 


PCT/IB02/03921 


-77- 

affects the concentration of the vector (R1) by its specific strength or 
activity, the Hill's coefficient (r) and the constant of resistance of the 
reaction or process (/c). 

The constant of resistance k 
5 The concept of k is analogous to those of dissociation, kinetics, 

equilibrium or affinity constants concepts for simple chemical and 
biological reactions, k is a feature of the process (reaction) and of the 
biological system tested (cell type), k is a key parameter for the 
characterization of the assay system and the assessment of its 

10 performance as a test for the reaction under study. 

k measures the internal resistance offered by the process or 
reaction triggered by the biological agent, to proceed to P. k is specific to 
a particular process or reaction tested. In addition k is specific to the 
particular biological system tested. Different cell lines and types will 

15 display different k for the same reaction. Moreover, factors affecting the 
performance of a cell to accomplish the reaction (like contaminants, toxic 
agents, etc.) affect k in that cell. 

Variations in k affect equation (2) by shifting the curve to the right 
or to the left, according to whether the value of k increases or decreases, 

20 respectively. All curves differing only in k are parallel each other. 

k finds its direct and practical application in i) assay development and 
validation and ii) assessment of the susceptibility or sensitivity of different 
cell types or tissues to undertake the reaction under study and to be 
affected by it. 

25 The potency n 

n measures the intrinsic potency of the biological agent to 
accomplish P against the resistance (k) offered by the reaction process. 
For every infectious virus particle (R1 ) added to the assay, the actual 
activity of the virus added is given by /rR1 . In order to report an output P, 
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the potency n has to push forward the reaction inside the cell against k. 
tt is specific to the particular biological agent for the reaction under study. 
tt is a feature of the biological agent. 

Different versions or variants of the biological agent will display 
5 different tt for the same reaction. Thus, mutations, conformational 
changes or other modifications on the biological agent are expected to 
change its tt for a given reaction process. 

The concept of tt is analogous to that of chemical activity by 
opposite to concentration for simple compounds, tt is a correction factor 
10 that affects the concentration (R1) of the biological agent to indicate its 
actual strength or activity for a given reaction process. 

Variations in tt affect equation (2) by shifting the curve to the right 
or to the left, according to whether the value of tt decreases or increases, 
respectively. Curves differing only in tt are not parallel each other. The 
15 slope of the curve given by equation (2) increases as tt increases. 

tt is a key parameter for the characterization of the biological agent 
and the assessment of its performance to accomplish the reaction under 
study, n finds its direct and practical application in i) biological agent 
optimization and development as it allows to compare the relative potency 
20 of variants of the agent. 

tt is a valuable tool in the field of vaccine, gene transfer vector and 
antibody development, for the comparison between two or more different 
agents or different versions of the same agent, for performance. Two 
agents, for instance, may elicit equivalent potencies for gene transfer, 
25 while their potencies for immunogenicity be different. The use of rr, a 
quantitative and accurate parameter for assessing potency, will allow for 
ranking the candidates according to their potency {i.e., for gene transfer, 
gene expression, immunogenicity and other such properties and activities) 
and to make rational decisions about the relative value of the agent leads. 
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The efficiency € 

" € measures the maximal global efficiency of the reaction process 
when a biological agent characterized by a given n value interacts with a 
biological system characterized by a given k. e is specific to the particular 
5 couple biological agent (n) I biological system (k) for the reaction under 
study. € is a feature of the global reaction process and intervening 
reagents. Changes in either n , k, or both, will lead to changes in e. 

The efficiency of the reaction process described by equation (2) is 
given by the increase in the output P that can be obtained by increasing 
10 the input Rl . Thus, the first derivative of P with respect to Rl, or the 

slope of the curve described by equation (2), gives the global efficiency of 
the reaction at every R1 input. The maximal global efficiency, or e, is 
given directly by either the slope at the inflection point of the curve 
described by equation (2) or by the maximum of its derivative <yP/JR1. 
15 The slope of the curve given by equation (2) and the maximum of rfP/£R1 
increase as € increases. 

€ is a key parameter for the characterization of the efficiency global 
process, considering the assay conditions and reagents all together. It is 
therefore useful for assay optimization once n and k have been fixed and 
20 to detect changes in n when k is kept constant or, inversely, changes in k 
while 77 is kept unchanged. 

The heterogeneity index rj 
n measures the internal heterogeneity of the reaction process under 
study. Complex processes include a huge chain of individual and causal 
25 events inside a multidimensional network of interrelated and interregulated 
biological reactions. Thus, the constant of resistance (k) for the particular 
reaction process under study is a macroscopic indicator of the global 
resistance of that process {k = a1/c1 +a2/c2 + ...an/cn/n). If the 
contribution of the individual microscopic constants of resistance (a 1*1, 
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a2*2, ...an/cn) for the individual steps involved in the process were 
homogeneous and no thresholds were present from one step to the next, 
then , no discontinuities in the increase of the Hill coefficient (i.e. in the 
change of k) with R1 should be observed. The existence of a major 
5 heterogeneity among the /ci values corresponding to the microscopic 
individual steps (i.e. the existence of thresholds for the intermediate 
steps) might lead to a macroscopic discontinuity in the system. 
Heterogeneity would cause a change in the rate of variation of the Hill 
coefficient and, which would require a jump in the macroscopic value of k 

10 in order for equation (2) to fit the data. 

The presence of internal heterogeneity in the reaction process can 
be detected by the appearance of steps in the rate of change of the Hill 
coefficient, corresponding to the Hill curve that fits the experimental data. 
/7_is defined as an index of heterogeneity and its value corresponds to the 

15 number of steps in the rate of variation of the Hill coefficient (one step, 
/7 = 1; two steps, /7 = 2; n steps, /7= n). 

q is a key parameter for the dissection and detailed analysis of the 
reaction process. It is useful for the independent optimization and 
development of every one of the steps identified by r?. 

20 As mentioned, the presence of steps in the rate of change of q 

translates in an abrupt discontinuity in k. Therefore, every step is 
determined by a different macroscopic constant of resistance *. Systems 
with q = 2, can thus be described by a Hill equation in which k takes two 
different values (/c1 and k2), according to the R1 interval considered. One 

25 part of the curve is described by *1 and the other by k2. 

Hill curves describing reaction processes characterized by /7 = 2, are 
hybrids generated from two parallel Hill curves differing only in /c. The 
transition from one curve to the other may alter the smooth change in the 
slope of the resulting Hill curve. 
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The apparent titer t 

In the Hill equation (2), when R1 increases, r increases from 1 to 
2,3,4... and P approaches its P max value. On the other direction, on the 
contrary, R1 can only decrease up to a minimal point (R1 min ), at which r 
5 and P reach their minimal values. The Hill sigmoidal curve is not 

symmetric, only the right arm is asymptotic (towards P max ). On the left 
arm, the curve has an origin at R1 min ; the empirical curve does not fit the 
data for values below R1 min . 

From a biological point of view, the fact that P does not exist for 
10 R1 below R1 min , means that there is no 'product' when the concentration 
of 'substrate' is lower than R1 min ; e.g. that the system is not responsive 
to concentrations below R1 mln . The minimal concentration of R1 that the 
system can detect and report is R1 mln . 

In terms of biological agents, R1 mln represents the minimal amount 
15 that can elicit a response in a given reporter system, and it is represented 
by t. The titer defined this way, is neither an asymptote value nor a value 
approached by extrapolation, but a precise parameter of the Hill equation, 
at the very mathematical origin of the curve. 

t measures the limiting dilution or apparent titer of the biological 
20 reagent. The value of r is determined by the limit of sensitivity of the 
biological assay system and of method used for the measurement of the 
product P; that is why it is said to be apparent titer. 

t is specific to the batch or stock of the biological reagent tested, r 
represents the apparent concentration of the biological agent and is 
25 expressed in units per volume, e.g. the maximal dilution of the biological 
agent that leads to the production of P. r is given by the maximal R1 for 
which the Hill coefficient reaches its minimal value (the Hill coefficient 
becomes constant at a value equal or close to 1). The concept of r 
corresponds to that of titer, of general use for viruses, antibodies and 
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vectors. Variations in r affect equation (2) by shifting the curve to the 
right or to the left, according to whether the value of r decreases or 
increases, respectively. 

r is a key parameter that measures the 'apparent' concentration of 
5 a stock of the biological agent, which is necessary for whatever use it will 
be given. 

The absolute titer 0 

9 is a the parameter that measures the absolute concentration 
(titer) of a stock or batch of the biological agent. The value of 8 is not 
10 determined by nor dependent on the limit of sensitivity of the biological 
assay system or of the method used for the measurement of the product 
P; that is why it is said to be absolute titer. 9 is specific to the batch or 
stock of the biological reagent tested. It represents the real physical 
concentration of the biological agent and is expressed in units per volume, 
15 e.g. the maximal dilution of the biological agent that leads to the 
production of P. 

6 is given by the following equation 
0/r=r/s (3),. 
where s is the sensitivity of the detection method. Therefore, for agents 
20 detected using the same method, the following expression is valid: 
91/7-1 /r1 = einllrl = 0n/7n/m = constant (4) 
Using equation (4), the ratio of the absolute titer 8, corresponding 
to two biological agent preparations, can be obtained from their 
respective rr and r. Variations in B affect the equation (2) by shifting the 
25 curve to the right or to the left and/or by changing its slope. 
Compensation between n and *r 

rr and k may appear to compensate to generate two different Hill 
curves (one differing in rr and the other one differing in k) that would 
apparently fit with the same experimental data. As rr and k have opposite 
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effects, two Hill curves; in which the increase in rr is compensated by the 
decrease in /c, and vice versa, may seem to represent the same curve, 
which could make it difficult to determine whether two Hill curves are 
different because a change in rr or in /c. 
5 Detailed analysis of the Hill curves indicates that rr and k do not 

compensate very well. Although curves differing in compensatory values 
of either rr or k may vary close each other, they do not fit exactly in any 
of the two regions of highest curvature (before and after the inflection 
point). This dispersion is caused by the fact that rr, but not k, changes 
10 the slope at the inflection point of the Hill curve. Therefore, e, which is 
the slope of the Hill curve at the inflection point, can be used to easily 
differentiate between two Hill curves that apparently compensate for rr 
and k. 

Conclusions 

15 The application of the Hill analysis to resolve complex biological 

processes is effective for the precise and objective understanding of 
processes like virus or vaccine action, entry, genome replication, 
transgene expression, vector/transgene immunogenicity, cytotoxicity and 
other such parameters. The analysis is independent of the virus vaccine, 

20 vector and protein type involved and from the output parameter and 

variable measured, such as the internalized vector DNA, transgene mRNA 
level and transgene product activity. 

As in the field of chemical pharmaceuticals, the structure of the 
potential drug (in this case the biological agent) must be optimized to a 

25 maximal possible intrinsic potency. In analytical development, the goal is 
to search for better performing reporter systems (the lowest possible k), 
as analytical tool. Two different systems characterized by constants kA 
and kB, respectively, can be compared (using the same biological agent) 
for their relative resistance or performance. 
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Complex systems involving the interaction of biological agents, 
such as viruses, vaccines, gene transfer vectors, antibodies proteins and 
living cells (either in vitro or in vivo) can be analyzed using the Hill 
equation. A complex succession of unitary processes, each of them 
5 susceptible to be individually analyzed by the Hill equation, as a global 
process, can be also described by the same equation as its constitutive 
steps. 

EXAMPLE 4 

Materials and Methods 
10 Cells: 

293 human embryo kidney (HEK) cells, obtained from ATCC, were 
cultured in Dulbecco's modified Eagle's medium containing 4.5 g/l 
glucose (DMEM; GIBGO-BRL) 10 % fetal bovine serum (FBS, Hyclone). 
Hela rep-cap 32 cells, described above, were obtained from Anna Salvetti 
15 (CHU, Nantes) and cultured in the medium described above. 
Plasmids: 

pNB-Adeno, which encodes the entire E2A and E4 regions and VA 
RNA I and II genes of Adenovirus type 5, was constructed by ligating into 
the polylinker of multiple cloning site of pBSII KS ( + /-) (Stratagene, San 

20 Diego, USA) the Sall-Hindlll fragment (9842-1 1555 nt) of Adenovirus 
type 5 and the BamHI-Clal fragment (21563- 35950) of pBR325. All 
fragments of adenovirus gene were obtained from the plasmid pBHG-10 
(Microbix, Ontario, Canada). pNB-AAV encodes the genes rep and cap of 
AAV-2 and was constructed by ligation of Xbal-Xbal PCR fragment 

25 containing the genome of AAV-2 from nucleotide 200 to 4480 into Xbal 
site of polylinker MCS of pBSIIKS( + /-). The PCR fragment was obtained 
from pAV1 (ATCC, USA). Plasmid pNB-AAV was derived from plasmid 
pVA1l, which contains the AAV genomic region, rep and cap. pNB-AAV 
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does not contain the AAV ITR's present in pAV1 . pAAV-CMV(nls)LacZ 
was provided by Dr Anna S^lvetti (CHU, Nantes). 

pCMV(nls)LacZ (rAAV vector plasmid) and pNB-Adeno were 
prepared in DH5a E. coli and purified by Nucleobond AX PC500 Kit 
5 (Macherey-Nagel), according to standard procedures. Plasmid pAAV- 
CMV(nls)LacZ is derived from plasmid psub201 by deleting the rep-cap 
region with SnaB I and replacing it with an expression cassette harboring 
the cytomegalovirus (CMV) immediate early promoter (407 bp), the 
nuclear localized £-ga!actosidase gene and the bovine growth hormone 
10 polyA signal (324 bp) (see, Chadeuf et al. (2000) J. Gene Med. 2:260- 
268. pAAV-CMV(nls)LacZ was provided by Dr Anna Salvetti. 
Virus: 

Wild type adenovirus (AV) type 5 stock, originally provided by Dr 
Philippe Moullier (CHU, Nantes), was produced accordingly to standard 
15 procedures. 

Construction of Rep mutant libraries 

25 pmol of each mutagenic primer was placed into a 96 PCR well 
plate. 15 //I of reaction mix (0.25 pmol of pNB-AAV), 25 pmol of the 
selection primer (changing one non-essential unique restriction site to a 

20 new restriction site), 2 //I of 10 X mutagenesis buffer (100 mM Tris- 
acetate pH7.5, 100 mM MgOAc and 500 mM KOAc pH7.5) was added 
into each well. The samples were incubated at 98°C for 5 minutes and 
then immediately incubated for 5 minutes on ice. Finally, the plate was 
placed at room temperature for 30 minutes. 

25 The primer extension and ligation reactions of the new strands 

were completed by adding to each sample: 7 pi of nucleotide mix (2.86 
mM each nucleotide and 1 .43 X mutagenesis buffer) and 3//I of a fresh 
1:10 enzyme dilution mix (0.025 U///I of native T7 DNA polymerase and 
1 U///I of T4 DNA ligase were diluted in 20 mM Tris HCI pH7.5, 10 mM 
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KCI, 10 mM 0- mercaptoethanol, 1 mM DTT, 0.1 mM EDTA and 50% 
glycerol). Samples were incubated at 37°C for 1 hour. The T4 DNA ligase 
was inactivated by incubating the reactions at 72°C for 15 minutes to 
prevent re-ligation of the digested strands during the digestion of the 
5 parental plasmid (pNB-AAV). 

Each mutagenesis reaction was digested with restriction enzyme to 
eliminate parental plasmids: 30 //I solution containing 3 //I of 10X enzyme 
digestion buffer and 10 units of restriction enzyme were added to each 
mutagenesis reaction and incubated at 37 °C for at least 3 hours. 
10 90 //I of the £. coli XLmutS competent cells (Stratagene, San Diego 

CA; supplemented with 1.5 /yl of £-mercaptoethanol to a final concen- 
tration of 25 mM) were aliquoted into prechilled deep-well plates. The 
plates were incubated on ice for 10 minutes and swirling gently every 2 
minutes. 

15 A fraction of the reactions that had been digested with restriction 

enzyme (1/10 of the total volume) was added to the deep well plates. The 
plates were swirled gently prior to incubation on ice for 30 minutes. A 
heat pulse was performed in a 42°C water bath for 45 seconds, the 
transformation mixture was incubated on ice for 2 minutes and 0.45 ml of 

20 preheated SOC medium (2% (w/v) tryptone, 0.5% (w/v) yeast extract, 
8.5 mM NaCI, 2.5 mM KCI, 10 mM MgCI 2 and 20 mM glucose at pH 7) 
was added. The plates were incubated at 37°C for 1 hour with shaking. 

To enrich for mutant plasmids, 1 ml of 2X YT broth medium (YT 
medium is 0.5% yeast extract, 0.5% NaCI, 0.8% bacto-tryptone), 

25 supplemented with 100 /yg/ml of ampicillin, was added to each 

transformation mixture and the cultures were grown overnight at 37 °C 
with shaking. Plasmid DNA isolation was performed from each mutant 
culture using standard procedure described in Nucleospin Multi-96 Plus 
Plasmid Kit (Macherey-Nagel). Five hundred jjg of the resulting isolated 
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DNA was digested with 1 0 units of the selection restriction enzyme in a 
total volume of 30//I containing 3 fj\ of 10X enzyme digestion buffer for 
overnight at 37 °C. 

A fraction of the digested reactions (1/10 of the total volume) were 
5 transformed into 40 p\ of Epicurian coli XL1 -Blue competent cells 

supplemented with 0.68 //I of jff-mercaptoethanol to a final concentration 
of 25 mM. After heat pulse, 0.45 ml of SOC was added and the 
transformation mixtures were incubated for 1 hour at 37 D C with shaking^ 
before to be plate on LB-ampicillin agar plates. The agar plates were 
10 incubated overnight at 37°C and the colonies obtained were picked up 
and grown overnight at 37 °C into deep-well plates. 

Four clones per reaction were screened for the presence of the 
mutation using restriction enzyme specific to the new restriction site 
introduced into the mutated plasmid with the selection primer. The cDNA 
15 from selected clones was also sequenced to confirm the presence of the 
expected mutation. 

Monitoring rAAV Production 

rAAV from each of the above wells, were produced by triple 
transfection on 293 HEK cells. 3 x 1 0 4 cells were seeded into each well 

20 of 96 micro-well plate and cultured for 24 hours before transfection. 
Transfection was made on cells at about 70% confluency. 25 kDa PEI 
(poly-ethylene-imine, Sigma-Aldrich) was used for the triple transfection 
step. Equimolar amounts of the three plasmids AV helper plasmid (pNB- 
Adeno), AAV helper plasmid (pNB-AAV or a mutant clone rep plasmid) 

25 and vector plasmid (pAAV-CMV(nls)LacZ) were mixed with 10 mM PEI by 
gently shaking. The mixture was the added to the medium culture on the 
cells. 60 hours after transfection, the culture medium was replaced with 
100 //I of lysis buffer {50 mM Hepes, pH 7.4; 150 mM NaCI; 1 m MgCI 2 ; 
1 mM CaCI 2 ; 0.01 % CHAPS). After one cycle of freeze-thawing the 
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cellular lysate was filtered through a millipore filter 96 well plate and 
stored at -80°C. 

rAAV infection particles (ip) 

Titers of rAAV vector particles were determined on HeLa rep/cap 
5 32 cells using standard dRA (serial dilution replication assay) test. Cells 
were plated 24 hours before infection at a density of 1 x 10 4 cells in 96- 
well plates. Serial dilutions of the rAAV preparation were made between 1 
and 1 x 10 6 //l and used for co-infection of the HeLa rep/cap 32 cells 
together with wt-AV type 5 (MOI 25). 48 hours after infection the ip 
10 were measured by real time PCR or by the quantification of biological 
activity of the transgene. 
Real Time PCR 

Infected HeLa rep/cap 32 cells were lysed with 50 //I of solution 
(50 mM Hepes, pH 7.4; 150 mM NaCI). After one cycle of freeze-thawing 

15 50 //I of Proteinase K (10 mg/ml) and the lysate were incubated one hour 
at 55 °C. The enzyme was inactivated by incubation 10 min at 96°C. 

For real time PCR, 0.2 //I of lysate was taken. Final volume of the 
reaction was 1 0 pi in 384 well plate using an Applied Biosystem Prism 
7900. The primers and fluorescence probe set corresponding to the CMV 

20 promoter were as follows: CMV 1 primer 5'- 

TGCCAAGTACGCCCCCTAT-3' (SEQ ID No. 733) (0.2 //M) and CMV 2 
primer 5 ' - A G GTC ATG T A CTG G G C AT A ATG C -3' (SEQ ID No. 734) (0.2 
//M) ; probe VIC-Tamra 5'-TCAATGACGGTAAATGGCCCGCCT-3' (SEQ 
ID No. 735) (0.1 //M). dRA plots were obtained by plotting the DNA copy 

25 number (obtained by real time PCR) vs. the dilution of the rAAV 
preparation. 
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/?-Galactosidase activity 

After 48 hours of infection, cells were treated with trypsin, and 
100 jjI of reaction solution (GalScreen Kit, Tropix) was added and 
incubated for one hour at 26 °C. Luminescence was measured in 
5 NorthStar (Tropix) HTS station. dRA plots were obtained plotting the 
intensity of £-Galatosidase activity vs. the dilution of the rAAV 
preparation. 

Mathematical Model for results analysis: 

Results were analyzed using the Hill equation-based analysis 

10 (designated NautScan™; see, Patent n° 9915884, 1999, France; 

published as International PCT application No. WO 01/44809 (PCT n° 
PCT/FR00/03503, Dec, 2000, see EXAMPLE above). Briefly, data were 
processed using a Hill equation-based model that allows extraction of key 
feature indicators of performance for each individual mutant. Mutants 

15 were ranked based on the values of their individual performance and 
those at the top of the ranking list were selected as Leads. 
Results 

Generation of diversity. 

To identify candidate amino acid (aa) positions on the rep protein 

20 involved in rep protein activity an Ala-scan was performed on the rep 
sequence. For this, each amino acid in the rep protein sequence was 
replaced with Alanine. To do this sets of rAAV that encode mutant rep 
proteins in which each differs from wild type by replacement of one 
amino acid with Ala, were generated. Each set of rAAV was individually 

25 introduced into cells in a well of a microtiter plate, under conditions for 
expression of the rep protein. The amount of virus that could be 
produced from each variant was measured as described below. Briefly, 
activity of Rep was assessed by determining the amount of AAV or rAAV 
produced using infection assays on HeLa Rep-cap 32 cells and by 
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measurement of AAV DNA replication using Real Time PCR, or by 
assessing transgene (jff-galactosidase) expression. The relative activity of 
each individual mutant compared to the native protein was assessed and 
"hits" identified. Hit positions are the positions in the mutant proteins 
5 that resulted in an alteration (selected to be at least about 20%), in this 
instance all resulted in a decrease, in the amount of virus produced 
compared to the activity of the native (wildtype) gene (see Fig. 3A). 

The hits were then used for identification of leads (see, Fig. 3B). 
Assays for Rep activity were performed as described for identification of 
10 the hit positions. Hit positions on Rep proteins and the effect of specific 
amino acids on the productivity of AAV-2 summarized in the following 
table: 


Hit position 

replacing amino acid (effect) 

4 (ttt) F 

(get) A (decrease) 


10 (aag) K 

(gcg) A (decrease) 


20 (ccc) P 

(gec) A (decrease) 


22 (att) 1 

(get) A (decrease) 


28 (tgg) W 

(gcg) A (decrease) 


32 (gag) E 

(gcg) A (decrease) 


38 (ccg) P 

(gcg) A (decrease) 


39 (cca) P 

(gca) A (decrease) 


54 (ctg) L 

(get) A (decrease) 


59 (ctg) L 

(gcg) A (decrease) 


64 (ctg) L 

(gcg) A (decrease) 


74 (ccg) P 

(gcg) A (decrease) 


86 (gag) E 

(gcg) A (decrease) 


! 88 (tac) Y 

(gec) A (decease) 


101 (aaa) K 

(gca) A (decrease) 
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Hit position 

replacing amino acid (effect) 

1 24 (ate) I 

(gee) A (decrease) 


125 (gag) E 

(gcg) A (decrease) 


127 (act) T 

(get) A (decrease) 


132 (ttc) F 

(gee) A (decrease) 


140 (ggc) G 

(gee) A (decrease) 


161 (acc) T 

(gee) A (decrease) ; 


1 63 (cct) P 

(get) A (decrease) 


175 (tat) Y 

(get) A (decrease) 


193 (ctg) L 

(gcg) A (decrease) j 


196 (gtg) V 

(gcg) A (decrease) 


197 (teg) S 

(gee) A (decrease) 


221 (tea) S 

(gca) A (decrease) 


228 (gtc) V 

(gcg) A (decrease) 

I 

231 (etc) L 

(gee) A (decrease) 

II 

234 (aag) K 

(gcg) A (decrease) 

I 

237 (acc) T 

(gee) A (decrease) 


250 (tac) Y 

(gee) A (decrease) 


258 (aac) N 

(gee) A (decrease) 


260 (egg) R 

(gcg) A (decrease) 


263 (ate) I 

(gec) A (decrease) 


264 (aag) K 

(gcg) A (decrease) 

— 

334 (ggg) G 

(gcg) A (decrease) 


335 (cct) V 

(get) A (decrease) 


337 (act) T 

(get) A (decrease) 


341 (acc) T 

(gee) A (decrease) 


342 (aac) N 

(gee) A (decrease) 
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Hit position 

replacing amino acid (effect) 


347 (ata) I 

(gca) A (decrease) 



350 (act) T 

(get) A (decrease 

(aat) N (increase) 


354 (tac) Y 

(gec) A (decrease) 



363 (aac) N 

(gec) A (decrease) 


5 

364 (ttt) F 

(get) A (decrease) 



367 (aac) N 

(gec) A (decrease) 



370 (gtc) V 

(gec) A (decrease) 



376 (tgg) W 

(gcg) A (decrease) 



381 (aag) K 

(gcg) A (decrease) 


10 

382 (atg) M 

(gcg) A (decrease) 



389 (teg) S 

(gcg) A (decrease) 



407 (tec) S 

(gec) A (decrease) 



41 1 (ata) I 

(gca) A (decrease) 



414 (act) T 

(get) A (decrease) 


15 

420 (tec) S 

(get) A (decrease) 



421 (aac) N 

(gec) A (decrease) 



422 (acc) T 

(gec) A (decrease) 



424 (atg) M 

(gcg) A (decrease) 



428 (att) I 

(get) A (decrease) 


20 

429 (gac) D 

(gec) A (decrease) 



438 (cag) Q 

(gcg) A (decrease) 



440 (ccg) P 

(gcg) A (decrease) 



451 (acc) T 

(gec) A (decrease) 



460 (aag) K 

(gcg) A (decrease) 


25 

1 462 (acc) T 

(gec) A (decrease) 

(ata) 1 (increase) 


484 (ttc) F 

(gec) A (decrease) 
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Hit position 

replacing amino acid (effect) 

488 (aag) K 

(gcg) A (decrease) 


495 (ccc) P 

(gec) A (decrease) 


497 (ccc) P 

(gecj A (decrease) 

(cga) R (increase) ! 

497 (ccc) P 

(gec) A (decrease) 

(etc) L (increase) 

497 (ccc) P 

(gec) A (decrease) 

(tac) Y (increase) 

498 (agt) S 

(get) A (decrease) 


499 (gac) D 

(gec) A (decrease) 


503 (agt) S 

(gcg) A (decrease) 


511 (tea) S 

(gca) A (decrease) 


512 (gtt) V 

(get) A (decrease) 


516 (teg) S 

(gcg) A (decrease) 


517 (acg) T 

(get) A (decrease) 

(aac) N (increase) 

518 (tea) S 

(gca) A (decrease) 


519 (gac) D 

(gcg) A (decrease) 


542 (ctg) L 

(gcg) A (decrease) 

(teg) S (increase) 


(gca) A (decrease) 

(age) S (increase) 

598 (gga) G 

(gca) A (decrease) 

(gac) D (increase) 

598 (gga) G 

(gca) A (decrease) 

(age) S (increase) 

600 (gtg) V 

(gcg) A (decrease) 

(ccg) P (increase) 

601 (cca) P 

(gca) A (decrease) 


Hit position 
(within intron) 

replacing sequence (effect) 

630 (tgc) 

gcg (decrease) 

cgc or tea or cct 
(increase) || 


25 The hits in other AAV serotypes (see, also FIGS. 5A and 5B) are as 

follows: 

HIT POSITION 
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428 

428 

428 

a no 
428 




429 

429 

429 

429 

429 

**-29 


438 

438 

438 

438 

A OQ 

H-OO 

4*34. 

440 

440 

440 

440 

a An 



451 

451 

451 

a e «i 

451 

A C 1 

as; 1 

4.4.7 

460 

460 

460 

460 

A CO 

4bO 

*fOU 

too 

462 

462 

462 

462 

462 

/f A O 

4.c»ff 

484 

484 

484 

484 

A O A 

484 

y| O A 

tou 

488 

488 

488 

488 

A O O 

4oo 

4QQ 

*frO*r 

495 

495 

495 

49 o 

A Q C 
4-9 O 


4Q1 

497 

497 

497 

49 / 

A Q"7 



498 

498 

498 

a no 

498 

*f9o 



499 

499 

499 

499 

499 

*l-99 


503 

. 503 

503 

503 

503 



511 

511 

511 

51 1 

51 1 

31 1 


512 

512 

512 

512 

CIO 

512 

CI o 
D U 


51 6 

516 

516 

516 

516 

516 

534 

517 

517 

517 

517 

517 

517 

535 

518 

518 

518 

518 

518 

518 

536 

519 

519 

519 

519 

519 

519 

537 

542 

543 

542 

542 

542 

543 

561 

548 

549 

548 

548 

548 

549 

567 

598 

599 

600 

600 

599 

599 


600 

602 

603 

603 

602 

602 

589 

601 

603 

604 

604 

603 

603 

590 


Sets of nucleic acids encoding the rep protein were generated. The 
rep proteins encoded by these sets of nucleic acid molecules were those 
in which each amino acid position identified as a "hit" in the ala-scan 
30 step, were each sequentially replaced by all remaining 1 8 amino acids 
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using site directed mutagenesis. Each mutant was designed, generated, 
processed and analyzed physically separated from the others in 
addressable arrays. No mixtures, pools, nor combinatorial processing 
were used. 

5 As in the first round (alanine scan), a library of mutant rAAV was 

generated in which each individual mutant was independently and 
individually generated in a independent reaction and such that each 
mutant contains only a single amino acid change and this for each amino 
acid residue. Again, each resulting mutant rep protein was then 

10 expressed and the amount of virus produced in cells assessed and 
compared to the native protein. 
Lead identification 

Since rep proteins that result in increased virus production are of 
interest, those mutants that lead to an increase in the amount of virus 

15 produced (2 to 10 times the native activity), were selected as "leads." 
Ten such mutants were identified. 

Based on the results obtained from the assays described above (i.e. 
titer of virus produced by each rep variant), each individual rep variant 
was assigned a specific activity. Those variant proteins displaying the 

20 highest titers were selected as leads (see Table above). Leads include: 
amino acid replacement of T by N at Hit position 350; T by I at Hit 
position 462; P by R at Hit position 497; P by L at Hit position 497; P by 
Y at Hit position 497; T by N at Hit position 517; G by S at Hit position 
598; G by D at Hit position 598; V by P at Hit position 600. 

25 Also provided are combinations of the above mutant Rep 78, 68, 

52 and 40 proteins, nucleic acids encoding the proteins, and recombinant 
AAV (any serotype) containing the mutation at the indicated position or 
corresponding position for serotypes other than AAV-2, including any set 
forth in the following table and corresponding SEQ ID Nos. Each amino 
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acid sequence is set forth in a separate sequence ID listing; for each 
mutation or combination thereof there is a single SEQ ID setting forth the 
unspliced nucleic acid sequence for Rep78/68, which for all mutations 
from amino acid 228 on, includes the corresponding Rep 52 and Rep 40 
5 encoding sequence as well. 

Amino acid sequences of exemplary mutant Rep proteins 



Seq no. 

gene 

position(s) 

codon(s) 


seq.l 

rep78 

4 

uL 1 


seq. 2 

rep68 

4 

ul 1 

10 

seq. 3 

rep78 

10 

GCG 


seq.4 

rep68 

10 

utu 


seq. 5 

rep78 

20 



seq. 6 

rep68 

20 

GCC 


seq. 7 

rep78 

22 

GCT 

15 

seq. 8 

rep68 

22 

GCT 


seq. 9 

rep78 

29 

GCG 


seq. 10 

rep68 

29 

GCG 


seq.1 1 

rep78 

38 

GCG 


seq. 12 

rep68 

38 

GCG 

20 

seq.1 3 

rep78 

39 

GCA 


seq. 14 

rep68 

39 



seq.1 5 

rep78 

53 

GCT 


seq. 16 

rep68 

53 

GCT 


seq. 17 

rep78 

59 

GCG 

25 

seq. 18 

rep68 

59 

GCG 


seq. 19 

rep78 

64 

GCT 


seq. 20 

rep68 

64 • 

GCT 


seq. 21 

rep78 

74 

GCG 


seq. 22 

rep68 

74 

GCG 

30 

seq. 23 

rep78 

86 

GCG 


seq. 24 

rep68 

86 

GCG 


seq. 25 

rep78 

88 

GCC 


seq. 26- 

rep68 

88 

GCC 


seq. 27 

rep78 

101 

GCA 

35 

seq. 28 

rep68 

101 

GCA 


seq. 29 

rep78 

124 

GCC 


seq. 30 

rep68 

124 

GCC 


seq. 31 

rep78 

125 

GCG 


seq. 32 

rep68 

125 

GCG 

40 

seq. 33 

rep78 

127 

GCT 


seq. 34 

rep68 

127 

GCT 


seq. 35 

rep78 

132 

GCC 


seq. 36 

rep68 

132 

GCC 


seq. 37 

rep78 

140 

GCC 

45 

seq. 38 

rep68 

140 

GCC 
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seq.39 

rep78 


seq.40 

rep68 


seq/41 

rep78 


seq.42 

rep68 

5 

seq.43 

rep78 


seq.44 

rep68 


seq.45 

rep78 


seq.46 

rep68 


seq.47 

rep78 

10 

seq.48 

rep68 


seq.49 

rep78 


seq.50 

rep68 


seq.51 

rep78 


seq.52 

rep68 

15 

seq.53 

rep78 


seq.54 

rep52 


seq.55 

rep68 


seq.56 

rep40 


seq.57 

rep78 

20 

seq.58 

rep52 


seq.59 

rep68 


seq.60 

rep40 


seq.61 

rep78 


seq.62 

rep52 

25 

seq.63 

rep68 


seq.64 

rep40 


seq.65 

rep78 


seq.66 

rep52 


seq.67 

rep68 

30 

seq.68 

rep40 


seq.69 

rep78 


seq.70 

rep52 


seq.71 

rep68 


seq.72 

rep40 

35 

seq.73 

rep78 


seq.74 

rep52 


seq.75 

rep68 


seq.76 

rep40 


seq.77 

rep78 

40 

seq.78 

rep52 


seq.79 

rep68 


seq.80 

rep40 


seq.81 

rep78 


seq.82 

rep52 

45 

seq.83 

rep68 


seq.84 

rep40 


seq.85 

rep78 


seq.86 

rep52 


seq.87 

rep68 

50 

seq.88 

rep40 


seq.89 

rep78 


-99- 


161 

GCC 

161 

GCC 

163 

GCT 

163 

GCT 

175 

GCT 

175 

GCT 

193 

GCG 

193 

GCG 

196 

GCC 

196 

GCC 

197 

GCC 

197 

GCC 

221 

GCA 

221 

GCA 

228 

GCG 

228 

GCG 

228 

GCG 

228 

GCG 

231 

GCC 

231 

GCC 

231 

GCC 

231 

GCC 

234 

GCG 

234 

GCG 

234 

GCG 

234 

GCG 

237 

GCC 

237 

GCC 

237 

GCC 

237 

GCC 

250 

GCC 

250 

GCC 

250 

GCC 

250 

GCC 

258 

GCC 

258 

GCC 

258 

GCC 

258 

GCC 

260 

GCG 

260 

GCG 

260 

GCG 

260 

GCG 

263 

GCC 

263 

GCC 

263 

GCC 

263 

GCC 

264 

GCG 

264 

GCG 

264 

GCG 

264 

GCG 

334 

GCG 
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seq.90 

rep52 

334 

GCG 


seq.91 

rep68 

334 

GCG 


seq.92 

rep40 

334 

GCG 


seq.93 

rep78 

335 

GCT 

5 

seq.94 

rep52 

335 

GCT 


seq.95 

rep68 

335 

GCT 


seq.96 

rep40 

335 

GCT 


seq.97 

rep78 

337 

GCT 


seq.98 

rep52 

337 

GCT 

10 

seq.99 

rep68 

337 

GCT 


seq.100 

rep40 

337 

GCT 


seq.101 

rep78 

341 

GCC 


seq.102 

rep52 

341 

GCC 


seq.103 

rep68 

341 

GCC 

15 

seq.104 

rep40 

341 

GCC 


seq.105 

rep78 

342 

GCC 


seq.106 

rep52 

342 

GCC 


seq.107 

rep68 

342 

GCC 


seq.108 

rep40 

342 

GCC 

20 

seq.109 

rep78 

347 

GCA 


seq.110 

rep52 

347 

GCA 


seq.111 

rep68 

347 

GCA 


seq.112 

rep40 

347 

GCA 


seq.113 

rep78 

350 

AAT 

25 

seq.114 

rep52 

350 

AAT 


seq.1 1 5 

rep68 

350 

AAT 


seq.116 

rep40 

350 

AAT 


seq.117 

rep78 

350 

GCT 


seq.118 

rep52 

350 

GCT 

30 

seq."M9 

rep68 

350 

GCT 


seq.120 

rep40 

350 

GCT 


seq.121 

rep78 

354 

GCC 


seq.122 

rep52 

354 

GCC 


seq.123 

rep68 

354 

GCC 

35 

seq.124 

rep40 

354 

GCC 


seq.125 

rep78 

363 

GCC 


seq.126 

rep52 

363 

GCC 


seq.127 

rep68 

363 

GCC 


seq.128 

rep40 

363 

GCC 

40 

seq.129 

rep78 

364 

GCT 


seq.130 

rep52 

364 

GCT 


seq.131 

rep68 

364 

GCT 


seq.132 

rep40 

364 

GCT 


seq.133 

rep78 

367 

GCC 

45 

seq.134 

rep52 

367 

GCC 


seq.135 

rep68 

367 

GCC 


seq.136 

rep40 

367 

GCC 


seq.137 

rep78 

370 

GCC 


seq.138 

rep52 

370 

GCC 

50 

seq.139 

rep68 

370 

GCC 


seq.140 

rep40 

370 

GCC 
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seq.141 rep78 376 GCG 

seq.142 rep52 376 GCG 

seq.143 rep68 376 GCG 

seq.144 rep40 376 GCG 

5 seq.145 rep78 381 GCG 

seq.146 rep52 381 GCG 

seq.147 rep68 381 GCG 

seq.148 rep40 381 GCG 

seq.149 rep78 382 GCG 

10 seq.150 rep52 382 GCG 

seq.151 rep68 382 GCG 

seq.152 rep40 382 GCG 

seq.153 rep78. 389 GCG 

seq.154 rep52 389 GCG 

15 seq.155 rep68 389 GCG 

seq.156 rep40 389 GCG 

seq.157 rep78 407 GCC 

seq.158 rep52 407 GCC 

seq.159 rep68 407 GCC 

20 seq.160 rep40 407 GCC 

seq.161 rep78 411 GCA 

seq.162 rep52 411 GCA 

seq.163 rep68 411 GCA 

seq.164 rep40 411 GCA 

25 seq.165 rep78 414 GCT 

seq.166 rep52 414 GCT 

seq.167 rep68 414 GCT 

seq.168 rep40 414 GCT 

seq.169 rep78 420 GCT 

30 seq.170 rep52 420 GCT 

seq.171 rep68 420 GCT 

seq.172 rep40 420 GCT 

seq.173 rep78 421 GCC 

seq.174 rep52 421 GCC 

35 seq.175 rep68 421 GCC 

seq.176 rep40 421 GCC 

seq.177 rep78 422 GCC 

seq.178 rep52 422 GCC 

seq.179 rep68 422 GCC 

40 seq.180 rep40 422 GCC 

seq.181 rep78 424 GCG 

seq.182 rep52 424 GCG 

seq.183 rep68 424 GCG 

seq.184 rep40 424 GCG 

45 seq.185 rep78 428 GCT 

seq.186 rep52 428 GCT 

seq.187 rep68 428 GCT 

seq.188 rep40 428 GCT 

seq.189 -rep78 429 GCC 

50 seq.190 rep52 429 GCC 

seq.191 rep68 429 GCC 
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seq.192 

rep40 


seq.193 

rep78 


seq.194 

rep 5 2 


seq.195 

rep68 

5 

seq.196 

rep40 


seq.197 

rep78 


seq.198 

rep 5 2 


seq.199 

rep68 


seq.200 

rep40 

10 

seq.201 

rep78 


seq.202 

rep52 


seq.203 

rep68 


seq.204 

rep40 


seq.205 

rep78 

15 

seq.206 

rep 5 2 


seq.207 

rep68 


seq.208 

rep40 


seq.209 

rep78 


seq.210 

rep52 

20 

seq.211 

rep68 


seq.212 

rep40 


seq.213 

rep78 


seq.214 

rep52 


seq.215 

rep68 

25 

seq.216 

rep40 


seq.217 

rep78 


seq.218 

rep52 


seq.219 

rep68 


seq.220 

rep40 

30 

seq.221 

rep78 


seq.222 

rep52 


seq.223 

rep68 


seq.224 

rep40 


seq.225 

rep78 

35 

seq.226 

rep52 


seq.227 

rep68 


seq.228 

rep40 


seq.229 

rep78 


seq.230 

rep52 

40 

seq.231 

rep68 


seq.232 

rep40 


seq.233 

rep78 


seq.234 

repS2 


seq.235 

rep68 

45 

seq.236 

rep40 


seq.237 

rep78 


seq.238 

rep 5 2 


seq.239 

rep68 


seq.240 

rep40 

50 

seq.241 

rep78 


seq.242 

rep52 


-102- 


429 

GCC 

438 

GCG 

438 

GCG 

438 

GCG 

438 

GCG 

440 

GCG 

440 

GCG 

440 

GCG 

440 

GCG 

451 

GCC 

451 

GCC 

451 

GCC 

451 

GCC 

460 

GCG 

460 

GCG 

460 

GCG 

460 

GCG 

462 

GCC 

462 

GCC 

462 

GCC 

462 

GCC 

462 

ATA 

462 

ATA 

462 

ATA 

462 

ATA 

484 

GCC 

484 

GCC 

484 

GCC 

484 

GCC 

488 

GCG 

488 

GCG 

488 

GCG 

488 

GCG 

495 

GCC 

495 

GCC 

495 

GCC 

495 

GCC 

497 

GCC 

497 

GCC 

497 

GCC 

497 

GCC 

497 

CGA 

497 

CGA 

497 

CGA 

497 

CGA 

497 

CTC 

497 

CTC 

497 

CTC 

497 

CTC 

497 

TAC 

497 

TAC 
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seq.243 

rep68 


seq.244 

rep40 


seq.245 

rep78 


seq.246 

rep52 

5 

seq.247 

rep68 


seq.248 

rep40 


seq.249 

rep78 


seq.250 

rep52 


seq.251 

rep68 

10 

seq.252 

rep40 


seq.253 

rep78 


seq.254 

rep52 


seq.255 

rep68 


seq.256 

rep40 

15 

seq.257 

rep78 


seq.258 

rep52 


seq.259 

rep68 


seq.260 

rep40 


seq.261 

rep78 

20 

seq.262 

rep 5 2 


seq.263 

rep68 


seq.264 

rep40 


seq.265 

rep78 


seq.266 

rep52 

25 

seq.267 

rep68 


seq.268 

rep40 


seq.269 

rep78 


seq.270 

rep 5 2 


seq.271 

rep68 

30 

seq.272 

rep40 


seq.273 

rep78 


seq.274 

rep 5 2 


seq.275 

rep68 


seq.276 

rep40 

35 

seq.277 

rep78 


seq.278 

rep52 


seq.279 

rep68 


seq.280 

rep40 


seq.281 

rep78 

40 

seq.282 

rep52 


seq.283 

rep68 


seq.284 

rep40 


seq.285 

rep78 


seq.286 

rep52 

45 

seq.287 

rep 6 8 


seq.288 

rep40 


seq.289 

rep78 


seq.290 

rep 5 2 


seq.291 

rep78 

50 

seq.292 

rep 5 2 


seq.293 

rep78 
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497 

TAC 

497 

TAC 

498 

GCT 

498 

GCT 

498 

GCT 

498 

GCT 

499 

GCC 

499 

GCC 

499 

GCC 

499 

GCC 

503 

GCG 

503 

GCG 

503 

GCG 

503 

GCG 

510 

GCA 

510 

GCA 

510 

GCA 

510 

GCA 

511 

GCA 

511 

GCA 

511 

GCA 

511 

GCA 

512 

GCT 

512 

GCT 

512 

GCT 

512 

GCT 

516 

GCG 

516 

GCG 

516 

GCG 

516 

GCG 

517 

GCT 

517 

GCT 

517 

GCT 

517 

GCT 

517 

AAC 

517 

AAC 

517 

AAC 

517 

AAC 

518 

GCA 

518 

GCA 

518 

GCA 

618 

GCA 

519 

GCG 

519 

GCG 

519 

GCG 

519 

GCG 

598 

GCA 

598 

GCA 

598 

GAC 

598 

GAC 

598 

AGC 
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5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


seq.294 

seq.295 

seq.296 

seq.297 

seq.298 

seq.299 

seq.300 

seq.301 

seq.302 

seq.303 

seq.304 

seq.305 

seq.306 

seq.307 

seq.308 

seq.309 

seq.310 

seq.311 

seq.312 

seq.313 

seq.314 

seq.315 

seq.316 

seq.317 

seq.318 

seq.319 

seq.320 

seq.321 

seq.322 

seq.323 

seq.324 

seq.325 

seq.326 

seq.327 

seq.328 

seq.329 

seq.330 

seq.331 

seq.332 

seq.333 

seq.334 

seq.335 

seq.336 

seq.337 

seq.338 

seq.339 

seq.340 

seq.341 

seq.342 

seq.343 

seq.344 


rep52 

rep78 

rep52 

rep78 

rep52 

rep78 

rep52 

rep78 

rep 5 2 

rep68 

rep40 

rep78 

rep68 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep 5 2 

rep68 

rep40 

rep78 

rep68 

rep78 

rep52 

rep68 

rep40 

rep78 

rep 5 2 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep 5 2 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 


598 
600 
600 
600 
600 
601 
601 

335 420 495 
335 420 495 
335 420 495 
335 420 495 
39 140 
39 140 
279 428 451 
279 428 451 
279 428 451 
279 428 451 
1 25 237 600 
1 25 237 600 
125 237 600 
1 25 237 600 
163 259 
163 259 
163 259 
163 259 
17 127 189 
17 127 189 
350 428 
350 428 
350 428 
350 428 
54 338 495 
54 338 495 
54 338 495 
54 338 495 
350 420 
350 420 
350 420 
350 420 
189 197 518 
189 197 518 
189 197 518 
189 197 518 
468 516 
468 516 
468 516 
468 516 

127 221 350 54 140 
127 221 350 54 140 
127 221 350 54 140 
127 221 350 54 140 


AGC 
GCG 
GCG 
CCG 
CCG 
GCA 
GCA 

GCT GCC GCC 
GCT GCC GCC 
GCT GCC GCC 
GCT GCC GCC 
GCA GCC 
GCA GCC 
GCC GCT GCC 
GCC GCT GCC 
GCC GCT GCC 
GCC GCT GCC 
GCG GCC GCG 
GCG GCC GCG 
GCG GCC GCG 
GCG GCC GCG 
GCT GCG 
GCT GCG 
GCT GCG 
GCT GCG 
GCG GCT GCG 
GCG GCT GCG 
GCT GCT 
GCT GCT 
GCT GCT 
GCT GCT 
GCC GCC GCC 
GCC GCC GCC 
GCC GCC GCC 
GCC GCC GCC 
GCT GCC 
GCT GCC 
GCT GCC 
GCT GCC 
GCG GCG GCA 
GCG GCG GCA 
GCG GCG GCA 
GCG GCG GCA 
GCC GCG 
GCC GCG 
GCC GCG 
GCC GCG 

GCT GCA GCT GCC GCC 
GCT GCA GCT GCC GCC 
GCT GCA GCT GCC GCC 
GCT GCA GCT GCC GCC 
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seq.345 
seq.346 
seq.347 
seq.348 
5 seq.349 
seq.350 
seq.351 
seq.352 
seq.353 

10 seq.354 
seq.355 
seq.356 
seq.357 
seq.358 

15 seq.359 
seq.360 
seq.361 
seq.362 
seq.363 

20 seq.364 
seq.365 
seq.366 
seq.367 
seq.368 

25 seq.369 
seq.370 
seq.371 
seq.372 
seq.373 

30 seq.374 
seq.375 
seq.376 
seq.377 
seq.378 

35 seq.379 
seq.380 
seq.381 
seq.382 
seq.383 

40 seq.384 
seq.385 
seq.386 
seq.387 
seq.388 

45 seq.389 
seq.390 
seq.391 
seq.392 
seq.393 

50 seq.394 
seq.395 


rep78 221 285 

rep52 221 285 

rep68 221 285 

rep40 221 285 

rep78 23 495 

rep52 23 495 

rep68 23 495 

rep40 23 495 

rep78 20 54 420 495 

rep52 20 54 420 495 

rep68 20 54 420 495 

rep40 20 54 420 495 

rep78 412 612 

rep52 412 612 

rep68 412 612 

rep40 412 612 

rep78 197 412 

rep52 197 412 

rep68 197 412 

rep40 197 412 

rep78 412 495 511 

rep52 412 495 511 

rep68 412 495 511 

rep40 412 495 511 

rep78 98 422 

rep52 98 422 

rep68 98 422 

rep40 98 422 

rep78 1 7 1 27 1 89 

rep68 17 127 189 

rep78 20 54 495 

rep52 20 54 495 

rep68 20 54 495 

rep40 20 54 495 

rep78 259 54 

rep52 259 54 

rep68 259 54 

rep40 259 54 

rep78 335 399 

rep52 335 399 

rep68 335 399 

rep40 335 399 

rep78 221 432 

rep52 221 432 

rep68 221 432 

rep40 221 432 

rep78 259 516 

rep52 259 516 

rep68 259 516 

rep40 259 516 

rep78 495 516 


GCA GCG 
GCA GCG 
GCA GCG 
GCA GCG 
GCT GCC 
GCT GCC 
GCT GCC 
GCT GCC 

GCC GCC GCC GCC 

GCC GCC GCC GCC 

GCC GCC GCC GCC 

GCC GCC GCC GCC 

GCC GCG 

GCC GCG 

GCC GCG 

GCC GCG 

GCG GCC 

GCG GCC 

GCG GCC 

GCG GCC 

GCC GCC GCA 

GCC GCC GCA 

GCC GCC GCA 

GCC GCC GCA 

GCC GCC 

GCC GCC 

GCC GCC 

GCC GCC 

GCG GCT GCG 

GCG GCT GCG 

GCC GCC GCC 

GCC GCC GCC 

GCC GCC GCC 

GCC GCC GCC 

GCG GCC 

GCG GCC 

GCG GCC 

GCG GCC 

GCT GCG 

GCT GCG 

GCT GCG 

GCT GCG 

GCA GCA 

GCA GCA 

GCA GCA 

GCA GCA 

GCG GCG 

GCG GCG 

GCG GCG 

GCG GCG 

GCC GCG 
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seq.396 rep52 495 516 GCC GCG 

seq.397 rep68 495 516 GCC GCG 

seq.398 rep40 495 516 GCC GCG 

seq.399 rep78 414 14 GCT GCC 

5 seq.400 rep52 414 14 GCT GCC 

seq.401 rep68 414 14 GCT GCC 

seq.402 rep40 414 14 GCT GCC 

seq.403 rep78 74 402 495 GCG GCC GCC 

seq.404 rep52 74 402 495 GCG GCC GCC 

10 seq.405 rep68 74 402 495 GCG GCC GCC 

seq.406 rep40 74 402 495 GCG GCC GCC 

seq.407 rep78 228 462 497 GCC GCC GCC 

seq.408 rep52 228 462 497 GCC GCC GCC 

seq.409 rep68 228 462 497 GCC GCC GCC 

15 seq.41 0 rep40 228 462 497 GCC GCC GCC 

seq.41 1 rep78 290 338 GCG GCC 

seq.41 2 rep52 290 338 GCG GCC 

seq.41 3 rep68 290 338 GCG GCC 

seq.41 4 rep40 290 338 GCG GCC 

20 seq.415 rep78 140 511 GCC GCA 

seq.416 rep52 140 511 GCC GCA 

seq.417 rep68 140 511 GCC GCA 

seq.418 rep40 140 511 GCC GCA 

seq.41 9 rep78 86 378 GCG GCG 

25 seq.420 rep52 86 378 GCG GCG 

seq.421 rep68 86 378 GCG GCG 

seq.422 rep40 86 378 GCG GCG 

seq.423 rep78 54 86 GCC GCG 

seq.424 rep68 54 86 GCC GCG 

30 seq.425 rep78 54 86 GCC GCG 

seq.426 rep68 54 86 GCC GCG 

seq.427 rep78 214 495 140 GCG GCC GCC 

seq.428 rep52 214 495 140 GCG GCC GCC 

seq.429 rep68 214 495 140 GCG GCC GCC 

35 seq.430 rep40 214 495 140 GCG GCC GCC 

seq.431 rep78 495 511 .GCC GCA 

seq.432 rep52 495 511 GCC GCA 

seq.433 rep68 495 51 1 GCC GCA 

seq.434 rep40 495 511 GCC GCA 

40 seq.435 rep78 495 54 GCC GCC 

seq.436 rep52 495 54 GCC GCC 

seq.437 rep68 495 54 GCC GCC 

seq.438 rep40 495 54 GCC GCC 

seq.439 rep78 197 495 GCG GCC 

45 seq.440 rep52 197 495 GCG GCC 

seq.441 rep68 197 495 GCG GCC 

seq.442 rep40 197 495 GCG GCC 

seq.443 rep78 261 20 GCC GCC 

seq.444 rep52 261 20 GCC GCC 

50 seq.445 rep68 261 20 GCC GCC 

seq.446 rep40 261 20 GCC GCC 
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seq.447 

rep78 

54 20 


seq.448 

rep68 

54 20 


seq.449 

rep78 

197 420 


seq.450 

rep 5 2 

197 420 

5 

seq.451 

rep68 

197 420 


seq.452 

rep40 

197 420 


seq.453 

rep78 

54 338 495 


seq.454 

rep52 

54 338 495 


seq.455 

rep68 

54 338 495 

10 

seq.456 

rep40 

54 338 495 


seq.457 

rep78 

197 427 


seq.458 

rep 5 2 

197 427 


seq.459 

rep68 

197 427 


seq.460 

rep40 

197 427 

15 

seq.461 

rep78 

54 228 370 387 


seq.462 

rep 5 2 

54 228 370 387 


seq.463 

rep68 

54 228 370 387 


seq.464 

rep40 

54 228 370 387 


seq.465 

rep78 

221 289 

20 

seq.466 

rep52 

221 289 


seq.467 

rep68 

221 289 


saq.468 

rep40 

221 289 


seq.469 

rep78 

54 163 


seq.470 

rep68 

54 163 

25 

seq.471 

rep78 

341 407 420 


seq.472 

rep 5 2 

341 407 420 


seq.473 

rep68 

341 407 420 


seq.474 

rep40 

341 407 420 


seq.475 

rep78 

54 228 

30 

seq.476 

rep52 

54 228 


sea. 477 

rep68 

54 228 


seq.478 

rep40 

54 228 


seq.479 

rep78 

96 125 511 


seq.480 

rep 5 2 

96 125 511 

35 

seq.481 

rep68 

96 125 511 


seq.482 

rep40 

96 125 511 


seq.483 

rep78 

54 1 63 


seq.484 

rep68 

54 163 


seq.485 

rep78 

197 420 

40 

seq.486 

rep52 

197 420 


seq.487 

rep68 

197 420 


seq.488 

rep40 

197 420 


seq.489 

rep78 

334 428 499 


seq.490 

rep52 

334 428 499 

45 

seq.491 

rep68 

334 428 499 


seq.492 

rep40 

334 428 499 


seq.493 

rep78 

197 414 


seq.494 

rep 5 2 

197 414 


seq.495 

rep68 

197 414 

50 

seq.496 

rep40 

197 414 


seq.497 

rep78 

30 54 127 


GCC GCC 
GCC GCC 
GCG GCC 
GCG GCC 
GCG GCC 
GCG GCC 
GCC GCC GCC 
GCC GCC GCC 
GCC GCC GCC 
GCC GCC GCC 
GCG GCG 
GCG GCG 
GCG GCG 
GCG GCG 

GCC GCC GCC GCG 
GCC GCC GCC GCG 
GCC GCC GCC GCG 
GCC GCC GCC GCG 
GCA GCC 
GCA GCC 
GCA GCC 
GCA GCC 
GCC GCT 
GCC GCT 
GCC GCC GCC 
GCC GCC GCC 
GCC GCC GCC 
GCC GCC GCC 
GCC GCC 
GCC GCC 
GCC GCC 
GCC GCC 
GCA GCG GCA 
GCA GCG GCA 
GCA GCG GCA 
GCA GCG GCA 
GCC GCT 
GCC GCT 
GCG GCC 
GCG GCC 
GCG GCC . 
GCG GCC 
GCG GCT GCC 
GCG GCT GCC 
GCG GCT GCC 
GCG GCT GCC 
GCG GCT 
GCG GCT 
GCG GCT 
GCG GCT 
GCG GCC GCT 
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seq.498 rep68 30 54 127 GCG GCC GCT 

seq.499 rep78 29 260 GCG GCG 

seq.500 rep52 29 260 GCG GCG 

seq.501 rep68 29 260 GCG GCG 

5 seq.502 rep40 29 260 GCG GCG 

seq.503 re P 78 4 484 GCT GCC 

seq.504 re P 52 4 484 GCT GCC 

seq.505 re P 68 4 484 GCT GCC 

seq.506 rep40 4 484 GCT GCC 

10 seq.507 rep78 258 124 132 GCC GCC GCC 

seS.508 rep52 258 124 132 GCC GCC GCC 

!S.509 rep68 258 124 132 GCC GCC GCC 

seq.510 rep40 258 124 132 GCC GCC GCC 

seq.51 1 rep78 231 497 GCC GCC 

15 seq.51 2 re P 52 231 497 6CCQCC 

seq.51 3 rep68 231 497 GCC GCC 

seq.514 rep40 231 497 GCC GCC 

seq.51 5 re P 78 221 258 GCA GCC 

seq.51 6 re P 52 221 258 GCA GCC 

20 seq.517 re P 68 221 258 GCA GCC 

seq.51 8 rep40 221 258 GCA GCC 

sea 519 rep78 234 264 326 GCG GCG GCC 

seq^O re P 52 234 264 326 GCG GCG GCC 

seq.521 re P 68 234 264 326 GCG GCG GCC 

25 seq.522 rep40 234 264 326 GCG GCG GCC 

seq.523 re P 78 153 398 AGC GCG 

seq.524 rep52 153 398 AGC GCG 

seq.525 rep68 153 398 AGC GCG 

seq.526 rep40 153 398 AGC GCG 

30 seq.527 rep78 53 216 GCG GCC 

seq.528 re P 68 53 216 GCG GCC 

seq.529 re P 78 22 382 GCT GCG 

seq.530 rep52 22 382 GCT GCG 

seq.531 rep68 22 382 GCT GCG 

35 seq.532 rep40 22 382 GCT GCG 

seq.533 rep78 231 41 1 GCC GCA 

seq.534 rep52 231 411 GCC GCA 

seq.535 rep68 231 411 GCC GCA 

seq.536 rep40 231 411 GCC GCA 

40 seq.537 rep78 59 305 GCG GCC 

seq.538 rep52 59 305 GCG GCC 

seq.539 rep68 59 305 GCG GCC 

seq.540 rep40 59 305 GCG GCC 

seq.541 rep78 53 231 GCG GCC 

45 seq.542 rep52 53 231 GCG GCC 

seq.543 rep68 53 231 GCG GCC 

seq.544 rep40 53 231 GCG GCC 

seq.545 rep78 258 498 GCC GCT 

seq.546 rep52 258 498 GCC GCT 

50 seq.547 rep68 258 498 GCC GCT 

seq.548 rep40 258 498 GCC GCT 
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seq.549 
seq.550 
seq.551 
seq.552 
5 seq.553 
seq.554 
seq.555 
seq.556 
seq.557 
10 seq.558 
seq.559 
seq.560 
seq.561 
seq.562 

1 5 DNA Sequences 

Sequence 

seq.563 

seq.564 

seq.565 
20 seq.566 

seq.567 

seq.568 

seq.569 

seq.570 
25 seq.571 

seq.572 

seq.573 

seq.574 

seq.575 
30 seq.576 

seq.577 

seq.578 

seq.579 

seq.580 
35 seq.581 

seq.582 

seq.583 

seq.584 

seq.585 
40 seq.586 

seq.587 

seq.588 

seq.589 


45 


50 


rep 7 8 

88 231 

rep 5 2 

88 231 

rep68 

88 231 

rep40 

88 231 

rep78 

101 363 

rep 5 2 

101 363 

rep68 

101 363 

rep40 

101 363 

rep78 

354 132 

rep 5 2 

354 132 

rep68 

354 132 

rep40 

354 132 

rep78 

10 132 

rep68 

10 132 


seq.590 


seq.591 


aa position 

codon 

4 

GCT 

10 

GCG 

20 

GCC 

22 

GCT 

29 

GCG 

38 

GCG 

39 

GCA 

53 

GCT 

59 

GCG 

64 

GCT 

74 

GCG 

86 

GCG 

88 

GCC 

101 

GCA 

124 

GCC 

125 

GCG 

127 

GCT 

132 

GCC 

140 

GCC 

161 

GCC 

163 

GCT 

175 

GCT 

193 

GCG 

196 

GCC 

197 

GCC 

221 

GCA 

228 (Rep78/68) 

GCG 

228 <Rep52) 

GCG 

228 (Rep 40) 

GCG 

231 (Rep78/68) 

GCC 

231 (Rep 52) 

GCC 

231 (Rep 40) 

GCC 

234 (Rep78/68) 

GCG 

234 (Rep 52) 

GCG 


GCC 
GCC 
GCC 
GCC 
GCA 
GCA 
GCA 
GCA 
GCC 
GCC 
GCC 
GCC 
GCG 
GCG 


GCC 
GCC 
GCC 
GCC 
GCC 
GCC 
GCC 
GCC 
GCC 
GCC 
GCC 
GCC 
GCC 
GCC 
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10 


15 


seq.592 
5 seq.593 
seq.594 
seq.595 
seq.596 
seq.597 
20 seq.598 
seq.599 
seq.600 
seq.601 
seq.602 
35 seq.603 
seq.604 
seq.605 
seq.606 
seq.607 
50 seq.608 


25 


30 


40 


45 
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234 (Rep 40) 

GCG 

237 (Rep78/68) 

GCC 

237 (Rep 52) 

GCC 

237 (Rep 40) 

GCC 

250 (Rep78/68) 

GCC 

250 

GCC 

250 

GCC 

258 (Rep78/68) 

GCC 

258 

GCC 

258 

GCC 

260 (Rep78/68) 

GCG 

260 

GCG 

260 

GCG 

263 (Rep78/68) 

GCC 

263 

GCC 

263 

GCC 

264 (Rep78/68) 

GCG 

264 

GCG 

264 

. GCG 

334 (Rep78/68) 

GCG 

334 

GCG 

334 

GCG 

335 (Rep78/68) 

GCT 

335 

GCT 

335 

GCT 

337 (Rep78/68) 

GCT 

337 

GCT 

337 

GCT 

341 (Rep78/68) 

GCC 

341 

GCC 

341 

GCC 

342 (Rep78/68) 

GCC 

342 

GCC 

342 

GCC 

347 (Rep78/68) 

GCA 

347 

GCA 

347 

GCA 

350 (Rep78/68) 

AAT 

350 

AAT 

350 

AAT 

350 (Rep78/68) 

GCT 

350 

GCT 

350 

GCT 

354 (Rep78/68) 

GCC 

354 

GCC 

354 

GCC 

363 (Rep78/68) 

GCC 

363 

GCC 

363 

GCC 

364 (Rep78/68) 

GCT 

364 

GCT 
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seq.609 
5 seq.610 
seq.61 1 

10 

seq.61 2 
seq.61 3 

15 

seq.614 
20 seq.61 5 
seq.61 6 

25 

seq.61 7 
seq.61 8 

30 

seq.61 9 
35 seq.620 
seq.621 

40 

seq.622 
seq.623 

45 

seq.624 
50 seq.625 
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364 

GCT 

367 (Rep78/68) 

GCC 

367 

GCC 

367 

GCC 

370 (Rep78/68) 

GCC 

370 

GCC 

370 

GCC 

376 (Rep78/68) 

GCG 

376 

GCG 

376 

GCG 

381 (Rep78/68) 

GCG 

381 

GCG 

381 

GCG 

382 (Rep78/68) 

GCG 

382 

GCG 

382 

GCG 

389 (Rep78/68) 

GCG 

389 

GCG 

389 

GCG 

407 (Rep78/68) 

GCC 

407 

GCC 

407 

GCC 

411 (Rep78/68) 

GCA 

411 

GCA 

411 

GCA 

414 (Rep78/68) 

GCT 

414 

GCT 

414 

GCT 

420 (Rep78/68) 

GCT 

420 

GCT 

420 

GCT 

421 (Rep78/68) 

GCC 

421 

GCC 

421 

GCC 

422 (Rep78/68) 

GCC 

422 

GCC 

422 

GCC 

424 (Rep78/68) 

GCG 

424 

GCG 

424 

GCG 

428 (Rep78/68) 

GCT 

428 

GCT 

428 

GCT 

429 (Rep78/68) 

GCC 

429 

GCC 

429 

GCC 

438 (Rep78/68) 

GCG 

438 

GCG 

438 

GCG 

440 (Rep78/68) 

GCG 

440 

GCG 
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440 GCG 

seq.626 451 <Rep78/68) GCC 

451 GCC 

451 GCC 

S seq.627 460 (Rep78/68> GCG 

460 GCG 

460 GCG 

seq.628 462 (Rep78/68) GCC 

462 GCC 

10 462 GCC 

seq.629 462 (Rep78/68) ATA 

462 ATA 

462 ATA 

seq.630 484 (Rep78/68) GCC 

15 484 GCC 

484 GCC 

seq.631 488 (Rep78/68) GCG 

488 GCG 

488 GCG 

20 seq.632 495 (Rep78/68) GCC 

495 GCC 

495 GCC 

seq.633 497 (Rep78/68) GCC 

497 GCC 

25 497 GCC 

seq.634 497 (Rep78/68) CGA 

497 CGA 

497 CGA 
seq.635 497 (Rep78/68) . CTC 

30 497 CTC 

497 CTC 

seq.636 497 (Rep78/68) TAC 

497 TAC 

497 TAC 
35 seq.637 498 (Rep78/68) w GCT 

498 GCT 

498 GCT 
seq.638 499 (Rep78/68) GCC 

499 GCC 
40 499 GCC 

seq.639 503 (Rep78/68) GCG 

503 GCG 

503 GCG 

seq.640 510 (Rep78/68) GCA 

45 510 GCA 

510 GCA 
seq.641 511 |Rep78/68) GCA 

511 GCA 

511 GCA 
50 seq.642 512 (Rep78/68) GCT 

512 GCT 
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seq.643 


5 seq.644 


seq.645 

10 

seq.646 


seq.647 

15 

seq.648 
seq.649 
seq.650 
20 seq.651 


seq.652 
seq.653 

25 

seq.654 


30 seq.655 


seq.656 
seq.657 

35 

seq.658 


40 seq.659 


seq.660 

45 

seq.661 


seq.662 

50 
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512 

GCT 

516 (Rep78/68) 

GCG 

516 

GCG 

516 

GCG 

517 {Rep78/68) 

GCT 

517 

w 1 / 

GCT 

517 

GCT 

517 (R6D78/68) 

J ] / \l * W# W W J 

AAC 

517 

AAC 

517 

AAC 

518 (Reo78/68) 

GCA 

518 

9 1 W 

GCA 

51 8 

GCA 

519 (ReD78/68) 

w 1 W \ 1 lO \J 1 w# ww/ 

GCG 

519 

GCG 

519 

GCG 

598 (Rep78/68) 

GCA 

600 (ReD78/68) 

GCG 

601 (ReD78/68) 

GCA 

335 420 495 

GCT GCC GCC 

335 420 495 

GCT GCC GCC 

335 420 495 

GCT GCC GCC 

39 140 

GCA GCC 

579 4.28 451 

*TfcW> *^W 1 

GCC GCT GCC 

979 428 451 

/ w ^*«W ^TW • 

GCC GCT GCC 

279 428 451 

GCC GCT GCC 

1 25 237 600 

1 W A.W » www 

GCG GCC GCG 

1 95 937 600 

I ^ w &w / www 

GCG GCC GCG 

1 95 237 600 

1 ^w AW / WWW 

GCG GCC GCG 

1 fi^ 95Q 

1 OJ w w 

GCT GCG 

163 259 

1 U w A ww 

GCT GCG 

163 259 

1 w O Aww 

GCT GCG 

17 127 189 

1 / \ 4- 1 1 WW 

GCG GCT GCG 

O ww tZO 

GCT GCT 


GCT GCT 


GCT GCT 

54 338 495 

GCC GCC GCC 

54 338 495 

GCC GCC GCC 

54 338 495 

GCC GCC GCC 

350 420 

GCT GCC 

350 420 

GCT GCC 

350 420 

GCT GCC 

189 197 518 

GCG GCG GCA 

189 197 518 

GCG GCG GCA 

189 197 518 

GCG GCG GCA 

468 516 

GCC GCG 

468 516 

GCC GCG 

468 516 

GCC GCG 

127 221 350 54 140 

GCT GCA GCT GCC GCC 

127 221 350 54 140 

GCT GCA GCT GCC GCC 

127 221 350 54 140 

GCT GCA GCT GCC GCC 
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seq.663 
seq.664 

5 

seq.665 
10 seq.666 
seq.667 

15 

seq.668 
seq.669 

20 

seq.670 
seq.671 

25 

seq.672 
seq.673 

30 seq.674 

seq.675 

35 

seq.676 
seq.677 

40 

seq.678 
45 seq.679 
seq.680 

50 

seq.681 
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221 285 

GCA GCG 

221 285 

GCA GCG 

221 285 

GCA GCG 

23 495 

GCT GCC 

23 495 

GCT GCC 

23 495 

GCT GCC 

20 54 420 495 

GCC GCC GCC GCC 

20 54 420 495 

GCC GCC GCC GCC 

20 54 420 495 

GCC GCC GCC GCC 

412 612 

GCC GCG 

412 612 

GCC GCG 

412 612 

GCC GCG 

197 412 

GCG GCC 

197 412 

GCG GCC 

197 412 

GCG GCC 

412 495 511 

GCC GCC GCA 

412 495 511 

GCC GCC GCA 

412 495 51 1 

GCC GCC GCA 

98 422 

GCC GCC 

98 422 

GCC GCC 

98 422 

GCC GCC 

17 127 189 

GCG GCT GCG 

20 54 495 

GCC GCC GCC 

20 54 495 

GCC GCC GCC 

20 54 495 

GCC GCC GCC 

54 163 

GCC GCT 

259 54 

GCG GCC 

259 54 

GCG GCC 

259 54 

GCG GCC 

335 399 

GCT GCG 

335 399 

GCT GCG 

335 399 

www www 

GCT GCG 

221 432 

GCA GCA 

221 432 

GCA GCA 

221 432 

GCA GCA 

259 516 

GCG GCG 

259 516 

GCG GCG 

259 516 

GCG GCG 

495 516 

GCC GCG 

495 516 

GCC GCG 

495 516 

GCC GCG 

414 14 

GCT GCC 

414 14 

GCT GCC 

414 14 

GCT GCC 

74 402 495 

GCG GCC GCC 

74 402 495 

GCG GCC GCC 

74 402 495 

GCG GCC GCC 

228 462 497 

GCC GCC GCC 

228 462 497 

GCC GCC GCC 

228 462 497 

GCC GCC GCC 

290 338 

GCG GCC 
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seq.682 

5 

seq.683 
seq.684 

10 

seq.685 
15 seq.686 
seq.687 

20 

seq.688 
seq.689 

25 

seq.690 
seq.691 

30 

seq.692 
seq.693 

35 

seq.694 

40 seq.695 

seq.696 
45 seq.697 

seq.698 

50 

seq.699 
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290 338 

GCG GCC 

290 338 

GCG GCC 

140 511 

GCC GCA 

140 511 

GCC GCA 

140 511 

GCC GCA 

86 378 

GCG GCG 

86 378 

GCG GCG 

86 378 

GCG GCG 

54 86 

GCC GCG 

54 86 

GCC GCG 

54 86 

GCC GCG 

214 495 140 

GCG GCC GCC 

214 495 140 

GCG GCC GCC 

214 495 140 

GCG GCC GCC 

495 511 

GCC GCA 

495 511 

GCC GCA 

495 511 

GCC GCA 

495 54 

GCC GCC 

495 54 

GCC GCC 

495 54 

GCC GCC 

197 495 

GCG GCC 

197 495 

GCG GCC 

197 495 

GCG GCC 

261 20 

GCC GCC 

261 20 

GCC GCC 

261 20 

GCC GCC 

54 20 

GCC GCC 

1 97 420 

GCG GCC 

1 97 420 

GCG GCC 

197 420 

GCG GCC 

54 338 495 

GCC GCC GCC 

54 338 495 

GCC GCC GCC 

54 338 495 

GCC GCC GCC 

197 427 

GCG GCG 

197 427 

GCG GCG 

197 427 

GCG GCG 

54 228 370 387 

GCC GCC GCC GCG 

54 228 370 387 

GCC GCC GCC GCG 

54 228 370 387 

GCC GCC GCC GCG 

221 289 

GCA GCC 

221 289 

GCA GCC 

221 289 

GCA GCC 

54 163 

GCC GCT 

54 163 

GCC GCT 

341 407 420 

GCC GCC GCC 

341 407 420 

GCC GCC GCC 

341 407 420 

GCC GCC GCC 

54 228 

GCC GCC 

54 228 

GCC GCC 

54 228 

GCC GCC 

96 125 511 

GCA GCG GCA 
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96 125 511 GCAGCGGCA 

96 1 25 51 1 GCA GCG GCA 

seq.700 197 420 GCG GCC 

197 420 GCG GCC 

5 197 420 GCG GCC 

seq.701 334 428 499 GCG GCT GCC 

334 428 499 GCG GCT GCC 

334 428 499 GCG GCT GCC 

seq.702 197 414 GCG GCT 

10 197 414 GCG GCT 

197 414 GCG GCT 

seq.703 30 54 127 GCG GCC GCT 

seq.704 29 260 GCG GCG 

29 260 GCG GCG 

15 29 260 GCG GCG 

seq.706 4 484 GCT GCC 

4 484 GCT GCC 

4 484 GCT GCC 

seq.707 258 1 24 1 32 GCC GCC GCC 

20 258 124 132 GCC GCC GCC 

258 1 24 1 32 GCC GCC GCC 

seq.708 231 497 GCC GCC 

231 497 GCC GCC 

231 497 GCC GCC 

25 seq.709 221 258 GCA GCC 

221 258 GCA GCC 

221 258 GCA GCC 

seq.710 234 264 326 GCG GCG GCC 

234 264 326 GCG GCG GCC 

30 234 264 326 GCG GCG GCC 

seq.71 1 1 53 398 AGC GCG 

153 398 AGC GCG 

153 398 AGC GCG 

seq.71 2 53 216 GCG GCC 

35 seq.71 3 22 382 GCT GCG 

22 382 GCT GCG 

22 382 GCT GCG 

seq.714 231 411 GCC GCA 

231 411 GCC GCA 

40 231 41 1 GCC GCA 

seq.71 5 59 305 GCG GCC 

59 305 GCG GCC 

59 305 GCG GCC 

seq.71 6 53 231 GCG GCC 

45 53 231 GCG GCC 

53 231 GCG GCC 

seq.71 7 258 498 GCC GCT 

258 498 GCC GCT . 

258 498 GCC GCT 

50 seq.71 8 88 231 GCC GCC 

88 231 GCC GCC 
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88 231 

GCC GCC 


seq.719 

101 363 

GCA GCC 


101 363 

GCA GCC 



101 363 

GCA GCC 

5 

seq.720 

354 132 

GCC GCC 


354 132 

GCC GCC 



354 132 

GCC GCC 


seq.726 

598 

GAC 


seq.727 

598 

AGC 

10 

seq.728 

600 

CCG 


The above nucleic acid molecules are provided in plasmids, which 
are introduced into cells to produce the encoded proteins. The analysis 
revealed the amino acid positions that affect Rep proteins activities. 
Changes of amino acids at any of the hit positions result in altered protein 

15 activity. Hit positions are numbered and referenced starting from amino 
acid 1 (nucleotide 321 in AAV-2 genome), also codon 1 of the protein 
Rep78 coding sequence under control of p5 promoter of AAV-2: 4, 20, 
22, 29, 32, 38, 39, 54, 59, 124, 125, 127, 132, 140, 161, 163, 193, 
196, 197, 221, 228, 231, 234, 258, 260, 263, 264, 334, 335, 337, 

20 342, 347, 350, 354, 363, 364, 367, 370, 376, 381, 389, 407, 41 1, 
414, 420, 421, 422, 424, 428, 438, 440, 451, 460, 462, 484, 488, 
495, 497, 498, 499, 503, 511, 512, 516, 517, 518, 542, 548, 598, 
600 and 601. The encoded Rep78, Rep68, Rep 52 and Rep 40 proteins 
and rAAV encoding the mutant proteins are provided. The corresponding 

25 nucleic acid molecules, Rep proteins, rAAV and cells containing the 
nucleic acid molecules or rAAV in which the native proteins are from 
other AAV serotypes, including, but are not limited to, AAV-1, AAV-3, 
AAV-3B, AAV-4, AAV-5 and AAV-6. 

Other hit positions identified include: 10, 64, 74, 86, 88, 101, 

30 175, 237, 250, 334, 429 and 519. 

Also provided are nucleic acid molecules, the rAAV, and the 
encoded proteins in which the native amino acid at each hit position is 
replaced with another amino acid, or is deleted, or contains additional 
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amino acids at or adjacent to or near the hit positions. In particular the 
following nucleic acid molecules and rAAV that encode proteins 
containing the following amino acid replacements or combinations 
thereof: T by N at Hit position 350; T by I at Hit position 462; P by R at 
5 Hit position 497; P by L at Hit position 497; P by Y at Hit position 497; T 
by N at Hit position 517; L by S at hit position 542; R by S at hit position 
548; G by D at Hit position 598; G by S at Hit position 598; V by P at Hit 
position 600; in order to increase Rep proteins activities in terms on AAV 
or rAAV productivity. The corresponding nucleic acid molecules, 
10 recombinant Rep proteins from the other serotypes and the resulting 
rAAV are also provided (see Figs. 5 and the above Table for the 
corresponding position in AAV-1, AAV-3, AAV-3B, AAV-4, AAV-5 and 
AAV-6). 

Mutant adeno-associate virus (AAV) Rep proteins and viruses 
15 encoding such proteins that include mutations at one or more of residues 
64, 74, 88, 175, 237, 250 and 429, where residue 1 corresponds to 
residue 1 of the Rep78 protein encoded by nucleotides 321-323 of the 
AAV-2 genome, and where the amino acids are replaced as follows: L by 
A at position 64; P by A at position 74; Y by A at position 88; Y by A at 
20 position 175; T by A at position 237; T by A at position 250; D by A at 
position 429 are provided. Nucleic acid molecules encoding these 
viruses and the mutant proteins are also provided. 

Also provided are nucleic acid molecules produced from any of the 
above-noted nucleic acid molecules by any directed evolution method, 
25 including, but are not limited to, re-synthesis, mutagenesis, recombination 
and gene shuffling and any way by combining any combination of the 

molecules, i.e., one, two by one, two by two, n by n, where n is the 

number of molecules to be combined ( i.e., combining all together). The 
resulting recombinant AAV and encoded proteins are also provided. 
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Also provided are nucleic acid molecule in which additional amino 
acids surrounding each hit, such as one, two, three ... ten or more, 
amino acids are systematically replaced, such that the resulting Rep 
protein(s) has increased or decreased activity. Increased activity as 
5 assessed by increased recombinant virus production in suitable cells is of 
particular interest for production of recombinant viruses for use, for 
example, in gene therapy. 

Also provided are combinations of the above noted mutants in 
which several of the noted amino acids are changed and optionally 
10 additional amino acids surrounding each hit, such as one, two, three . . • 
ten or more, are replaced. 

The nucleic acid molecules of SEQ ID Nos. 563-725 and the 
encoded proteins (SEQ ID Nos. 1-562 and 726-728) are also provided. 
Recombinant AAV and cells containing the encoding nucleic acids are 
15 provided, as are the AAV produced upon replication of the AAV in the 
cells. 

Methods of in vivo or in vitro production of AAV or rAAV using 
any of the above nucleic acid molecules or cells for intracellular 
expression of rep proteins or the rep gene mutants are provided. In vitro 

20 production is effected using cell free systems, expression or replication 
and/ or virus assembly. In vivo production is effected in mammalian cells 
that also contain any requisite cis acting elements required for packaging. 

Also provided are nucleic acid molecules and rAAV (any serotype) 
in which position 630 (or the corresponding position in another serotype; 

25 see Figs. 5 and the table above) has been changed. Changes at this 

position and the region around it lead to changes in the activity or in the 
quantities of the Rep or Cap proteins and/or the amount of AAV or rAAV 
produced in cells transduced with AAV encoding such mutants. Such 
mutations include tgc to gcg change (SEQ ID No. 721). Mutations at any 
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position surrounding the codon position 630 that increase or decrease the 
Rep or Cap proteins quantities or activities are also provided. Methods 
using the rAAV (any serotype) that contain nucleic acid molecules with a 
mutation at position 630 or within 1, 2, 3 ... .10 or more bases thereof 
5 for the intracellular expression rep proteins or the rep gene mutants 

covered by claims 1 0 to 13, for the production of AAV or rAAV (either in 
vitro, in vivo or ex vivo) are provided. In vitro methods include cell free 
systems, expression or replication and/or virus assembly. 

Also provided are rAAV (and other serotypes with corresponding 

10 changes) and nucleic acid molecules encoding an amino acid replacement 
by N at Hit position 350 of AAV- 1 , AAV-3, AAV-3B, AAV-4 and AAV-6 
or at Hit position 346 of AAV-5; by I at Hit position 462 of AAV-1 , 
AAV-3, AAV-3B, AAV-4 and AAV-6 or at Hit position 458 of AAV-5; by 
either R, L or Y at Hit position 497 of AAV-1, AAV-3, AAV-3B, AAV-4 

15 and AAV-6 or at Hit position 493 of AAV-5; by N at Hit position 517 of 
AAV-1, AAV-3, AAV-3B, AAV-4 and AAV-6 or at Hit position 535 of 
AAV-5; by S at hit position 543 of AAV-1 and AAV-6 or at hit position 
542 of AAV-3, AAV-3B and AAV-4 or at hit position 561 of AAV-5; by S 
at hit position 549 of AAV-1 and AAV-6 or at hit position 548 of AAV-3, 

20 AAV-3B and AAV-4 or at hit position 567 of AAV-5; by either D or S at 
Hit position 599 of AAV-1, AAV-4 and AAV-6 or at Hit position 600 of 
AAV-3 and AAV-3B; by P at Hit position 602 of AAV-1 , AAV-4 and AAV- 
6 or at hit position 603 of AAV-3 and AAV-3B or at hit position 589 of 
AAV-5 in order to increase Rep proteins activities as assessed by AAV or 

25 rAAV productivity. Methods using such AAV for expression of the 
encoded proteins and production of AAV are also provided. 
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Since modifications will be apparent to those of skill in this art, it is 
intended that this invention be limited only by the scope of the appended 
claims. 
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WHAT IS CLAIMED IS: 

1 . A process for the production of a peptide, polypeptide, or 
protein having a predetermined property, comprising: 

(a) producing a population of sets of nucleic acid molecules that 
5 encode modified forms of a target protein; 

(b) introducing each set of nucleic acid molecules into host cells 
and expressing the encoded protein, wherein the host cells are present in 
an addressable array; and 

(c) individually screening the sets of encoded proteins to identify 
10 one or more proteins that have activity that differs from the target 

protein, wherein each such protein is designated a hit, 

2. The process of claim 1, wherein each set of nucleic acid 
molecules is individually designed and synthesized. 

3. The process of claim 2, wherein each set is deposited at a 
15 locus in an addressable array. 

4. The process of claim 1 , wherein each polynucleotide in a set 
encodes a protein that differs by at least one amino acid from the target 
protein. 

5. The process of claim 1, wherein the array comprises a solid 
20 support with loci for containing or retaining cells; and each locus contains 

one set of cells. 

6. The process of claim 1 , wherein the array comprises a solid 
support with wells for containing or retaining cells; and each well contains 
one set of cells. 

25 7. The process of claim 1, wherein the nucleic acid molecules 

comprise viral vectors; and the cells are eukaryotic cells that are 
transduced with the vectors. 

8. The process of claim 1, wherein the nucleic acid molecules 
comprise plasmids and the cells are bacterial cells. 
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9. The method of claim 1 , further comprising: 

(d) modifying the nucleic acid molecules that encode the hits, to 
produce a set of nucleic acid molecules that encode modified hits; 

(e) introducing the each set of nucleic acids that encode the 
5 modified hits into cells; and 

(f) individually screening the sets of cells that contain the nucleic 
acid molecules that encode the modified hits to identify one or more ceils 
that encodes a protein that has activity that differs from the target protein 
and has properties that differ from the original hits, wherein each such 

10 protein is designated a lead. 

1 0. The process of claim 9, wherein each set nucleic acid 
molecules in step (d) is individually designed and synthesized. 

1 1 . The method of claim 1 , wherein the nucleic acid molecules in 
step (a) are produced by a method selected from among nucleic acid 

15 shuffling, recombination, site-directed or random mutagenesis and de 
novo synthesis. 

12. The method of claim 9, wherein the nucleic acid molecules in 
step (d) are produced by a method selected from among nucleic acid 
shuffling, recombination, site-directed or random mutagenesis, and de 

20 novo synthesis. 

13. The method of claim 2, wherein the nucleic acid molecules in 
step (a) are produced by systematically changing each codon in the target 
protein to a pre-selected codon. 

14. The method of claim 13, wherein the codon is selected from 
25 a codon encoding Ala (A), Ser (S), Pro (P) or Gly (G). 

15. The method of claim 13, wherein the codon is selected from 
a codon encoding Arg (R), Asn (N), Asp (D), Cys (C), Gin (Q), GIu (E), His 
(H), He <l), Leu (L), Lys (K), Met (M), Phe (F), Thr (T), Trp (W), Tyr (Y) or 
Val (V). 


WO 03/023032 


PCT/IB02/03921 


-124- 

16. The method of claim 9, wherein the nucleic acids of step (d) 
are produced by systematically replacing each codon that is a hit, with a 
codon encoding the remaining amino acids, to produce nucleic acid 
molecules each differing by at least one codon and encoding modified hits 

5 to identify leads. 

17. The method of claim 9, further comprising: 
recombining the nucleic acid molecules encoding the leads; 
introducing those nucleic acid molecules into cells; and 
screening the cells to identify nucleic acid molecules that encode 

10 optimized leads. 

18. The method of claim 17, wherein the recombining is two, 
three or more up to all of the nucleic acids encoding the leads. 

19. The method of claim 17, wherein the recombining is effected 
by a method selected from among nucleic acid shuffling, recombination, 

15 site-directed or random mutagenesis and de novo synthesis. 

20. The method of claim 1, wherein the modifications are 
effected in a selected domain of the target protein. 

21 . The method of claim 1 , wherein the modifications are 
effected alqng the full length of the target protein. 

20 22. The method of claim 1, wherein the change in activity is at 

least about 10%, 20%, 30%, 40% or 50%. 

23. The method of claim 1, wherein the change in activity is at 
least about 75%, 100%, 200%, 500% or 1000%. 

24. The method of claim 7, wherein at step (b) the titer of the 
25 viral vectors in each set of cells is assessed. 

25. The method of claim 24, wherein titering is effected by real 
time virus titering, comprising: 

(i) incubating the nucleic acid molecules or a vector 
(biological agent) comprising the nucleic acid molecules at an initial 
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concentration C, which is the unknown titer, with the host cells at a 
constant known concentration, D; 

(ii) measuring at successive times, an output signal, i; 

(Hi) determining the time X0, wherein: 
5 tf3 corresponds to i-/?; 

Pmxty < fi < Anax*' 

fi min and £ max correspond to values of i at the 
inflection point of the curve i = f(t), for the minimal and maximal values, 
respectively, of the concentrations of a reference biological agent for 
10 which the curve t£ = f(c) is predetermined; and 

(iv) determining the initial concentration C. 

26. The method of claim 24, wherein titering is effected by 
Tagged Replication and Expression Enhancement, comprising: 

(i) incubating with host cells a reporter virus vector 
15 with a titering virus of unknown titer, wherein the titering virus increases 

or decreases the output signal from the reporter virus; and 

(ii) measuring the output signal of the reporter virus 
and determining the titer of the reporter virus; 

(ii) determining the titer of the interfering virus by 
20 comparing the titer of the reporter virus in the presence and absence of 
the interfering virus. 

27. The process of claim 9, wherein the nucleic acid molecules 
comprise viral vectors; and the cells are eukaryotic cells that are 
transduced with the vectors. 

25 28. The method of claim 27, wherein at step (f) the titer of the 

viral vectors in each set of cells is determined. 

29. The method of claim 28, wherein the target protein is a 
protein involved in viral replication. 
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30. The method of claim 1, wherein the performance of the 
screened proteins is evaluated by a Hill analysis or by fitting the output 
signal to a curve representative of the interaction of the target protein and 
a test compound. 

5 31. The method of claim 30, wherein the Hill analysis, comprises: 

(a) preparing a sample of each nucleic acid molecule or a plasmid or 
vector that comprises each nucleic acid molecule (biological agent), 
wherein each sample is "obtained by a serial dilution of the molecules or 
vector or plasmid at a concentration R1, 
10 (b) incubating each sample of the dilution obtained in (a) with the 

host cells (target cells) at a constant concentration R2, 

(c) determining a P product from the reaction R1 + R2, at a t 
moment, in each the sample; and 

(d) preparing a theoretical curve H from the experimental points R1 
15 and P, for each biological agent by iterative approximation of parameters 

of the reaction R1 + R2 -* P, at the t moment, in accordance with the 
equation: 

P - Pmax (TR1) r / (k + (/rR1) r ) r=1„..,n (2) 

in which: 

20 R1 represents the biological agent concentration in a sample 

from the scale; 

R2 is concentration of target cells {in vitro or in vivo) 
P (output) represents the product from the reaction R1 + 
R2 at a t moment; 
25 P max represents the reaction maximal capacity; 

k represents, at a constant R2 concentration, the biological 
system for responding to the biological agent (resistance constant R2); 

r represents a dependent coefficient of R1 and corresponds 
to the Hill coefficient; and 
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n represents the intrinsic power of the Rl biological agent to 
induce a response in the biological system {P production at the t 
moment), and 

(e) sorting the k and n values obtained in (d) for each protein 
5 encoded by the nucleic acid molecules or plasmids or vectors and the 
cells, and then ranking according to the values thereof. 

32. The process of claim 1 that is automated. 

33. The process of claim 32 that is computer-controlled. 

34. A non-random method for generating proteins with a desired 
10 property, comprising: 

identifying a target protein; 

preparing sets of variant nucleic acid molecules that each encode a 
protein that differs by one amino acid from the target protein; 

screening and selecting the sets of variant nucleic molecules to 

15 identify those that encode proteins that have activity that differs by a 
predetermined amount from the activity of the target protein, thereby 
identifying proteins that are hits; 

identifying the residues in the hit proteins encoded by the variant 
nucleic acid molecules that differ from the target proteins; 

20 preparing further sets of variant nucleic acid molecules in which 

each codon in the nucleic acid molecule encoding each of the identified 
residues in each of the hits is replaced with codons encoding each of the 
remaining 1 8 amino acids to produce the further sets of variant nucleic 
acids, wherein in each set differs from each other set by one codon; and 

25 screening the further sets of nucleic acid molecules to identify 

those that encode proteins that have activity that differs from the activity 
of the hits, thereby identifying nucleic acid molecules that encode leads. 

35. The method of claim 34, wherein the replaced amino acid 
positions comprise a functional domain of the protein. 
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36. The method of claim 34, wherein the positions in the protein 
in which amino acids are replaced comprise at least about 50% of the 
amino acids in the protein. 

37. The method of claim 34, wherein the positions in the protein 
5 in which amino acids are replaced comprise at least 90% of the amino 

acids in the protein. 

38. The method of claim 34, wherein the positions in the protein 
in which amino acids are replaced comprise at least 95% of the amino 
acids in the protein. 

10 39. The method of claim 34, wherein the positions in the protein 

in which amino acids are replaced comprise all of the amino acids in the 
protein. 

40. The method of claim 34, wherein each set of nucleic acid 
molecules is generated, processed and screened separately or in parallel. 
15 41. A method for producing a protein having modified properties, 

comprising: 

(a) preparing a population of nucleic acid molecules that encode 
rationally modified proteins; 

(b) inserting the population into expression vectors; 

20 (c) introducing each vector into host cells therefor, and expressing 

the modified proteins, 

(d) screening each modified protein, and selecting one or more that 
has (have) a modified property. 
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10 20 30 40 50 60 

1 MPCIFTEIVIJCVT?SDLDEHLPGTSDSFVSOTA£ 60 

2 MPGl^EIVIitVPSDIXJEHI^GXSDSFVNWA^ 60 

3 MPGFYSIVXJCVPSDl^SHI^GXSNSFVNWArK^^ 60 

4 MPGFYEIVIJCVPSDLDEHI^GXSNSFVNWVAFK^ 60 

5 MPGFYEIVXJCVPSDLDEHI^GXSDSFVSP/VATKZVre^ 60 

6 MPGi^IVIJCVPSDLDGKIJKirSDSFVNOTASKEWE^ 60 

7 MATf?fSVIVRVPFDVEEHLPGXSDSFVDWVTGQIWELPPESDLNLTLVSQPQLrVADRIR 60 
C M**FYE** : *VP*D"*HLPGXS+SFV : WV* ** *WELPP*SD»* + *L*SQ W *LTVA* * * * 

70 BO 90 100 110 120 

1 RDFLVQWRRVSKAPEALFFVQFEKGJSTSyFHLKILVETTGViCSMVLGRFLSQIRDKLVQTI 120 

2 RD FL VQWRRVS KAPEALFFVQFEKGES YFHLH I LVETTGVJCS MVLGRF LS Q IRD KLVQT I 120 

3 RE FX VEWRRVS KAPEALF FVQFEKGETYFHLHVL I ET I GVJCS MWGRYV S Q I KS KLVTR I 120 

4 RE FL VE WRRVS KAP E ALF FVQ FE KGFTYFKLHVL I ET I GVJCS MWGRYV SQIKE KLVTR I 120 

5 REFL VEWRR VS KAPE AL F FVQ FE KGDS YFHLHI LVETVG VKS MWGR YVS Q I KS KL VTRI 120 

6 RDFLTEWRRVSKAPEALFFVQFEKGBSyFHMHVLVETTGVXSMVLGRFLSQIREKLIQRI 120 

7 RVFLYEWNKFSKQ-ESKFFVQFEKGSEyFKLHTLVETSGISSMVLGRYVSQIRAQLVKVV 119 
C R:FL++W***SK**E**FFVQFEKG+:YFH*K;L+ET:G**SMV:GR: :SQI: : *L* : :* 

130 140 150 160 170 180 

1 YRGXFPrLPNWFAVTKTRNGAGGGNKVVDECYI PNYLL PKTQPELQWAWTNMEE YISACL 180 

2 YRGXEPTLPNWFAVTKTRNGAGGGNKWDECY I PNYLL P KTQPE LQWAWTNMEE YI S ACL 180 

3 YRGVEPQLPNWFAVTKTRNGAGGGNKWDDCYI PNYLLPKTQPELQWAWTNMDQYLSACL 180 

4 YRGVFP Q L PNWFAVTKTRNGAGGGNKWDD C Y I PNYLL PKTQPSLQWAWTNMDQ YLS ACL 180 

5 YRGVEPQLPNWFAVTKTRNGAGGGNKWDDCYI PNYLLPKTQFELQWAWTNMDQ YI SACL 180 

6 YRGXFPrLPNWFAVTKTRNGAGGGNKVVDECYIPNYLLPKrQPELQWAWTNMEQ YLSACL 180 

7 FQGXFPQINDWVAITKVKK- - GGANKVVDSGYI PAYLLPKVQPSLQWAWTNLDE YKLAAL 177 
C **G:£P : * * *W*A*TK* * * * *GG*NKVVD : *YIP*YLLPK*QPELQWAWTN* : : Y:*A*L 

190 200 210 220 230 240 

1 NLAERKRLVAQHLTHVSQTQEQNKENLNPNSDAPVIRS KTSARYMELVGWL VDRG I TS EK 24 0 

2 NIJ^£RKRLVAHDLTHVSQTQEQNKE^LNPNSDAPVIRSKTSARYMELVGWiVI)RGirSEK 240 

3 NIAERKRLVAQHLTH VSQTQEQNKENQNPNSDAPVIRSKTSARYMELVGWLVDRG I TS EK 2 4 0 

4 NLAERKRLVAQHLTHVSQTQEQNKENQNPNSDAPVIRSKT5ARYMELTOWLVDRGirSEK 24 0 

5 NLAERKRLVAQKLTHVSQTQEQNKENQNPNSDAPVIRS KTS AR YME L VGWL VDRG I TS E K 24 0 

6 NLTERKRLVAQHLTHVSQTQEQNKENQNPNSDAPVIRSKTSARYMELVGWLVDiCGirSEK 24 0 

7 NLEERKRLVAQFLAESSQRS - QEAASQRSFSADPVIKSKTSQKYMALVNWLVEHGITSEK 236 
C NL+ERKRLVA*+L***SQ***Q**** + ***S**PVI*SKTS**YM*LV*WLV*+GirSEK 

250 260 270 280 290 300 

1 QWIQEDQASYISFNAASNSJ^SQXiCAALDNAGKIMALTKSAPDYLVGPAPPADIKTNRI YR 300 

2 QW I QEDQAS YI S FNAASitfSLRSQXJCAALDNAGKIMALTKS APD YLVG PAP PAD I KTNR. I YR 300 

3 QWIQEDQASYTSFNAASJ/SRSQXiCAALDNASKIMSLTKTAPDYLVGSNPPEDITKNRIYQ 3 00 

4 QWI QEDQAS YISFNAASJ/SRSQXXAALDNASKIMSLTKTAPDYLVGSNPPEDITKNRIYQ 300 

5 QWI QEDQAS YT S FNAASNSRS QXJCAALDNAS KI MS LTKTAPD YLVGQNP PED I S SNR I YR 3 00 

6 QW I QEDQAS YIS FNAASN'S.RSQXiCAALDNAGKIMS LTKTAPD YLVGQQFVED I S SNR I YK 300 

7 QWIQENQESYLSFNSTG^SJ^SQXiCAALDNATKIMSLTKSAVDYLVGSSVPEDISKNRIWQ 2 96 
C QWIQE*Q*SY*SFN***2tfSRSQXJCAALDNA: KIM+LTK: A*DYLVG : :**+DI: :NRI* : 

310 320 330 340 350 360 

1 ILELNGYEPAYAGSVFLGWAQK3*FGKRNTIWLFGPATTGKTNIAEAXAHAVPFYGCVN^ 360 

2 I LELNGYDPAYAGS VFLGWAQKRFGKRNTI WLFGPArTGKTtfl AEAXAHAVPFYGCVNWT 360 

3 I LELNGYDPQYAASVFLGWAQKKFGKRNTIWLFGPAITGKTNI AEAXAHAVPFYGCVNWT 360 

4 I LELNGYDPQYAASVFLGWAQKKFGKRNTIWLFGPATTGKrtfl AEAXAHAVPFYGCVNWT 360 

5 ILEMNGYDPQYAASVFLGWAQKKFGKRNTIWLFGPArTGKTNIAEAXAHAVPFYGCVNWT 3 60 

6 I LELNGYDPQYAAS VFLGWATKKFGKRNT I WLFGFA2TGKTJ7I AEAXAHTVP F YGCVNWT 3 60 
-7 I FEMNGYDPAYAGS ILYGWCQRSFNKRNTVWLYGPAXTGICrNI AEAXAHTVPF YGCVNWT 356 

C I*E+NGY»P: YA:S»**GW*** : F*KRNT»WL*G PATTGKTOTAEAXAH + VP F YGCVNWT 
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370 380 390 400 410 420 

1 NEtfFPFNDCVTDKMVIFWEEGiO^ 420 

2 NENFPFOTCVTDKMVIWWEEGiOrcA^ 420 

3 NENFPFlTOCVDKMVIWWEEGi^AK^ 420 

4 NENFPFJroCVDKKVIWWSEGJafTAKVVESAKAILGGS 420 

5 NEJ^PFNDCVDKMVIWWEEGJCMTAKVVESAKAILGGSKV^ 420 

6 NEZ^PFi^CVDKMVIWV^EGiCMTAKVVESAXAIL^ 420 

7 NE^FPFZTDCVDKMLIPWEEGiOifrNKVVESAKAl LGGSKVRVDQKCKSSVQXDSTPVIVTS 416 
C NEMFPFNDCVDKM* IWWEEGKMT* KWES AKAI LGGS KVRVDQKC KS S * QX + *TPVIVTS 

430 440 450 460 470 480 

1 NT^^MCAVXDGNSTTFEHOQPLQDRMFKFEL^RRLEHDFGiCV^KQEVKEFFRWAQDHVTEV 480 

2 imWCAVXDGNSTTFEHQQPLQDRMFKTCLrRRLEHDFG^ 480 

3 imCtfCAVXDGNSTTFEHQQPLQDRMFKFELrRRI^HDFGiCVrKQEVKDFFRWASDHVT 4 80 

4 JTOBtfCAVXDGNSTTFEHQQPLQDRMFKFELrRRI^^ 4 S 0 

5 /miMCAVXDGNSTTFEHOQPLQDRMFKFELrKRLEHDFGiCVrKQEVKDFFRWASDHVTEV 4 80 

6 imCtfCAVXDGNSTTFEHCQPLQDRMFKFEL 4 80 . 

7 NTNWCVVVDGNSTTFEHOQPLEDRMFKFELrKRLPPDFGJCirKQEVKDFFAWAKVNQVPV 476 

C NTNMC*V*PGNSTTFEHOQPL*DRMFKFELr+RL : *DFGJC*TKQEVK+FF*WA :***+: V 

490 500 510 520 

1 AHEFYVRJtGGANKRPAPDDADKSEPKRA CPSVADPSTSDAEG 522 

2 AHEJFYVRJCGGANKRPAPDDADKSE PKRA - CPSVADPSTSDAEG 522 

3 AHELFYVRJCGGAKKRPASNDAD VSE PKRQ - CTSLAQPTTSDAEA 522 

4 AKRFYVRXGGAKKRJPASND AD VSE P KRQ CTSLAQ PTTSD AEA 522 

5 THEfYVRJCGGARKRP APND AD ISE P KRA - CPSVAQPSTSDAEA 522 

6 EHEFYVKJCGGAKKRPAPSDADISEPKRV.- - - - -RESVAQPSTSDAEA 522 

7 THEFKVPRELAGTKGAEKSLKRPLGDVTNTSYKSLSKRARLSFVPETPRSSDVTVDPAPL 536 
C :HEF ir V+***A: ***A: :***.***** : +:*:*:*:** *A* : 

530 540 550 560 570 580 

1- APVDFADRYQNKCSRHAGMLQMLFPCKTCERMNQNFNICFTHGTRDCSECFP- -GVSESQ 580 

2 APVDFADRYQNKCSRHAGMLQMLFPCKTCERMNQNFNICFTHGTRDCSECFP - - GVSESQ 580 

3 p-ADYADRYQNKCSRHVGMNLMLFPCXTCERMNQISNVCFTHGQRDCGECFPGMSESQPV 581 

4 P-ADYADRYQNKCSRHVGMNLMLFPCKTCERMNQISNVCFTHGQRDCGECFPGMSESQPV 581 

5 P - VDYADRYQNKCSRKVGNttlLMLFPCRQCERMNQNVDICFTHGVMDCAECFP- VSESQPV 5 80 

6 S - INYADRYQNKCSRHVGKNLMLFPCRQCERMNQNSNICFTKGQKDCLECFP - - VSESQP 57 9 

7 RPLNWNSRYDCKCDY^QFDNISNKCDECEYLNRGKNGCICHNVTKCQICHG- - 588 

C : : : + : **RY* W KC**K : ** : : *» W *C: :CE**N* : : * :C wtr H* : : *C. *C** . . : : + : : 

590 600 610 620 

1 PVVRKRTYRKLCAIKHLLGRAPEIACSACDLVNVDLDDCVSEQ 623 

2 P WRKRTYRKLCAIKHLLGRAPE I ACS ACDLVNVDLDDCVS EQ 623 

3 S WKKKTYQKLC P IKKI LGRAPE I ACS ACDLANVDLDDC VSEQ 624 

4 SVVKKKTYQKLCPIHHILGRAPEIACSACDIANVDLDDCVSEO $ 24 

5 SVVRKRTYQKLCPIHHIMGRAPEVACSACELANVDLDDCDMEQ 623 

6 VSWKKAYQKLCYIHHIMG- KVPDACTACDLVNVDLDDCIFEQ €21 

7 IPPWEKENLSDFGDFDDANKEQ 610 

C : +*:*:+*:***:***::*++++:**+**+*:* *D*DD* :: EQ 
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